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Spinor Bose gases form a family of quantum fluids manifesting both magnetic order and superflu- 
idity. This article reviews experimental and theoretical progress in understanding the static and 
dynamic properties of these fluids. The connection between system properties and the rotational 
symmetry properties of the atomic states and their interactions are investigated. Following a 
review of the experimental techniques used for characterizing spinor gases, their mean-field and 
many-body ground states, both in isolation and under the application of symmetry-breaking ex- 
ternal fields, are discussed. These states serve as the starting point for understanding low-energy 
dynamics, spin textures and topological defects, effects of magnetic dipole interactions, and vari- 
ous non-equilibrium collective spin-mixing phenomena. The paper aims to form connections and 
establish coherence among the vast range of works on spinor Bose gases, so as to point to open 
questions and future research opportunities. 
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I. INTRODUCTION 



Nature has provided us with few quantum fluids, in 
which macroscopic characteristics of the fluid derive di- 
rectly from quantum coherences. Among such fluids, 
ones with a non-trivial internal degree of freedom are 
even rarer, examples being the d-wave and p-wave su- 
perconductors and certain phases of superfluid 3 He. Ul- 
tracold atomic physics has supplied us with a new family 
of such fluids, degenerate Bose gases with a spin degree 
of freedom. These so-called "spinor gases" are subject 
to the interplay of magnetism and superfluidity, both 
of which involve quantum phase coherence, long-range 
order and symmetry breaking. These fluids are inter- 
esting in their own right, as a newly conceived material 
now available for experimental investigation. In addition, 
since these atomic gases may be described starting from 
a simple theoretical framework, and since their proper- 
ties may be readily manipulated and measured in exper- 
iments, their study promises to shed insight on a range 
of topics such as the role of symmetry and topology in 
quantum-ordered materials, quantum phase transitions, 
non-equilibrium quantum dynamics, and the entangle- 
ment and squeezing of quantum fields. 

The study of degenerate spinor Bose gases w as sparked 
in 1998 by experiments on ultracold rubidium flHall et al 
1998a|bl |Matthews et al.\ |1998[ ) and sodium ( |Stenger 
et al.\ 1998), followed quickly by theoretical investiga- 



freedom, so that the state description of a single molecule 
or atom or of a system assembled of many such particles 
must describe both the external degrees of freedom, i.e. 
the center-of-mass motional state, and also the internal 
degrees of freedom. Quantum fluids composed of such 
particles may be described by a multi-component order 
parameter. In comparison, liquid 4 He in its superfluid 
state is described by a scalar order parameter. The in- 
terplay between the external and internal degrees of free- 
dom in the multi-component systems leads to a range of 
phenomena unfamiliar from studies of scalar quantum 
fluids. 

At first glance, it would seem that one has an enor- 
mous selection of multi-component quantum gases from 
which to choose. In fact, this plethora is severely depleted 
by the instability of most atomic and molecular internal 
states. Electronic excited states will typically decay due 
to spontaneous emission on a time scale much faster than 
that required for a gas to equilibrate kinetically through 
collisions. Short lifetimes are also expected (or observed) 
in many instances even within the electronic ground-state 
family of states, due to inelastic collisions that release 
large amounts of energy. 

However, in several cases, one does find one or more 
sets of internal states that may coexist for sufficiently 
long times to allow experimental studies of equilibrated 
multi-component gases. Let us discuss three examples: 



, 1. Hyperfine spin manifolds of hydrogen- 1 ike atoms 



tions of Bose-Einstein condensation in an interacting 
spin-1 gas pol 119981 lOhmi and Machidal 119981). A sig- 



nificant opus of literature has amassed since those works 
covering an extremely wide range of topics. The present 
review is intended to coalesce the understandings pre- 
sented in that work and to focus attention on core physi- 
cal concepts brought up in the study of spinor gases. We 
hope such a discussion will broaden the impact of spinor- 
gas research and draw new participants and perspective 
to this topic. 



A. Internal degrees of freedom of Bose gases 

By now, a wide selection of atoms and a few simple 
molecules have been produced in gaseous form in the ul- 
tracold regime 1 . Such objects carry internal degrees of 



1 The accepted definition of "ultracold" is time varying. Spurred 
by the quest for correlated atomic systems with ever-fewer de- 
fects and ever-subtler types of correlation, physicists are devising 
new methods to reach ever-lower temperature regimes. In this 
Review we will define relevant temperature regimes more specif- 
ically with respect to the properties of atomic i nteract ions (Sec. 
|III.A| >, the onset of quantum degeneracy (Sec. |VII.A| , and the 



The electronic ground states of hydrogen and alkali- 
metal atoms correspond to different orientations of the 
J = 1/2 electron and of the nuclear spin /. At low mag- 
netic fields, the hyperfine interaction between these spins 
dominates, and the ground-state subspace breaks into 
manifolds of states with definite total spin F — I db 1/2. 
At higher magnetic fields, the dominance of the electron 
magnetic moment causes the eigenstate structure to re- 
organize, as shown in Fig. [I] 

The collision properties of mixtures of different atomic 
hyperfine states are not always favorable to experi- 
ments on quantum gases. Atoms in the higher hyper- 
fine spi n manifold may unde rgo hyperfine-relaxation col- 
lisions ( jGorlitz et al.\ |2QQ3[ ) , leading to significant heat- 
ing and atom loss. In the case of 133 Cs, atoms in the 
lower hyperfine spin manifold also suffer from exother- 



mic c ollisions ( |Guery-Odelin et al. , 1998 Sod ing et al. 



1998) due to scattering resonances ( Arndt et al.\ 1997 

^nl/l/olmono o+ nl 1 QQSh Tn rvl-nor- r 



Kokkelmans et al. 1998[ ) . In other gases, such as 7 Li anc . 
8t) Rb, the s-wave scattering lengths for several collision 
channels are negative, implying that a Bose-condensed 
gas with large atom number acquires a negative corn- 



expected appearance of magnetic order in optical lattices (Sec. 
VII.C). 
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FIG. 1 (Color online) Ground-state hyperfine structure of 87 Rb. With a nuclear spin of I = 3/2, 87 Rb has a total angular 
momentum of either F — 1 or F — 2 in the J — 1/2 electronic ground state. At low magnetic fields, F is a good quantum 
number due to rotational invariance. The F — 1 and F — 2 spin manifolds are separated by the hyperfine interaction. The 
Zeeman sublevels, labeled by the magnetic quantum number tof, have a Zeeman energy that shows a linear variation with field 
at low magnetic fields (left), but deviates from the linear dependence strongly when the linear Zeeman shift is comparable to 
the hyperfine interaction energy (center). Right inset: the net quadratic Zeeman energy shift is exhibited after subtracting off 
the zero-field hyperfine splitting, the linear Zeeman term, and the energy of the \m z — 0) state; thus, the net quadratic shift 
of the \m z —) states is zero by definition. The quadratic shift due to a dc magnetic field is positive for the F — 1 levels and 
negative for the F — 2 levels. The hyperfine structure of other hydrogen-like atoms shows similar features. 



pressibility and collapses, leading to rapid decay (Donley 



et al. 2001 Sackett et al. 1998). 



Several compatible mixtures of internal states do re- 
main. The main focus of this Review is mixtures contain- 
ing all magnetic sublevels of a manifold of states with a 
single value of the hyperfine spin F. The internal state 
of an atom in such a manifold of states is represented as 
a spherical tensor of rank 2F in spin space; thus, such 
gases are called spinor gases. This spherical-tensor char- 
acter endows the spinor gas and the quantum fluid that 
it attains with an elegant phenomenology, as we describe 
in the remainder of this work. 



Experimentally, long-lived alkali spinor gases have 
been explored in three spin manifolds (Table the 
lower-energy F = 1 hyperfine spin manifol ds of both 23 Na 



( Stamper-Kurn et al. , 1998) and 87 Rb ( |Barrett et al. 
2001 ), and also the higher-energy F = 2 manifold of 87 Rb 



([Chang et al.\ |2004| |Kuwamoto et al.\ |2004| |Schmaljo- 



hann et al. 2004). The latter gas is moderately stable 
against hyperfine relaxation collisions on account of for- 
tuitous propertie s of the molecula r potential of the ru- 
bidium dimer (Burke et al. 1997). In comparison, the 



spinor gas in the F = 2 upper hyperfine levels of 23 Na 
was found to decay rapidly ( jGorlitz et al. ^003 ) . Several 
other promising systems still await experimental explo- 
ration, such as the F = 1 spinor gases of 7 Li and 41 K, 
which represent the lower of the two hyperfine spin man- 
ifolds for each atom, and perhaps also the higher hyper- 
fine spin manifolds of 41 K and of hydrogen. Experiments 



Stable 


Unstable 


(F z ) conserved 


(F z ) not conserved 


7 Li, F = 1 (f) 
23 Na, F = 1 (af) 
41 K, F = 1 (f) 
87 Rb, F = 1 (f) 
87 Rb, F = 2 (af or eye) 


52 Cr, F = 3 (not f) 
Dy, F = 8 (?) 
Er, F = 6 (?) 


7 Li, F — 2 
23 Na, F — 2 

39 K 

85 Rb 
133 Cs 


87 Rb pseudo-spin: 
|1,0>, |2,0) 
|1,±1>, |2,Tl> 


Tm, F = 4 (?) 



TABLE I Experimental candidates for the study of ultra- 
cold spinor Bose gases. Species are divided according to 
whether they are stable at zero magnetic field (information on 
thulium is lacking) , and whether the dipolar relaxation rate is 
small enough to allow the longitudinal magnetization ((F z )) 
to be conserved in an experiment. The nature of the spin- 
dependent contact interactions is indicated in parentheses (f: 
ferromagnetic, af: antiferromagnetic, eye: cyclic or tetrahe- 
dral, ?:unknown). Stable pseudo-spin- 1/2 gases of 87 Rb are 
indicated, with states labeled with quantum numbers |F, tuf) 
having the same low- field magnetic moment. 



with radioactive alkali isotopes are also possible 2 



2 In atoms with lifetimes longer than the typical experimental cycle 
time of around 10 s, the lower-spin manifolds of the electronic 
ground state offer F = 1 gases of 21 Na, 43 ' 45 K, 77,81,89,91^ 
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The collisional stability of 87 Rb also permits experi- 
ments on mixtures of atoms from both the F = 1 and 
F = 2 hyperfine spin manifolds (Myatt et al. 1997| ). In 
particular, the |F = l,m z = —1) and [F = 2, ra 2 = +1) 
states possess magnetic moments that are nearly equal 
at low magnetic field (identical at a field of 3.23 G). 
Both species can be confined in magnetic traps, in which 
the external potential is nearly state-independent. The 
magnetic-field insensitivity of the internal energy dif- 
ference between these states allows the coherence be- 
tween them to be probed over seconds-long evolution 
times ( [Deutsch et a/.l[20Tol|Harber et al. ||2002||Treutlein| 



et al., 2004). Similarly, mixtures of the magnetic-field- 
insensitive \F = l : m z = 0) and \F = 2, m z = 0) states of 
87 Rb a re stable in optical traps, an d show long coherence 
times ( |Kleine Biining et aZ.[ [2Qll] ). In contrast, for the 
spinor gases described above, the linear Zeeman energy 
shift between sublevels makes this coherence difficult to 



access, as discussed further in Sec. IV. C The two internal 



states of such pseudo- spin- 1/2 Bose gases are not related 
by rotational symmetry in spin space, and so one expects 
phenomena associated with ground-state degeneracy and 
symmetry breaking to differ from those of spinor gases. 
These differences are explored in Sec. III.C 

Other species of ultracold atoms comprise the 
metastable 3 Si states of helium and the noble gases, and 
the 3 Pj states of two-electron atoms such as the alkali- 
earth elements and Yb. For bosonic isotopes of these 
elements, all of which have a nuclear spin 1 = 0, the elec- 
tron angular momentum may be non-zero. Thus, held 
in optical traps, such atoms offer another realization of 
spinor Bose gases. However, the met ast ability of these 
atoms relies on the suppression of collisional decay to 
the electronic ground state due to spin polarization. Ex- 
periments on metastable 4 He demonstrated that s-wave 
collisions between a pair of atoms with total projected 
angular momentum of M p a j r = lead to strong inelas- 
tic losses ( |Partridge et al] 2010). Thus, it appears that 
stable spinor quantum fluids of such gases are precluded. 



2. Many-electron atoms: Spinor gases with strong dipolar 
interactions 

In condensed-matter systems, magnetic order is estab- 
lished both by local interactions, as in the Heisenberg- 
Dirac model of magnetism, and also by long-range inter- 
actions, such as the magnetic dipole interaction. Simi- 
larly, as we discuss below, magnetic order in degenerate 
spinor Bose gases is influenced both by the local spin- 
dependent contact interactions and by the long-range 
magnetic dipole-dipole interaction (MDDI). The strength 
of dipolar interactions increases with the square of the 



121 CS, and 223,225 Fr . p = 2 gaseg Qf 25 Na5 79,83 Rb? 131^ and 

22i Fr; F = 3 gageg of 135,137,139,141 Cs . anc i F = 4 gases of 119 Cs 
and 207 > 209 > 211 > 213 Fr 



magnetic moment fi of the atoms forming the spinor gas. 
For the ground-state alkali gases discussed above, ne- 
glecting the small nuclear contribution, the atomic mag- 
netic moment is \i = hbF/(I + 1/2), i.e. no larger than 
one Bohr magneton \±b- Dipolar effects are visible in 



such g ases (|Kawaguchi et al] |2006b[ 2007 Vengalattore 
Yi and Pu[ 2006| , although the local mag- 



et al. , 2008 



netic order is still selected predominantly by the contact 
interaction. 

Both the long-range and the local spin-dependent in- 
teractions are varied strongly in gases with larger mag- 
netic moments. For example, quantum degenerate gases 
of 52 Cr disp lay strong dip olar effects such as den sity de- 
forma t ions (|Lahaye et a/." , 2007[ ) and instability ( |Lahaye| 
et al, 2008) due to the anisotropic long-range dipolar 



interaction. The magnetic moment of this spin-3 gas is 
yielding dipolar energies that are at least 36 times 
higher than those in the alkali gases (144 times larger 
than in F = 1 87 Rb). Several lanthanide elements have 
been laser cooled, presenting new candidates for high- 
spin and highly magnetic spinor Bose gases. Thulium 
(Sukachev et al, 2010a|b[ ) carries a total spin of F = 4 
in its ground state, with the magnetic moment of 4/i#. 



Stronger dipolar interactions are realized in Er (Berglund 
20071 |McClelland and Hanssenl |2006[ ) and Dy 



et al. 



^Youn et al. 2010| ), in which the atomic angular mo- 
menta of bosonic species are J — 6 and J = 8, leading to 
magnetic moments of 1\±b and 10/i#, respectively. Bose- 
Einstein condensates of both Er and Dy have been pro- 
duced, and t hey show signs of strong magnetic dipolar 
interactions (Aikawa et al., 2012 Lu et al., 2011). 



The large electron spin of these gases also influences 
the short-range interactions between atoms. The pres- 
ence of several molecular potentials, split by the fine 
structure for different configurations of the total elec- 
tronic spin of the colliding atoms, and the strong dipo- 
lar interaction at molecular distances, altogether lead 
to much higher ra tes for dipolar relaxation among the 
Zeem an sublevels ( Hensler et al. , 2003| |Weinstein et al.\ 
1998) than in the alkali gases. The phenomenology of 



spinor gases for these highly magnetic elements is thus 
expected to be different. In the alkali gases, the total 
spin of the gas along the direction of an applied mag- 
netic field is typically constant over the duration of an 
experiment. In contrast, the longitudinal magnetiza- 
tion of high-magnetic- moment spinor gases is dynamic. 
The enhanced dipolar relaxation couples the spin and 
mass currents of spinor Bose gases, leading to phenom- 



ena such as the Einstein-de Haas effect ( Gawryluk et al. 



2007| |Kawaguchi et al] |2006a| |Santos and Pfau| |2006 ). 
The anisotropy of the dipolar interaction also favors low- 
energy states in which the spin-state of the condensate 
is not spatially uniform; rather, the spin order parame- 
ter varies spatially, producing what is known as a spin 
texture. T his spin texture spontaneously breaks chira l 
symmetry (|Kawaguchi et al] 2006b Yi and Pu, 2006). 



Strong dipolar relaxation also requires that spinor-gas 
experiments with these elements occur at precisely con- 
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trolled, weak magnetic fields. Such experimental condi- 
tions have recently been achieved, allowing fo r the first 
observation of spinor gas physics with 52 Cr (Pasquiou 



et al 
in Sec 



2011) 
. IVIT 



Dipolar effects in spinor gases are discussed 



3. Orbital degrees of freedom in optical lattices 

The physics of spinor Bose gases may also find re- 
alization in systems of atoms in optical lattice poten- 
tials. Similar to the orbital states of electrons within the 
unit cell of a periodic crystal, the quantum states of the 
center-of-mass motion within a unit cell of a periodic op- 
tical lattice can be regarded as an internal state space 
for a lattice-bound atomic gas. These motional states 
may have non-zero angular momentum. In deep opti- 
cal lattices where the lattice site potentials tend toward 
spherical symmetry, the system closely resembles a spinor 
Bose gas insofar as the spherical symmetry implies the 
presence of ground-state degeneracies and a particular 
structure to the interactions between atoms. As the lat- 
tice depth is reduced, the orbital degeneracy is strongly 
influenced by the structure and discrete symmetries of 
the lattice. 

A degenerate Bose gas in the superfluid state in such 
a shallow lattice may be regarded as analogous to a bulk 
spinor Bose-Einstein condensate. Both the kinetic and 
interaction energy in the lattice can depend strongly on 
the orbital quantum state. A simple example of this 
dependence arises from the anisotropic tunneling of p- 
orbital particles in a primitive lattice, with the tunneling 
rate being higher for the p-orbital aligned along the di 
rection of tunneling ( |Isacsson and Girvin , 2005 Liu and 
Wu 2006[) . In terms of orbital magnetism in a bosonic su- 



perfluid, such tunneling gives anisotropic stiffness of the 
magnetic order parameter, which may lead to ground- 
state "orbital textures" with circulating orbital currents, 
reminiscent of the spin textures that are expected to arise 
from the anisotropic spin interactions of dipolar spinor 
condensates. Such lat tice-gas systems have been realize d 
experimentally (iMuller et al] 120071 iWirth et al] 12011). 



The interplay of orbital dynamics and spin dynamics in 
non-Bravais lattices has also been considered (Wagner 



et al. 2011). However, one must note that orbital mag- 



netism comes about from the excitation of high-energy 
bands, and, at least for Bose gases, is generally not pro- 
tected against decay to the ground band; in contrast, 
at low magnetic fields or in the absence of dipolar relax- 
ation, non-trivial spin ordering can persist for long times. 

B. General properties of a quantum Bose liquid with 
internal symmetry 

Consider a scalar Bose-Einstein gas in a box potential. 
As the temperature of the gas is lowered into the regime 
of quantum degeneracy, why does the condensate form 
in the motional ground state of the box? 



The answer lies in the influence of bosonic particle- 
exchange symmetry on the enumeration of microstates 
available to the system. In particular, in comparison with 
Maxwell-Boltzmann statistics of distinguishable parti- 
cles, Bose-Einstein statistics greatly reduces, by a factor 
of TV!, the number of nondegenerate quantum states avail- 
able to an iV-body system. Thus, Bose-Einstein statis- 
tics makes states with macroscopic populations in single- 
particle states more likely. Once the system temperature 
T is reduced so that the number of single-particle states 
M accessible to the system (roughly the number of states 
with energies below fe^T) becomes smaller than A/", mul- 
tiply occupied single-particle states become the norm. In 
a d-dimensional box potential, this degeneracy occurs at 
ksT c ~ EiN 2 / d where E\ is the energy of the first ex- 
cited single-particle state in the box. Note that a huge 
number of single-particle states exist below the energy 
scale /cbT c , yet a macroscopic number of particles con- 
dense into the lowest single-particle state due to bosonic 
enhancement. 

To address the issue posed above, we consider now the 
population of atoms in the first excited single-particle 
state once the system is cooled below T c . This popula- 
tion is given as Ni < k B T c /Ei = N 2 / d , so that, in three 
dimensions, with the scaling exponent 2/d < 1, only the 
ground state is macroscopically occupied. In contrast, 
for a Maxwell-Boltzmann gas, the dominant macroscopic 
occupation of the ground state occurs only at the tem- 
perature Tmb = Ei/ks- Comparing the two situations, 
we see that Bose-Einstein statistics have the effect of 
greatly enhancing the effects of extremely small energy 
differences (here, E\) so as to allow the ground state to 
be distinguished already at temperatures that are N 2 / d 
higher than Ei/ks- 

In degenerate spinor Bose gases, one finds that the 
spin-dependent contact and dipolar interactions distin- 
guish between different spin states of the spinor order 
parameter. The energy differences per particle are of- 
ten minute, far smaller than the temperatures accessed 
in experiments; e.g. for the F = 1 gas of 87 Rb, the 
spin-dependent contact interaction contributes roughly 
1 nK of energy per particle, far lower than the several 
100 nK critical temperature for Bose-Einstein condensa- 
tion. And yet, due to the quantum-statistical Bose en- 
hancement effect described above, the gas becomes mag- 
netically ordered just below the Bose-Einstein conden- 
sation temperature. Such Bose-enhanced magnetization 
is a novel mechanism for magnetism, distinct from the 
more conventional mechanisms for localized or itinerant 
electrons in solids. 

However, such Bose enhancement still does not deter- 
mine uniquely the state of a quantum degenerate spinor 
Bose gas. In the absence of explicit symmetry break- 
ing, e.g. due to applied magnetic fields or non-spherically 
symmetric trapping potentials, the internal rotational 
symmetry of the spinor Bose gas guarantees that sev- 
eral single-particle spinor states achieve the same low- 
est energy of the interacting system. The set of spinor 
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ground states for the spin-F gas may not retain the full 
SU(2F + 1) symmetry owing to spin-dependent interac- 
tions, but the lower-symmetry manifold of spinor ground 
states is still guaranteed to be degenerate over a contin- 
uous symmetry group. 

In condensed-matter magnetic systems, such degenera- 
cies are often resolved by experimental actualities: ex- 
traneous magnetic bias fields, magnetic defects, stress 
anisotropy, etc. In spinor Bose gases, many of these 
artifacts are absent. Does the Bose condensed system 
now break symmetry spontaneously and choose a con- 
ventional Bose-Einstein condensate, with just a single 
macroscopically occupied single-particle state (Penrose 



and Onsager, 1956) but with random magnetization? Or 
does the system abide by the postulate of molecular chaos 
- assuming there are no constraints that bar the sys- 
tem from thermalizing - and form an unbroken symme- 
try state with macro scopic populations in severa l single - 
particle spin states ( jNozieres and Saint James} |1982[ ) ? 



This question is taken up in Sec. V.B 



The presence of a degenerate ground-state manifold 
also influences the dynamics and structure of the degen- 
erate spinor Bose gas. In a broken-symmetry state, the 
gas will possess one or more gapless Goldstone modes. 
The symmetry group of the ground-state manifold deter- 
mines the range of topological defects that may be found 



in low-energy spin textures, as discussed in Sec. VI. D 



C. Outline of this review 

The central characteristic of spinor Bose gases is their 
response to geometric transformations such as rotations 
and inversions. To visualize the response of a high-spin 
object to such transformations, in Sec. [II] we outline a 
method, first developed by Majorana and subsequently 
re-introduced for the description of spinor gases, in which 
quantum states are represented by the positions of sev- 
eral points on the unit sphere; for a spin- 1/2 state, this 
method reproduces the familiar Bloch-sphere representa- 
tion. We also introduce a family of high-symmetry spin 
states, known as inert states for the reason that they 
represent extremal energy states for generic rotationally 
symmetric spin Hamiltonians. 

We then begin the physical description of spinor Bose 
gases. Section [TTT] illustrates how rotational symmetry is 
used to simplify enormously the characterization of inter- 
particle interactions. After specifying several approxima- 
tions, we find that the spin-dependent interactions within 
the spinor gas are determined by the relative strength of 
the s-wave scattering lengths for collisions among parti- 
cles with a given value of the total spin F pa i r of the col- 
liding pair. This simplification sets the spinor gas apart 
from other multi-component Bose gases, as we discuss for 
the experimentally relevant pseudo- spin- 1/2 system. 

is 



The description of atomic interactions in Sec. Ill 



tions under which spinor gases are prepared and studied 
depart from such ideal conditions. Thus, before contin- 
uing with our exploration of the physics of spinor Bose 



gases, we use Sec. IV to describe these experimental con- 
ditions and techniques. 

In Sec. [V| we describe the spatially uniform ground 
states of spinor Bose gases as determined by their spin- 
dependent s-wave interaction energies. This discussion 
includes a review of experimental observations that char- 
acterize the ground states of the F = 1 23 Na spinor gas 
(polar), the F = 1 87 Rb spinor gas (ferromagnetic), and 
the F = 2 87 Rb spinor gas (either antiferromagnetic or 
cyclic/tetrahedral). Ground states are determined not 
only via the mean-field (Hartree) approximation, but also 
via a many-body quantum treatment that admits ground 
states that are fragmented. 

For a non- uniform degenerate spinor Bose gas, long 
length-scale variations of the local magnetic order param- 
eters produce low-energy spin textures. In Sec. (Vl} we 



discuss salient properties of these textures. We highlight 
the connections between spatial variation of the spin- 
state (which may represent spin currents) and the super- 
fluid velocity (i.e. mass currents), which are required due 
to the invariance of certain order-parameter manifolds to 
combinations of geometric rotations and gauge transfor- 
mations. We also summarize our knowledge of topologi- 
cal excitations that are either observed or are predicted 
to occur in such spin textures. The coupling of spin- and 
mass-currents is also inherent to several hydrodynamic 
descriptions that have been developed to describe the 
low-energy dynamics of spin textures. 

Another key characteristic of spinor Bose gases is the 
onset of magnetic ordering due to Bose-Einstein statistics 
and the enhancement of small single-particle energy dif- 
ferences in the quantum degenerate gas. What happens 
to such magnetic ordering when the proclivity to Bose- 
Einstein condensation is weakened, such as at non-zero 
temperature, in two- or one-dimensional systems where 
finite temperature Bose-Einstein condensation does not 
occur, or in periodic potentials where the role of particle 
interactions is enhanced? These questions are considered 
in Sec. EII1 

VIII[ we consider the influence of 



Next, 



Sec. 



strictly valid only when the full rotational symmetry of 
the gas is preserved. In fact, the true experimental condi- 



the magnetic dipole-dipole interaction on the low- 
temperature spinor Bose gas. Such interactions add new 
features to the spinor gas: a new source of spin-dependent 
interactions, which can compete with the contact inter- 
actions to change the local magnetic order of the gas; 
and a new long-range interaction that competes with the 
kinetic-energy-inducted magnetic stiffness of the spinor 
gas. Intricate spin textures are predicted to develop in 
the ground state to reduce the dipolar energy. Dipolar 
effects are expected to strongly influence the nature of 
spinor gases composed of atoms with large magnetic mo- 
ments. However, even for the less magnetic alkali spinor 
gases, the long-range nature of the dipolar interaction is 
expected to be significant, and we describe experimental 
evidence for its effects. 
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In the next two sections, we discuss the dynamics of 



degenerate spinor Bose gases. Sec.[TX]is devoted to coher- 
ent spin mixing. We summarize the observations of spin- 
mixing oscillations in both microscopic and macroscopic 
spinor gas systems. Spin mixing collisions also induce dy- 
namical instabilities that can be controlled by varying the 
single-particle energies (via the quadratic Zeeman shift) 
and the initial state of the degenerate spinor gas. These 
instabilities provide access to a rich variety of physical 
phenomena, such as symmetry breaking, quantum atom 
optics, and spin squeezing. 

Whereas the early-time dynamics of a non-equilibrium 
spinor condensate may exhibit well-controlled coherent 
spin mixing, at longer times one expects the gas to evolve 
toward thermal equilibrium via phase-ordering kinetics. 
In Sec. [Xj we summarize several areas where coarsening 
and equilibration have been examined. 

Following Sec. |XI[ in which we demonstrate how the 
properties of ultracold spinor Bose gases may be use- 
ful for applications in magnetic sensing and imaging, we 
conclude by suggesting several open questions to be ad- 
dressed in future research (Sec. XII). 

Several previous reviews are complementary to this 
work. The early round of experiments on spinor 
Bose-Einstein condensates of sodium is summarized by 



review is presented by 



Stamper-Kurn and Ketterle (2001). A more up-to-date 



Kawaguchi and Ueda (2012), and 



also briefly by Ueda (2012), with a focus on theoreti- 
cal frameworks. Spinor gases are also treated briefly in 



Lewenstein et al. (2007), with a focus on gases in op- 
tical lattices. In addition, several of the topics covered 
in this Review have been discussed separately. Several 
reviews on properties of Bose gases with strong dipolar 
interactions focus on gases with uniformly oriented dipole 



moments (Baranov et al. 2002 Baranov 2008 Lahaye 



et al] 2009). In contrast, the present focus is on dipolar 
interactions in spinor gases with a variable dipole orien- 
tation (Sec. VIII). Aspects of topological defects, partic- 
ularly of vortices, in spinor condensates are discussed in 
|Kasamatsu et al ( |2005[ ). Several reviews have summa- 
rized the explorations of the Kibble- Zurek mechanism in 
condensed-matter systems, both for thermal phase tran- 



sitions (Arodz et 



2003 



Zurek[ 1996 ) and, rece ntly, 
for quantum phase transitions (Dziarmaga, 2010). A 
broader range of non-e quilibrium phenomena is discussed 
in the Colloquium by Polkovnikov et al. (2011). 



II. QUANTUM STATES OF A SPIN-F OBJECT 

We begin by introducing and relating several methods 
for characterizing the states of a spin-F object. 



A. Geometric representation of spinor wavefunctions 

Majorana invented a geometric representation of a gen- 
eral spin-F state (Majorana, 1932), which helps illustrate 



the symmetries retained by a spinor wavefunction inde- 
pendent of the special coordinate system chosen. This 
representation can be used to picture the mean-field 
ground states, dynamics, and topological structures of 



degenerate spinor Bose gases (jBarnett et al.\ [2008 , 


12009 


2006, 


2007; Kawaguchi and Ueda[|2011| 


Lamacraft 


2010 


Make 


a and Suominen, 2007; Marcin anc 


Krzysztof 


2012 


Turner et al. t 2007 Turner and Demler 


2009). 



be represented pictorially as a point on the unit sphere 
(the Bloch sphere). That is, given a spin- 1/2 state |^), 
one can find a projection of the Pauli spin vector, — 
& • 7r where 7r is a unit vector, so that is the +1 
eigenstate of a n . Apart from an overall phase, the state 
|^) is then represented by a point on the Bloch sphere at 
the position tt. Alternately, one can instead consider the 
family of fully magnetized states, defined by the relation 
a -tz'^') = We then identify the fully magnetized 

state that is orthogonal to |^), thereby identifying a unit 
vector tt' that equivalently represents the state on the 
Bloch sphere. In this case we see trivially that 7r = — tz' . 

Now we extend this pictorial scheme to states of higher 
spin. First, let us imagine that a spin-F particle is ac- 
tually composed of 2F spin- 1/2 sub-particles, and that 
F is the quantum number of the total spin of these sub- 
particles. A spin-F wavefunction, written in the tensor 
space of these spin-1/2 sub-particles, is symmetric un- 
der the exchange of any two of the sub-particles. The 
state of each of these sub-particles defines a unit vector 
7r r (r G {1, 2, • • • 2F}) in the manner described above for 
spin-1/2 states. This set of vectors corresponds to 2F 
points on the unit sphere, providing a pictorial represen- 
tation of the spin-F state apart from an overall phase. 

Alternately, we ask whether the spin-F state is 
orthogonal to the maximum- value eigenstate of F • 7r', 
the projection of the spin vector along the unit vector 
7r f whose angular coordinates are (0' ', <fi'). Expressing the 
state in the F z eigenbasis as typ, i/jf-i, • • • ^-f) T 
(T denotes the transpose), the orthogonality condition 
produces the following polynomial equation: 



2F-1 



b 2F = 0, 



(1) 



where bj = ^p-jl V ~ an d where £ is the stere- 
ographic projection on the equatorial complex plane of 
the point (0 f ,(/)') from the south pole: 



tan ((972) e^ = ( 



(2) 



The 2F complex roots of the polynomial equation, 
Ci 5 C2 , • • * identify the set of unit vectors 7r' r . As 

in the spin-1/2 case, these vectors are antipodal to the 
unit vectors obtained from the spin-1/2 sub-particle ap- 
proach, i.e. 7r r = —iz' r . 

In this Review, we adopt the convention that a spinor 
is represented by the set 7r r . The state described 
by the unit vectors iz' r is generated from by time 
reversal, by which one reverses the orientation of each of 
the spin-1/2 sub-particles of which \ip) is composed. 
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B. Spin moments 

The state of a spin-F object, and of the magnetic or- 
dering produced in an ensemble of such objects, may 
also be described by its spin moments. Since the single- 
particle density matrix of a spin-F system is a (2F + 1) x 
(2F + 1) square matrix, its full description requires the 
specification of spin moments up to rank 2F. The rank-1 
moment is the expectation value of the spin vector oper- 
ator F. The rank-2 moment is the expectation value of 
the spin quadrupole tensor, defined as 



-F 2 5 
3 



(3) 



where fi and v specify Cartesian axes. This tensor, and 
also its density as defined below, are sometimes referred 
to as the nematicity. Throughout this Review, both the 
vector spin operator F and the spin quadrupole tensor 
N^ v are taken to be dimensionless, with the quantum 
of angular momentum h absorbed into the definitions of 
other terms. 

We shall also use the symbols T^{v) and A/^fr), writ- 
ten now as functions of the spatial position r, to define 
the densities of the spin vector and spin quadrupole mo- 
ments. That is, let us define the bosonic field operator 
i/j mF (r) for particles in the eigenstate of the F z operator 
with magnetic quantum number m^. The \i component 
of the spin vector density is then given as 



^(r) 



(4) 



where F^ is the matrix representation of the correspond- 
ing projection of the spin vector operator in the F z eigen- 
basis. The components of the spin quadrupole tensor 
density are written similarly. Both quantities have the 
units of number density. The magnetization M(r) is 
found by multiplying the spin vector density by the max- 
imum magnetic moment of a spin-F atom. 

The spin moments may be expressed in terms of the 
unit vectors 7r r defined in the Majorana representation, 
although such expressions beco me complicated for h igher 
spin and higher-rank moments (Barnett et al t 2009). For 
the simple case of F = 1, the vector spin is found to be 



(F) 



2(7ri +7T 2 ) 



3 + 7Ti • 7T 2 

and the spin-quadrupole tensor as 



(5) 



(N) = 



1 



(6) 



6 3 + 7Ti • 7T 2 

(7Ti + 7r 2 )(7Ti + 7T 2 ) + (7Ti X 7r 2 )(7Ti X 7T 2 ) 
1 + 7Ti • 7T 2 



where the matrices in the square brackets are written as 
dyadic tensors, and I is the identity matrix. 



Let us consider two examples of Majorana represen- 
tations for the simple case of F = 1. For the spinor 
(1/2, 1/V% 1/2) T , where T denotes the transpose, the 
polynomial ( 2 +2£+l = has a double root at £ = — 1, so 
that the state is represented by two overlapping points at 
the location (6 = 7r/2, <j> = 0) and has the SO (2) symme- 
try about the x axis. Geometric rotations of this state are 
represented by two overlapping points at the correspond- 
ing positions on the sphere. Majorana representations 
of longitudinal ferromagnetic states are shown in Fig. 
|5j States within this manifold maximize the vector spin 
|(F) |, and so they are known as "ferromagnetic" states. 
Note that the time-reversed spinor state is obtained by 
moving the points representing the state to their antipo- 
dal positions; indeed, the ferromagnetic states reverse 
their orientation under time reversal. 

Another manifold of states is obtained starting from 
the spinor ip = (1/V% 0, -l/<\/2) T , for which Eq. (l) 
gives roots at ( — 1 and ( = — 1. The corresponding 
points are located at antipodal points on the unit sphere; 
an example is shown in Fig. [5] For this state, we identify 
x as the nematic director n, i.e. the axis along which the 
spinor is the eigenstate of the spin projection F • n with 
eigenvalue zero. The state has both 50(2) symmetry 
about the x axis and also 7r rotation symmetry about 
any axis on the y-z plane. The states obtained from this 
starting state by rotations are known as "polar" states. 
They are characterized by zero magnetization, while their 
inversion symmetry is reflected in their nematicity. 

Representations for higher-spin states can be similarly 
obtained. Spin-2 states are represented by four points on 
the unit sphere. Examples of such states relevant to the 
zero-field mean-field ground states of F = 2 condensates 
(Fig. [7]) include the ferromagnetic state, represented by 
four overlapping dots; the uniaxial nematic state, repre- 
sented by two doubly degenerate dots at antipodal points 
on the sphere; the biaxial nematic state, represented by 
points lying on a great circle at the vertices of a regular 
square; and the cyclic/tetrahedral phase 3 , represented by 
points at the corners of a regular tetrahedron. 



Spin- 3 states are represented by six points on the uni t 
sphere ( |Barnett et a/.||2006 ; Kawaguchi and Ueda 2011 ). 
The possible symmetries constructed from six vertices are 
octahedron, tetrahedron, dihedral-2,3,4,5,6 and 02,3,4,5,6, 
where C n represents the cyclic symmetry of order n. The 
icosahedron symmetry is expected to appear for F > 6. 



The term "cyclic" was introduced by jCiobanu et a/.] ( |2000| , refer- 
ring to earlier classifications of the states of d-wave superconduc- 
tors. The term may also refer to the 3-fold rotational symmetry 
of the state. The Majorana representation of the state clarifies 
that its symmetry group is that of the tetrahedron, so a better 
name for the state may be "tetrahedral." We use both names in 
this Review in keeping with previous usage. 
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C. Inert states 



A. Quantum scattering under rotational symmetry 



There is a special class of states, called "inert states," 
which is defined to be stationary for generic rotation- 
ally symmetric, spin-dependent energy functionals and 



ton and Moore[ 19741 19751 Makela and Suominen 2007; 


Vollhardt and Wolfle 


1990 Yip 


2007). There is a 



close relationship between the symmetry of the order pa- 
ramete r and a statio nary point of the free-energy func- 
tional (Michel, 1980). This relationship greatly simplifies 



the otherwise rather involved calculation of stationary 
states. Let g G G be an element of the group of opera- 
tions on the order parameter ip of the system that leave 
the free-energy functional invariant: f[gip] = f[ip\- Let 
h G H be an element of the isotropy group H that leaves 
ip unaltered: hip = ip. The order parameter manifold M 
is then given by a coset M = G/H. Now take one ip and 
operate every g G G on it. The resulting set of gip is said 
to constitute an orbit. Now take another ip' and find the 
corresponding orbit by operating every element of G. If 
two orbits share the same isotropy group, they are said to 
belong to the same stratum which is defined as the union 
of orbits with the same isotropy group 



Michel 



(1980) 



asserts that if an orbit is isolated in its stratum, i.e., if 
it does not have another orbit in its neighborhood, the 
corresponding order parameter is a stationary function of 
the free-energy functional. This theorem is useful to find 
inert states because they depend only on the symmetry 
of the system and are independent of the details of the 
system's parameters. As an illustration, let us consider 
the following eigenvalue problem: e l ^e~ l ^ z ^ip = ip. For 
a spin-1 BEC, the solution is ip = (1,0,0) T , which im- 
plies that the ferromagnetic state is an extremum of the 
mean-field energy. 

The ground states of spinor Bose condensates in the 
absence of externally imposed symmetry-breaking fields 
are i nert, as are the A, Ai, and B phases of superfluid 



He (IVollhardt and Wolflel 19901. It can be shown that 



the ferromagnetic and polar states are the only two in- 
ert states of the spin-1 condensate. For the spin-2 con- 
densate, the ferromagnetic, uniaxial nematic, biaxial ne- 
matic, and cyclic/tetrahedral states are all inert states. 



III. BOSE GASES WITH INTERNAL ROTATIONAL 
SYMMETRY 



In high-energy, condensed- matter, atomic and molec- 
ular physics alike, the starting point for understanding a 
new system is the recognition of its underlying symme- 
tries, by which one can simplify its energetics, dynamics, 
topological structures, and so on. In this spirit, we eluci- 
date how rotational symmetry constrains and simplifies 
the interactions present in a spinor Bose gas. 



At the low densities of atomic gases, interactions re- 
sult almost exclusively from binary collisions and, even 
in a many-body system, such interactions can be under- 
stood through the treatment of the two-body quantum 
scattering problem. Neglecting the effects of particle con- 
finement, we treat such scattering in three spatial dimen- 
sions. Collisions are then characterized by the scattering 
matrix that connects the incoming and outgoing asymp- 
totic states. Such states are chosen as product states of 
the orbital motion and internal states. We specify the 
former in the spherical harmonic basis by the relative 
wave number fc, and the quantum numbers for the to- 
tal orbital angular momentum L pa i r and its projection 
on the quantization axis ra£ )Pa i r . We specify the latter 
using basis states of spinor wavefunctions ip for the col- 
liding particles. Time-reversal symmetry allows us to in- 
terchange the incident and outgoing quantum states. En- 
ergy conservation (invariance under translations in time) 
ensures that differences between the incident and outgo- 
ing kinetic energies, determined by ki and kf, must be 
compensated by changes in the internal energies of the 
spin states in the collision; e.g. dipolar relaxation in an 
applied magnetic field converts the Zeeman energy into 
the kinetic energy of the colliding particles. 

The next step in simplifying the scattering problem 
is to consider the character of the two-body interaction, 
which separates fairly cleanly into short-range and long- 
range parts, with the divide between these parts occur- 
ring at an internuclear distance r longer than the effective 
range of the potential tq. The long-range part is ulti- 
mately dominated by the magnetic dipole-dipole inter- 
action (MDDI), which falls off as 1/r 3 . This functional 
dependence requires that the dipole interaction be re- 
tained explicitly in the interaction Hamiltonian. For the 
sake of simplicity, we neglect the long-range portion at 
this point, deferring the treatment of dipolar interactions 
to Sec. IVlTTl 

The short-range portion of the potential is complicated 
as it comes about from the interaction among all the con- 
stituent particles within the colliding atoms. This por- 
tion also contains terms that mix the spin and orbital an- 
gular momenta of the colliding particles, such as through 
second-order spin-orbit coupling. The short-range inter- 
action also influences the strength of the MDDI; for ex- 
ample, a scattering resonance causes the two-body wave- 
function to be enhanced at short internuclear distances, 
enhancing the effect of dipolar interactions. 

To treat the quantum scattering due to the short-range 
potential, several approximations are commonly made. 
First, we assume that the incident collision energy is very 
low, in the sense that the thermal de Broglie length of 
the relative motion of the atoms obeys Xjb ^> ^o- I n this 
cold-collision approximation, only the lowest-order inci- 
dent partial waves (L pa ir,i = 0) undergo collisions (sub- 
script i indicates "initial" ) ; all other partial waves can be 
ignored. 
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Second, we make the "spinor gas collision" approxima- 
tion that the short-range potential is rotationally invari- 
ant - such rotations affect both the internal spin degrees 
of freedom and also the orbital spatial degrees of freedom. 
This approximation is exact in the absence of any exter- 
nal source of rotational symmetry breaking such as ap- 
plied magnetic fields, non-spherical trapping potentials, 
vector or tensor ac Stark shifts in the optical dipole force, 
etc. In practice, such influences are rarely completely 
absent; nevertheless, the spinor gas collision approxima- 
tion still applies at low magnetic fields, away from mag- 
netic or optical Feshbach resonances and in the normal 
Zeeman regime of the ground-state hyperfine structure. 
This approximation guarantees the conservation of the 
total angular momentum of the colliding pair, which is 
the sum of the total orbital angular momentum Lp a ir 
and the internal angular momentum (nuclear and elec- 
tronic) F pa i r . By the cold-collision approximation, the 
total angular momentum of a colliding pair is -F pa ir,i- 
However, each source of angular momentum need not 
be separately conserved. Notably, dipolar relaxation can 
convert internal to orbital angular momentum; to first 
order, such a collision promotes the colliding pair from 
an incident orbital s-wave (L pa ir,i = 0) to an outgoing 
d-wave (L pair j = 2) wavefunction, with a concomitant 
change in F pa i r j, where subscript / indicates "final." 

To simplify further, one commonly makes a third as- 
sumption, the "weak-dipolar approximation," whereby 
one neglects the spin-orbit coupling through the short- 
range molecular potential. Now, the orbital and internal 
angular momenta are separately conserved: for a col- 
liding atom pair, L pair ^ = L pair j = and F pair)i = 

-^pair,/ = -^pair- 

A final approximation is that s-wave collisions do not 
mix total hyperfine states of the colliding atoms. For 
example, in considering the properties of a spinor gas 
in an upper hyperfine manifold of the electronic ground 
state (say the F = 2 spinor gas of 87 Rb), we neglect 
collisions by which atoms undergo hyperfine relaxation 
to the lower manifold (F = 1). 

Only after all these approximations can we conclude 
that the myriad collisions among spinor-gas atoms are 
characterized simply by the s-wave scattering lengths 
QF, P air between two particles in the colli sion channel F pa i r 
po| |1998| |Ohmi and Machida[ |1998[ ). The parity of 
F pair is constrained further by quantum statistics. The 
many-body wave function of identical spin-F particles 
changes by a factor of (— 1) 2F under exchange of any two 
particles. By the same exchange, the spin and orbital 



parts of the wave function change by ( — 1) 



r+2F 



and 



(- 



(-i) 



d 2f 



respectively. To be consistent, we must have 



(-1) 



r + 2F 



Spinor gas 


Scattering lengths [as] 


7 Li, F = 1 
23 Na, F 1 

41 K, F=l 
87 Rb, F = 1 

87 Rb, F — 2 
52 Cr, F = 3 


a 23.9, a 2 6.8 
a = 50.0 ± 1.6, a 2 = 55.0 ± 1.7 
a = 68.5 ± 0.7, a 2 = 63.5 ± 0.6 
a = 101.8 ± 0.2, a 2 = 100.4 ± 0.1 
a = 87.4 ± 1.0, a 2 = 92.4 ± 1.0, 
a 4 = 100.5 ± 1.0 

a = 30 - 50, a 2 -7 ±20, 
a 4 = 58 ±6, a 6 = 102.5 ± 0.4 



TABLE II Scattering lengths a i? pair for several realizations 
of spinor Bose gases. These scattering lengths are calculated 
theoretically, using numerical models for the molecular po- 
tential with input parameters that are determined by com- 
parison with experimental measurements of the positions of 
Feshbach resonances, of the line strengths of photoassociation 
resonances, and of scattering length differences surmised from 
the dynamics of quantum gases. Values for lithium are from 
?, for sodium from Crubellier et al. ( 1999 ) , for potassium from 
|Falke et oT\ ( |2008| ) and |Lysebo and Veseth| fl2010|) (we guess 
errors at t he level of 1%J, for rubid ium from jKlausen et al' 
(2001 ) and|van Kempen et a l. (2002), and for chromium from 
Wern er it al\ ( |2005| ), |Pasquiou et al. (2010) and the value of 
ao is due to H. Sadeghpour (private communication). Val- 
ues and standard error estimates (where they are known) are 
given in units of the Bohr radius as- 



Ohmi and Machida, 1998) 



T *- r j) 2_. Papain (7) 



i,j evenF pair 



where M is the atomic mass and the sum is taken over 
particle pairs, labeled i and j. The operator Pp,pair 
projects a pair of atoms - implicitly the pair - into 
the total spin-Fp a i r state, where only even values of F pa ir 
are now considered. 

Table [II] summarizes the known scattering lengths 
&F,pair fc> r several spinor Bose gases that are experimen- 
tally accessible. As discussed below, and in Sec. V.A[ the 
different relations among the scattering lengths for differ- 
ent Fpair yield different predictions for the spin ordering 
of the ground state. 



(— l) Lpair ; hence F l 



pair 



L 



pair 



must be even. For s-wave collisions, F pair must be even. B - Spin-dependent s-wave interaction energy 

This conclusion holds for both bosons and fermions. 

Let us now construct the Hamiltonian describing the It is convenient to express the s-wave interaction en- 
rotationally symmetric s-wave interaction under the ap- ergy (Eq. ^) in less opaque nomenclature involving one- 
proximations listed above . By the pseudopotential particle and two-particle spin operators, as we demon- 
method (Huang and Yang 1957), we obtain (Ho, 1998 strate below for the spin-1 and spin-2 cases. 
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1. Spin-1 case 



2. Spin-2 case 



As a starting point, we consider the identity operator 
for a system of two spin-F particles: 



(8) 



all Fr> ; 



Fi • F2 



all F Dai 



^ -^pair V 



(F pair + 1) - F(F + 1) 



F,pai: 



ir.(9) 



where f\ is the one-particle vector spin for particle i. 

We recall that the restriction to s-wave interactions 
allows us to consider only the two-particle spin states that 
are symmetric under the exchange of two particles, i.e. for 
all two-particle spin operators multiplied by the spatial 
(5-function in Eq. ([7]), we may eliminate terms with odd 
Fp a i r from the summations. Restricting to such states, 



we then identify, for F = 1, 



Po and 



Fi • F 2 



P 2 — 2Pq, where the S subscript reminds 



us of the symmetric-state restriction. Combined with 
Eq. 0, we obtain the interaction Hamiltonian of spin-1 
bosons as 





c (1) 






s 2 


Fi • F2 



with 



c (1) 
2 



(i) _ 4ttH 2 a + 2a 2 
C ° = ^ 3 ' 
(i) _ 47r^ 2 a 2 — ao 



c 



M 



(10) 

(11) 
(12) 



In second quantization, each identity operator 
becomes the number-density operator n(r) = 
^Z m =-F V ; m( r )V ; m(r) and the v component of the 
spin vector becomes F v = Em,n=-F ^m( r )(^)mn^n(r), 
where ^ mF (r) is the Bose field operator for the Zeeman 
state labeled by mp- Thus, the second-quantized 
interaction operator is given by 



: Fi • F 2 



, (13) 



where :: denotes normal ordering. The density and spin 
operators are evaluated locally, enforcing the restriction 
to symmetric spin states and allowing us to drop the 
subscript S. 



For F = 2, Eq. (9) gives 



4P 4 - Combined with Eqs. Q 



Fi <8>F 2 
and 



= -6P - 3P 2 

s 

leads to 



y(2) ^0 



(2) 



where we sum over all P pa ir £ {0, 1, 2, • • • , 2P}. Here, 1{ where 
is the identity operator of the ith particle, and (g) denotes 
the tensor product. From the composition law of spin 
operators, we obtain also 



r (2) c r (2) 

^ Fl .F 2 + ^P ,(14) 



(15) 



4ttH 2 



J 4a 2 + 3a 4 a 4 — a 2 7a - 10a 2 + 3a 4 1 
\ 7 ' 7 ' 35 J ' 



The second-quantized form of Eq. (14) is given by 



1 



V® = \j d * [4 2) :n 2 : +cf } :Fi • F 2 : +<%> A^A 2 \{lG) 

where A 2 = Y^L=-2(~ l) m '^m'0-m/v / 5 is the annihila- 
tion operator of a spin-singlet pair of F = 2 atoms. 



C. Pseudo-spin systems 

Any mixture of TV condensates may, in principle, be 
described as a fictitious spin- (N — l)/2 system. For ex- 
ample, a mixture of condensates occupying the two hy- 
perfine st ates \F = l,m z = — 1) and \F = 2,m z = +1) 
of 87 Rb (Myatt et al 1997[ ) constitutes a pseudo-spin- 



1/2 system. Unlike spinor gases, pseudo-spin systems do 
not possess rotational symmetry in spin space, and so 
many of the simplifications attained above for describing 
the interactions and structures of quantum gases do not 
apply. 

Consider a Bose-Einstein condensate composed of two 
internal states of the same atom. Its energetics are de- 
scribed by the following Gross-Pitaevskii (mean-field) en- 
ergy functional: 



E 



dr 



E 

i=l,2 



rfv 2 



2M 



2irh 2 a H , 



M 



47rfi 2 ai2 
M 



*i| 2 |* 2 | 



(17) 



where M is the atomic mass, Ui is the trapping po- 
tential for state \i = 1,2), and an and ai 2 are the 
intra- and inter- species s-wave scattering lengths. To 
draw connection to the characteristics of spinor Bose- 
Einstein condensates, this energy functional can be ex- 
pressed in the pseudo-spin representation as follows. For 
simplicity, we assume U\ = U2 = and decompose the 
condensate pseudo-spinor wavefunction as \E^(r) = 
yn(r)xi(r), where n(r) is the total density of atoms and 
I Xi ( r ) 1 2 + I X2 (r) 1 2 = 1 . Using x% •> we can introduce a unit 
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spin vector S as 
/ S x 

Sy 

where Re and Im denote the real and imaginary parts, 
respectively. Then Eq. (17) can be rewritten as 



/ 2Re( X iX2) 
2Im(xiX5) 
VlXil 



(18) 



1X2 | 



E 



dr 



U ■ 



m^ 2 



h 2 n /^^x9 Mn 9 

m £ (v^) 2 + — v e 2 ff 



-(co + ciS z + c 2 S 2 ) 



(19) 



where v eff = (fi/ m ) Im Z)i=i,2(x*VXi) J and c = 
47rfi 2 (ai + a 2 + 2a 12 )/M, c\ = 47rft 2 (ai - a 2 )/(2M), and 
C2 = 7rh 2 (ai + <22 — 2a\2)jM. This result shows that, 
even at zero applied magnetic field, the system experi- 
ences fictitious interaction-induced linear and quadratic 
Zeeman effects for nonzero c\ and C2, respectively. Differ- 
ing from the case of spinor Bose gases, these terms break 



the SU(2) sy mmetry of the nonlinear sigma model (Ra 

z (S7Si) 2 term 



1987) in which only the £\= 



jaraman 
exists. 

Another generic feature of spinor Bose gases that may 
be absent in the pseudo-spin system is the dynamics of 
interaction-induced spin mixing. For example, for the 
F = 1 spinor gas the s-wave interaction energy (Eq. (p~3|) ) 
contains the terms 



47rh 2 (a2 — do) 



3M 



+1^- 



(20) 

These terms describe spin-mixing collisions wherein two 
atoms in the \m z = 0) state collide to yield atoms in 
the \m z = +1) and \m z = —1) states, and vice versa. 
Such spin mixing is necessarily present if the s-wave in- 
teractions are spin dependent (a 2 ^ «o) 7 and leads to 
characteristic phenomena such as coherent spin oscilla- 
tions, spontaneous magnetization and symmetry break- 
ing, parametric amplification and spin squeezing, as dis- 
cussed in Sec. |IX| Such phenomena are physical conse- 
quences of the coherence between the different internal- 
state components of a Bose-condensed spinor gas. In 
contrast, in pseudo-spin Bose systems, such internal-state 
dynamics may be absent. However, the pseudo-spin Bose 
condensate may still show spontaneous phase coherence 
among its intern al-state populations, as demonstrated 
experimentally by Wheeler et al (2004). In other words, 



a spin-F spinor Bose-Einstein condensate cannot be con- 
sidered an incoherent overlap of N = 2F individual con- 
densates, whereas, in many regards, a pseudo-spin- F con- 
densate can. 



D. Exchange interactions and spin waves 

As argued above, a pair of atoms in a spinor gas will 
only experience a contact interaction if its total spin 



F pair is even. For two atoms in distinguishable mo- 
tional states, this particle-exchange effect leads to an 
interaction-energy splitting between states of even and 
odd 

-^pair? even when the s-wave scattering lengths are 
all equal, similar to exchange-energy effects within multi- 
electron atoms. The exchange- interaction effect is pro- 
nounced for a gas in the cold-collision regime, for which 
all odd partial wave collisions are frozen out. Interactions 
among atoms in a Bose-Einstein condensate are not af- 
fected by this phenomenon, as the particles all occupy the 
same state of motion. Rather, the exchange interaction 
affects the spin dynamics of nondegenerate gases. 

Spin dynamics driven by the exchange effect in the 
cold-collision regime were observed in experiments on 
trapp ed hydrogen ( Bigelow et al.\ |1989 | Johnson et al. , 
|1984[ ). The gas was comprised of atoms with a common 
electronic spin polarization, but differing nuclear spin ori- 
entations. For this particular pseudo-spin 1/2 system, 
because all collisions occur through the electronic-triplet 
molecular potential, all s-wave scattering lengths are ex- 
pected to be the same. However, nuclear magnetic res- 
onances detected in a gas cell placed in an inhomoge- 
neous magnetic field showed narrow resonances ascribed 
to standing spin- wave modes of this rarified gas. 

In a pseudo-spin- 1/2 gas, the exchange interaction 
yields a singlet-triplet splitting that can be understood 
simply as causing the spins of two distinguishable atoms 
to precess about one another. Similarly, in a mean-field 
picture, the spin of an atom passing through a gas that 
is spin polarized in the e direction will precess about e 
at a rate proportional to the product of the gas density 
and the s-wave scattering length. Such an interaction 
leads to a dynamical coupling between spin currents and 
spin polarizations. Solving the gas-kinetic equations in 
the presence of such coupling reveals a spectru m of spin 
waves, wh i ch are da mped by particle diffusion ( |Lhuiller 
and Laloe||l982a|bp . 



The physics of spin waves arose again in studies of 
pseudo-spin 1/2 rubidium gases for which the s-wave 
scattering lengths are, accidentally, all nearly equal. A 
nondegenerate gas was prepared in a superposition of two 
spin states, within a magnetic trapping potential that 
differed slightly for the two states. The slight spin de- 
pendence of the trap led to "anomalously" rapid spin 
dynamics, seen as a spatial segregation of the two spin 
components (Lewandowski et al.\ |2002[ ). Such dynam- 



ics were soon explained as a real-space signature of spin 



waves (Fuchs et al. 


20021 McGuirk et al 


and Levitov| 2002 


Williams et al\ 2002|). 



2002 Oktel 



have a lso been im aged in two-component ultracold Fermi 
gases (Du et al, 2009). Spin- wave dynamics, also called 



the identical spin rotation effect, have been shown to pre- 
vent the dephasing of spins ( "spin- locking" ) in tightly 



confined pseudo-spin- 1/2 g ases (Deutsch et al., 2010 



Kleine Biining et al., 2011| ). The interplay between the 



strong spin-wave dynamics of the thermal fraction and 
the otherwise languid spin dynamics of the condensed 
fraction of a partly condensed pseudo-spin- 1/2 gas was 
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shown t o accelerate the conde nsate spin dynamics signif- 
icantly flMcGuirk erZ||2QQ3 ). 

These works on pseudo-spin gases suggest that ex- 
change effects in nonzero temperature spinor Bose gases 
should be significant, particular in situations where the 
differences among s-wave scattering lengths ai? pa i r are 
small, and also richer than those observed previously 
owing to the higher value of the spin F. Several dis- 
tinct spin-wave modes have been predicted, marked by 
damped oscillations of both the spin vector moments and 
also of the spin quadrupo l e mom ents due to exchange 
effects ( |Nikuni an d Endo, 2008). An examination by 



Natu and Mueller (2010) suggests that spin waves should 
lead to spin-mixing instabilities in nondegenerate F = 1 
spinor Bose gases. However, the full implications of ex- 
change effects in nonzero-temperature spinor Bose gases 
remain unexplored. 



IV. EXPERIMENTAL REALITIES 

A. Optical trapping 

Experiments on spinor gases commence typically with 
the preparation of a spin polarized, and thus scalar, gas 
at low temperature. This spin polarization is established 
in the preparation of the atomic gas, for example, by 
optical pumping and the selective trapping of a single 
spin state in a magnetic trap. The spinor properties of 
the gas become relevant once the atoms are held in a 
trap that accommodates all spin states and the atomic 
spin is no longer constrained to the state of maximum 
magnetization. 

A suitably spin-independent trap is provided by the 
optical dipole force of a focused beam of light. An op- 
tical trap confines an atom by dint of the electric dipole 
moment that is induc ed on the polarizable atom by the 
optical electric field (Metcalf and Straten, |1999 ). The 
trapping potential is expressed as a dyad of dipole vector 
operators, and is, thus, composed generally of a scalar, 
vector, or rank-2 tensor term. In the case of alkali gases, 
for light far detuned from atomic resonance, the rela- 
tive strength of these three contributions scale as 1/A ep , 
Af s /A^ and Ahf s /A^, respectively, where A eg is the de- 



tuning of the trapping light from the atomic resonance, 
Af s is the excited-state fine-structure splitting and Ahf s 
is the excited-state hyperfine-structure splitting. By us- 
ing a detuning A eg that greatly exceeds Af s and Ahf s , the 
scalar light shift can be made to dominate the vector and 
rank-2 tensor shifts 4 . Moreover, by time-reversal symme- 



try, the vector term is eliminated for linearly polarized 
trapping light. For such conditions, the optical trap po- 
tential is regarded as being insensitive to the atomic spin. 
Alternately, by reducing the laser detuning or using el- 
liptically polarized light, the optical dipole potential be- 
comes a tool for manipulating the atomic spin. 

In optical traps, the confinement strength, and hence 
the extent of the trapped gas, along different dimen- 
sions can be varied. An extreme elliptical focus pro- 
duces pancake-shaped gases, with the diameter being far 



small er along one axis than along the others ( |Sadler et al. 
2006). A single-beam red-detuned optical trap with a 



tight circular focus can produce prolate gases with high 



aspect ratios (Schmaljohann et al. , 2004| Stenger et al. 



1998). Forming light traps using two independent inter- 



secting beams allowed for tightly confining and nearly 



isotropic traps (Chang et al., 2004 Scherer et al. 2010). 



Under such confinement, spatial variations of the spin 
and also the density of the degenerate spinor gas can be 
made energetically prohibited along one, two, or three 
dimensions, respectively. Such a dimensional restric- 
tion will apply if the spatial dimension of the gas along 
a given direction is smaller than the healing lengths, 
£sp ~ (87rnAa) -1 / 2 , and £ n ~ (87rna) -1 / 2 , pertinent to 
the spin and density, respectively. Here, we take into 
account just the contact interactions, letting Aa charac- 
terize the relevant scattering-length difference that gives 
rise to the spin-dependent interaction term and a (> 0) 
the scattering length that gives rise to the compressibil- 
ity. Under the typical condition |Aa| < a, the spin dy- 
namics may be restricted to a lower dimensionality than 
the density dynamics. Similar considerations apply to 
determining whether thermal excitations or the MDDI 
can lead to spin variations in each spatial dimension. 



B. Effects of applied fields 

With the atoms in an optical trap, an experiment is 
conducted by setting the initial state of the spin distribu- 
tion and then allowing this state to evolve, either freely or 
under continued perturbation. The chief tool for manip- 
ulating atomic quantum gases, and spinor gases in partic- 
ular, is the application of external electric and magnetic 
fields. Below, we point out some uses and implications 
of these applied fields. 



1. Static magnetic fields: the constraint of spin conservation 
and spin-orbit coupling 



4 For example, for the 87 Rb atom, A^fg is on the order of 10 8 
Hz in the excited 5P states. The fine-structure splitting of 
Afs — 7 x 10 12 is determined from the difference in the reso- 
nant wavelengths of 795 and 780 nm for the Dl and D2 lines, 
respectively. An optical trap with light at a wavelength of 1 /im, 
or A eg = 8 x 10 13 Hz, satisfies the relation A eg ^> (Af s ,Ahf s ) 



The spin-dependent interactions of spinor Bose gases, 
defined by characteristic energies on the order of h x 100 
Hz, become far smaller than the linear Zeeman energy 



and is therefore dominated by the scalar contribution. 



14 



splitting between spin states already at magnetic fields of 
tens of fiG. Under the application of such fields, exother- 
mic dipolar relaxation collisions would allow a spinor 
Bose gas to evolve toward the trivial lowest-energy state 
with uniform magnetization along the field direction. Re- 



cent experiments on Cr spinor gases (Pasquiou et al 



20 l it, for which t he dipolar relaxation rate is substan- 

r 



tial ( jHensler et a/.[|2003||Weinstein et a/.| [l998), confirm 
this evolution toward a single Zeeman state. Only when 
the magnetic field is stabilized to very weak strengths, 
at which either the thermal energy ksT or the spin- 
dependent interaction energy of a condensate becomes 
comparable to the Zeeman energy, does the gas populate 
additional Zeeman states. Studies of non-trivial spinor- 
gas physics in high-magnetic-moment gases will therefore 
require fine magnetic field control and shielding. 



However, for alkali spinor gases (Boesten et al 1996 



Gerton et al. , 1999) (other than cesium (Soding et al 



1998)), dipolar relaxation collisions are rare, and thus, in 
a uniform bias field, the magnetization of the gas along 
the field direction can be considered a conserved quantity. 
In light of the (approximate) conservation of the longi- 
tudinal spin, the theoretical treatments for the ground 
states of the spinor Bose gas must be reconsidered. 

One approach is to introduce spin conservation into 
the many-body Hamiltonian via a Lagrange multiplier A 
( Stenger et a/.|[T998| ). Let the spinor gas in the absence 
of an applied magnetic field be described by the Hamil- 
tonian Hq. Considering also the linear Zeeman energy 
from a uniform magnetic field B = B z z, we substitute 



H Q ^H -(h 1 B z -\)Y J F z , 



(21) 



where 7 is the atomic gyromagnetic ratio, and the sum is 
taken over the atoms in the gas. Let us consider the case 
where F z ^) is constrained to be zero. By symmetry, 
we expect A = hjB z : the effect of the applied field has 
disappeared altogether! Under this global conservation 
constraint, the evolution of the spinor gas magnetization 
is still quite rich, so that meaningful studies of the spinor 
Bose gas are still possible under non-stringent experimen- 
tal conditions. 

The insertion of a Lagrange multiplier is equivalent to 
treating the evolution of the spinor Bose gas in a frame 
rotating at the Larmor precession frequency ujl — —jB z . 
In that rotating frame, the magnetic field appears to 
be absent, and the effects of spin-dependent interactions 
and/or spontaneous symmetry breaking are free to de- 
termine the ground-state or other configurations of the 
quantum fluid. Viewed in the lab frame, however, this 
configuration undergoes rapid Larmor precession, with 
uol typically greatly exceeding other dynamical frequen- 
cies of the gas. This rapid precession must be taken into 
account, for exam ple, in probing the magnetization of the 
gas (Sec. (1X^02 ). 

While the linear Zeeman effect of a uniform magnetic 
field is effectively eliminated through global spin conser- 
vation, that of a non-uniform field B(r) is not. If we 



may neglect spatial variations in the orientation of the 
magnetic field, e.g. in the case that the spinor gas is con- 
fined tightly along some symmetry axis of the magnetic 
field, the constrained Hamiltonian of Eq. (21) is simply 
modified as follows: 



(22) 



The global constraint provides a uniform adjustment to 
the applied magnetic field, but does not cancel out the 
field inhomogeneity. 

This remnant variation was utilized to reveal features 
of the mean-fi eld phase diagram of sodium spinor conden- 
sates ( |Stenger et a/. , 1998[ ). Considering just the gradi- 
ent B' z of the field magnitude along the long axis of their 
quasi-one-dimensional spinor gas and adopting the local 
density approximation, the spinor state at the axial posi- 
tion z was taken to represent the constrained equilibrium 
state for different values oip(z) = X — hjB z (z) = po+p'z, 
with po determined by the longitudinal magnetization of 
the sample. By comparison to the theoretical phase dia- 
gram (Fig. [5j in which one varies also the quadratic Zee- 
man shift; see below), the researchers determined that 
the F = 1 spinor gas of sodium experiences antiferro- 
magnetic interactions. Their approach is among the first 
applications of the now-common method of measuring 
slices through a phase diagram through local measure- 
ments of a gas in an inhomogeneous potential (Ho and 
Horikoshi etlU] [2010] |Navon et al\ |2010| 



Zhou! 2010 



Shin et al 



2008). This method relies on the validity of 



the local-density approximation. For the spinor Bose gas, 
this approximation holds when the locally determined 
equilibrium spin order varies on length scales larger than 
the spin healing length (on the order of microns). 

Accounting for the near-conservation of longitudinal 
spin in a magnetic field with inhomogeneous orientation 
is more subtle. It is convenient to describe the spinor 
wavefunction in a spatially dependent basis of spin states, 
by applying a spatially dependent spin rotation operator 
R(r) = exp[— iQ(r) • F] so that the basis states corre- 
spond to local eigenstates of the spin projection along 
the magnetic field. Applied in this basis, the constrained 
Hamiltonian has the form of Eq. ( 22 ) and also acquires a 



gauge field, with the momentum operator transforming 

as 



(23) 



with the sum taken over components of the spin vector. 
The gauge field is defined through the relation 



Y / ^F fl = ihV r R\ 



#(r) 



(24) 



The imposition of additional rotations about the local 
field direction in R(r) corresponds to a gauge transfor- 
mation of the gauge field A. This gauge field generally 
provides a form of spin-orbit coupling to the spinor gas. 
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As described by Ho and Shenoy (1996), such spin-orbit Santos et aZ. , 2007). This optically induced shift can be 



coupling may introduce vorticity to the ground state of 
a spinor Bose condensate. 

This means of synthesizing a gauge field for a neutral- 
atom gas is related to that realized experimentally by the 
Spielman group ( |Lin et al 2011 , |2009 ). In that work, a 
pair of intersecting laser beams drives a steady Raman 
transition between hyperfine spin states. In a rotating 
frame, this Raman coupling provides a spatially periodic 
effective magnetic field, which then leads to an effective 
spin-dependent gauge field by the mechanism described 
above. 

Even in the absence of dipolar relaxation collisions, the 
magnetization of the spinor gas is generally no longer con- 
served in such inhomogeneous fields. For example, atoms 
in magnetic spherical quadrupole traps may undergo spin 
flips when they cannot adiabatically follow the rapidly 
varying field orientation near the location of zero field. 
In the case of pure magnetic t rapping, such spi n flips are 
known as "Majorana losses," (Majorana 1932) since the 
spin-flipped atoms are expelled from the trap. For an 
optically trapped spinor gas, such Majorana transitions 
would provide a means of Zeeman relaxation. However, 
if the magnetic fields within the volume of the spinor gas 
do not vary rapidly in their orientation, the absence of 
dipolar relaxation collisions will still permit the longitu- 
dinal magnetization - properly defined at each locale - 
to be conserved. 



2. Quadratic shifts 

The Zeeman shift of hyperfine spin states at low mag- 
netic field is not strictly linear. For alkali gases, one can 
evaluate the next-order term in the Zeeman shift given 
the Breit-Rabi Hamiltonian (Breit and Rabi 1931[ ) (Fig. 
[I]). Measuring energies with respect to the \m z = 0) 
state, we obtain the net single-particle quadratic Zeeman 
energy 



Aiy(i + 2/) 2 



BzFz 



(25) 



where g s ~ 2 is the g-factor of the electron, gi^jh is 
the nuclear gyromagnetic ratio, AW is the hyperfine en- 
ergy splitting, B z is the magnitude of the magnetic field, 
assumed to be oriented along z, and the hyperfine spin 
is F = I± 1/2. 

A similar quadratic energy shift can be obtained by 
other means as well. For instance, a quadratic energy 
shift for the F = 1 spinor gas is obtained using the ac 
Zeeman effect from a microwave drive tuned near the 
\F — 1; m z = 0) — \F = 0, 2; m z = 0) hyperfine transi- 



tion (Gerbier et al t 2006). With this method, the sign 



of q can be varied by using either positive or negative 
detuning. Additionally, off-resonance linear polarized 
light will exert a quadratic shift - the tensor portion of 
the ac Stark shift - along the polar ization axis (|Cohen-| 
|Tannoudji and Dupont-Roc , 1972 Jensen et al] 2009; 



spatially tailored at high resolution and rapidly adjusted 
using electro-optics. 

This quadratic shift (which we will call a "Zeeman" 
shift regardless of its provenance) breaks rotational sym- 
metry by favoring one spatial axis, but still preserves 
the lower SO (2) symmetry about that axis. This shift 
gives a significant experimental handle on the ground- 
state phas e diagram of spinor Bose-Einstein condensates 
(Sec. V.A) and on spin-mixing dynamics (Sec. IX). 



3. Effects of radio- frequency, microwave, and Raman coupling 
between spin states 

Magnetic fields are also applied to manipulate the spin 
of the gas dynamically. Transitions among the mag- 
netic sublevels of the spinor gas are driven with rf mag- 
netic fields. Given the strong Zeeman sensitivity of the 
transition frequencies, the radio- frequency (rf) excitation 
is usually either pulsed or frequency-chirped across the 
spin resonance frequency for the purpose of rapid adia- 
batic passage. When the differences between the various 
Am = 1 transition frequencies are small compared to the 
Rabi frequency and the inverse of the pulse duration, e.g. 
at low magnetic fields where the quadratic Zeeman shift 
is small, the rf excitation effects a rotation of the spin. 
When the different Am = 1 transition frequencies are 
resolved, e.g. at high magnetic fields where the quadratic 
Zeeman shift is large, rf fields can be used for arbitrary 



unitary transformations of the atomic spin ( Giorda et al. 
20031. 



The internal state of the spinor gas may also be ma- 
nipulated by driving transitions among the different hy- 
perfine spin manifolds. Here, even at low magnetic fields, 
hyperfine transitions for the different magnetic sublevels 
of the spinor gas are readily resolved. Thus, chosen uni- 
tary transformations may be imposed without applying 
large magnetic fields and contending with the resulting 
eddy currents or hysteresis in the experimental chamber. 

Internal-state transitions of the spinor gas can also be 
achieved optically (Wright et al, 2008). This method 



has the benefit of allowing fine spatial control over the 
coupling between states. A striking use of this capabil- 
ity was the generation of spin-vortex-like structures in 
a sodium spinor condensate using Gauss-Laguerre light 
beams to impose the necessary phase variatio n across the 
gas ( |Leslie et al] |2009a| [Wright et a/.[|2009| ). 



C. Probes 

Here, we highlight measurement techniques tailored for 
studies of magnetism in spinor gases. 
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FIG. 2 Stern-Gerlach and time-of-flight analysis reveals the 
axial spin-state distribution of an elongated 23 Na spinor Bose- 
Einstein condensate, 50 ms after a quench from positive to 
negative quadratic Zeeman shift. The trapped condensate is 
prolate, with radial and axial dimensions of 4 and 270 /mi 
(Thomas-Fermi radii). The optical trap is extinguished, after 
which the condensate expands rapidly in the radial direction 
(vertical in image), so that the in situ axial spin distribu- 
tion (horizontal in image) is nearly preserved. During the 
free expansion, a magnetic field gradient separates the Zee- 
man populations onto separate portions of the image plane 
(a), allowing the column density (b) of each population to 
be quantified. Here, a condensate prepared initially in the 
\m z = 0) state undergoes a dynamical instability, evolving 
toward the transversely aligned polar state (as indicated by 
the correlation between the \m z = ±1) distributions). Fig- 
ure pr epared by Chandra Raman, based on |Bookjans et al.\ 
(2011bL 



gases, the time-of-flight images carry information on the 
spatial, momentum, and energy distribution of atoms in 
the trap ( |Ketterle et al. , 1999). The measurement resolu- 
tion for the in-trap spatial distribution is limited by the 
minimum time of flight needed to separate the spin com- 
ponents, a minimum established not only by the maxi- 
mum available magnetic field gradients, but also by the 
need to allow the gas to expand before applying field 
gradients so as to avoid interparticle scattering as the 
components separate. 

By scattering many photons off each atom being im- 
aged, the populations in the different sublevels can be 
determined with arbitrarily high accuracy. Sub-Poisson 
measurement uncertainty, i.e. with an imprecision AN < 
TV 1 / 2 for measurements o n N atoms, has been achieved 
by several groups flBookjans et al , 2011a Gross et al. 



2Q11| |Hamley et al\ 



2012 Liicke et al, 2011) to probe 



correlations predicted from solutions for the many-body 
spin Hamiltonian (Sec. V.B) and the spin squeezing 



expec ted to result from spin-mixing instabilities (Sec. 
IX.D). 



2. Dispersive optical measurements 

To achieve higher spatial resolution while imaging the 
trapped and usually optically dense gas, it is advanta- 
geous to turn to dispersive imaging using off-resonant 
light, for which the refractive index of the cold gas is 



dominantly real ( Andrews et al. 1996). The dispersion 



of light passing through the sample conveys information 
not only on the density of the atomic gas, but also on 
its spin polarization. The real part of the linear optical 
susceptibility \ (which is unitless) for monochromatic il- 
lumination is a tensor quantity related to the reduced 
one-body, local density matrix p(r) as 



1. Stern-Gerlach and time-of-flight analysis 



■di\e){e\d s \g)(g\ (26) 



Characterizing the magnetic order in a spinor gas re- 
quires knowledge of both the populations in and the 
coherences among basis states of the atomic spin. A 
straightforward method for measuring populations is re- 
leasing the gas from its optical confinement, using a mag- 
netic field gradient to separate atoms in different mag- 
netic sublevels into different spatial regions, and then 
imaging the separated portions of the gas; this is an amal- 



gam of the Stern-Gerlach experiment (Gerlach and Stern 
1924) and time-of-flight analysis (Fig. |2|. In typical ex- 



periments, the magnetic field is varied adiabatically be- 
tween the uniform field applied during optical trapping to 
the gradient field applied during the Stern-Gerlach sep- 
aration; thus, one obtains a projective measurement of 
the distribution of atoms among the Zeeman sublevels in 
the trap. 

Additionally, as in measurements of scalar quantum 



where ij denote Cartesian directions, d is the electric 
dipole operator, and eo is the vacuum permittivity. Here 
we assume the population is entirely in the electronic 
ground state manifold. The sum runs over the energy 
eigenstates of the ground-state manifold (labeled with g) 
and those of the electronic excited states (labeled with 
e), with A eg = uo — u eg indicating the optical detuning 
for each transition. 

To simplify this expression further, we assume the ex- 
cited states are divided into several manifolds of well 
defined angular momentum, and make the approxima- 
tion that the detuning for transitions between the ground 
states (taken as a single spin manifold) and each of the 
excited-state manifolds is independent of the magnetic 
quantum numbers of each. This approximation is well 
suited to spinor gases, studied under conditions of small 
Zeeman and Stark shifts compared to the optical detun- 
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ing. The susceptibility now becomes 



Xij(r) 



Aire, 



-Tr p(r)di 



\e)(e\ 

hA P „ 



(27) 



where the operator in square brackets is scalar. As the 
expectation value of a dyadic tensor, the optical suscep- 
tibility is identified as a sum of spherical tensors of ranks 
0, 1 and 2, i.e. 

Xij (r) = x (0) ^nW+X (1) ^^W + X ( %W. (28) 

Thus, the scalar, rank-1 and rank-2 components of the 
optical susceptibility tensor relate directly to the den- 
sity n(r), the vector magnetization F(r) and nematicity 
Afij(r), respectively, of the spinor gas. The vector and 
tensor contributions give rise to circular and linear bire- 
fringence, respectively. This optical sensitivity summa- 
rizes magneto-optical properties of atomic gases ( |Budker 



et al. ,2002 Suter 2005 [ ), and was expressed in relation t o 
studies of spinor gases by Carusotto and Mueller ( 2004 ) . 



One means of measuring the susceptibility tensor is to 
send collimated monochromatic probe light through the 
spinor gas. In passing through a thin slice of the gas, 
the electric field of probe light propagating along the z 
direction is modified as 



E(x,y,z+dz) = 



1 + ik ( 1 + y, z) I dz 



-E(x,y,z) 



(29) 

where k is the wave number in vacuum. The electric 
field lies in the x-y plane, and only the x, y elements 
of the susceptibility tensor are retained in x±, owing to 
the transverse nature of the optical field. When the dis- 
persive phase shifts are small, and applying the thin- lens 
approximation, the field of the light after passing through 
the entire gas is approximated as 



E ont (x,y) 



-k 



(tfxx(x, y) - Myy(x, ?/)) i X W F z (x, y) + X {2) ^xy{x, y) 

-iX {1) Fz(x,y) J cX {2) N xy (x,y) (Nxx{x,y) - ti yy {x,y) 



(30) 



where the common phase shift is <j) = (ko/2)[x^fi + 
X^ 2 \Nxx +«A^/i/)/2] and the tilde denotes the column av- 
eraged quantity along the optical axis. 

Dispersive optical detection has been used to mea- 
sure the magnetizati o n dyna mics of F = 1 spinor gases 
of 23 Na ( |Liu et al] |2009b| ) and 87 Rb ( |Higbie et al 



2005). The sodium experiments highlighted the capacity 



for performing weak continuous measurements on a sin- 
gle trapped spinor gas. Strong optical confinement con- 
strained the sodium spinor condensate to a single spatial 
mode. Far detuned linear polarized light was passed once 
through the sample and then through a linear polarizer 
rotated 7r/4 with respect to the incident polarization in 
order to determine the circular birefringence. Through 
Eq. (30) and consistent with the stated assumptions, the 



measured signal is proportional to fcoX^^z? giving the 
magnetization component of the gas along the optical 
axis. By imposing a magnetic field transverse to this 
axis, the rapid Larmor precession of the magnetization al- 
lowed both components of the transverse magnetization 
to be detected as a temporal oscillation of the optical 
rotation signal. The amplitude of this oscillation (Fig. 
[3| gave a continuous measure o f spin mixing dyna mics 



in the polar spinor condensate (Liu et al. 2009b]). In 



experiments on nondegenerate gases, it has been demon 
strated that a continuous Faraday rotation measurement, 
supplemented with imposed linear and quadratic Zeeman 
effects, can be used for complete character ization of th e 
spin density matrix for arbitrary spin (Smith et al. , 2004). 



This powerful method may be extended to the study of 
spinor gases in future works. 



Applications of dispersive birefringent imaging to ru- 
bidium gases have highlighted the ability to measure the 
vector magnetization at high spatial resolution. Circu- 
lar polarized light was sent through the gas (along the 
y axis), and a polariz ation-independent phase-contrast 
method (Hecht, 1989) was used to convert the disper 



sive signal into an intensity image (Higbie et 



Following Eq. (30|), the signal is proportional to x 

x (1) ^ + x (2) (M 



2005). 

-M zz )/2. To isolate the magnetiza- 
tion signal, a series of images was taken while the atomic 
spins underwent Larmor precession in a transverse ap- 
plied magnetic field (along z), with each probe pulse last- 
ing a small fraction of the Larmor precession period. The 
pixel-by-pixel amplitude and phase of the signal varia- 
tion at the Larmor precession frequency (aliased by the 
slower frequency at which images were taken) gave a de- 
tailed spatial map of both the magnitude and orientation 
of the transverse magnetization (Fig. [4]). In later work, 
a spatial map of the longitudinal magnetization was ob- 
tained also by applyi ng a tt/2 rf pulse in the middle of 
the imaging sequence ( Vengalattore et al. , 2010). The ne- 



maticity signal, precessing at twice the Larmor frequency 
and with x^ being small for the imaging settings, was 
ignored. 
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FIG. 3 (Color) Faraday rotation measurement, revealing 
spin-mixing dynamics in a 23 Na spinor Bose condensate, (a) 
Magnetization along the probe light axis causes the optical 
linear polarization to rotate, as measured by a photocurrent 
difference in a polarization analyzer. With the condensate ex- 
posed to a magnetic field perpendicular to the imaging axis, 
the rotation signal measures the transverse magnetization, 
which oscillates at the Larmor precession frequency. (b,c) 
The Larmor precession amplitude is recorded for a conden- 
sate that is prepared in the transversely magnetized state 
and allowed to evolve at constant q. (b) For low \q\, at- 
tained at low magnetic field, the spin-dependent interactions 
cause the spin-mixing oscillations to follow an orbital trajec- 
tory around the transversely magnetized state; the condensate 
remains mostly magnetized, and, thus, the Larmor precession 
amplitude varies only slightly, (c) In contrast, for high \q\, 
the condensate undergoes complete orientation-to-alignment 
conversion as expected in the absence of interactions. This 
is indicated by the periodic oscillation between large preces- 
sion amplitude, characteristic of the oriented (magnetized) 
state, and zero precession amplitude, characteristic of the 
aligned (unmagnet ized) state. Figu re prepared by Yingmei 
Liu, adapted from Liu et al\ ( |2009b| ). 



3. Fundamental limits to dispersive measurements 

The back action of optical measurements of the atomic 
spin influences the gas being probed, necessitating a 
compromise between measurement precision and the fi- 
delity of repeated measurements. When measurements 
are made of several spin moments of the gas, the lim- 
itations to measurement fidelity are fundamental, aris- 
ing from the non-commutativity of the measured observ- 
ables. These limitations have been studied extensively 
for measurements of the vector spin for gases of spin- 1/2 
atoms (|Takano et al. 2009 ) and also higher-spin atoms 
( Kuzmich et al. , 1998, 2000 |, and for pseudo-spin- 1/2 sys- 
tems (lAppel et al. 1 12009; Esteve et a/.[|2008[|Gross et al. 



201Q||Hammerer et a/.||2010||Riedel et al4|2010||Schleier- 



Smith et al. 2010). For measurements on a spinor gas 



that is fully magnetized along the e axis, standard quan- 
tum limits apply to measurements of the two magnetiza- 
tion projections transverse to e, defining the condition 



when the requisite back action from measurements of 
one projection contributes significant fluctuations to the 
other. This limit to the measurement variance is given 
as (AMi) 2 = nVF/2 where n is the atomic density, F 
is the (dimensionless) quantum number of the total an- 
gular momentum, and V is the volume over which the 
magnetization is measured. 

Dispersive measurements on optically thick gases are 
capable, in principle, of surpassing the standard quantum 
limits for spin detection non-destructively. Dispersive 
measurements on spinor gases have yet to achieve this 
precision, being limited by c ollective optical scatterin g 
from the optically dense gas ( Vengalattore et al, 2007). 



However, the prospect of doing so is compelling. Analyz- 
ing the fluctuations in a low-noise dispersive optical mea- 
surement may distinguish between the many-body cor- 



related spin states of spinor gases (Eckert et al., 2007b). 



Similar to the cases cited above for spin- and pseudo-spin- 
1/2 atoms, quantum non-demolition measurements of the 
atomic spin provide a means to spin squeezing. Achieving 
such measurements in high-resolution dispersive imag- 
ing would allow one to transform a spinor Bose-Einstein 
condensat e into a non-class ical, spatially extended quan- 
tum field (Sau et al., 2010). It is interesting to consider 



how dispersive imaging, or other non-destructive detec- 
tion methods, can achieve squeezing and entanglement 
that is specific to ensembles of high-spin particles, cri- 
teria for which have been specified by Vitagli ano et al.\ 
(|2011D (see also Sec.|IX.D[). 



4. Cavity aided detection and manipulation 

Placing a dispersive medium inside an optical res- 
onator leads to an enhancement of the dispersive optical 
phase shift by the cavity finesse. The cavity- enhanced de- 
tection of the magnetization of a cold atomic gas enabled 
a recent implementation of magnetic resonance imag- 
ing for optically trapped ultracold atomic gases. The 
atomic magnetization was measured with an impreci- 
sion below the Poisson limit, although not as a quantum 
non-demolition measurem ent due to pe rturbations to the 
atomic motional states ( jBrahms et a/.|[2011[ ). This sen- 
sitivity is sufficient to detect quantum- correlated ground 
states of spinor Bose condensates according to recent pro- 
posals flCui et al] |2008| |Zhang et al] |2009[) . 



Optical cavities have also been examined as a means to 
influence the dynamics of spinor Bose gases. Variations 
in the atomic spin populations alter the optical suscep- 
tibility of the gas, changing the resonance frequency of 
differently polarized optical cavity modes. In turn, for 
an externally driven cavity, these changes modulate the 
cavity fields and thereby act back on the atomic spin dy- 
namics through the varying ac Stark shifts. This coher- 
ent back action is predicted to lead to nonlinearities and 
bist ability in spin- mixing dynamics ( Zhou et al] 2010| 
[2009] ). 

Missing from these theoretical descriptions is a full ac- 
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FIG. 4 (Color) Repeated polarization-contrast imaging is used to measure both components of the transverse magnetization 
of a spin- 1 87 Rb gas. (a) Several brief polarization-rotation images are taken while the transverse magnetization undergoes 
Larmor precession between image frames, (b) Such precession is seen in data selected from a common location in each image 
frame (yellow boxes in a), (c) Analyzing data from each image pixel, one extracts the spatially varying Larmor precession 
amplitude and phase, or, equivalently, the two components of the transverse magnetization, F x and F y . These are represented 
either in gray scale, or in a color scale with brightness indicating amplitude and hue indicating orientat ion (c olor wheel shown 
m Fig(6}). This image analyzes a degenerate gas equilibrating to a ferromagnetic state ( |Guzman et al.\\2011\ . 



counting for the effects of quantum noise and quantum 
measurement. Aside from influencing the atomic spin 
dynamics, the cavity field also effects a continuous mea- 
surement of the atomic spin. Such measurement requires 
a back action on the atomic spin, which comes typically 
from the quantum-optical (shot noise) fluctuations of the 
cavity field. Some aspects of this back-action, in terms 
of disrupting the otherwise coherent spin-mixing dynam- 
ics, are discussed by | Zhang et al. (2011). Cavity quan- 
tum electrodynamics offers a clean setting in which to 
explore such measurement back action effects in future 
experiments. 



given by 



pF z + qFl 



^m 2 (r), (32) 



where we have assumed that the magnetic field is oriented 
uniformly along z, and that the quadratic Zeeman shift 
also selects the z axis. 



V. GROUND-STATE PROPERTIES OF SPINOR BOSE 
GASES 

A. Mean-field ground states 

Gathering energy terms discussed in the previous sec- 
tions, and continuing to ignore the MDDI (see later dis- 
cussion in Sec. VIII. C), the total Hamiltonian of the 



spinor Bose gas is given by 



h 2 v 2 

' 2M 



U(r) 



+H Z + V, 



(31) 



where U is the trapping potential, assumed to be scalar, 
and Hz describes the linear and quadratic Zeeman effects 



To focus only on the spin-dependent energies, let us 



make the single-mode approximation (SMA) (|Law et al. 



1998) whereby we express field operators as ^(r) = 
< I ) (r)'0 m , with 3>(r) defining the spatial mode of the Bose- 
Einstein condensate. This approximation is valid so long 
as the energetic cost of spatial variations in the spin state 
of the gas exceed any consequent reductions in the spin- 
dependent energies, or equivalently, if the system's size 
is much smaller than the spin healing length £ sp . By 
the discussion of Sec. |IV.A[ such a condensate is zero- 
dimensional for spin degrees of freedom. When the sys- 
tem's size is larger than £ sp , the SMA may also be ap- 
plied locally within the local density approximation. The 
ground states obtained under the SMA can thus clarify 
how spin-dependent interactions are locally minimized, 
even in spatially extended samples. 
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1. Spin-1 case 

Within the SMA, the spin-dependent part of the 
Hamiltonian for a spin-1 Bose gas reads as 

(l) 

^ = - ^ ipl ni tpl n2 F mirn3 • F m2rn4 ip m3 ip rn4 



mi ,m2 

7713,7714 



of an F = 1 87 Rb condensate were measured after the 
gas was prepared with zero longitudinal magnetization 
and allowed to equilibrate for several seconds at variable 
q. At high g, the condensate equilibrated with its entire 
population in the \m z = 0) state, while in the range 
> q > go, this fractional population diminished to 
about one half, as predicted by the mean-field treatment 



(33) 



where n = J d 3 r |$ 4 |/ J d 3 r |<£ 2 | oc N is the average den- 
sity, proportional to the number of atoms in the gas N. 

If we ignore quantum fluctuations, and take expecta- 
tion values with respect to the mean values of the mode 
operators, we obtain the mean-field energy functional 



( Chang et al] 2004]). A second experiment provided ev- 
idence that this state of mixed Zeeman populations was 
indeed transversely magnetized. For this, the equilibrium 
state at low q was coherently manipulated by briefly puls- 
ing a higher value of g, launching spin mixing dynamics, 



as described in Sec. [IXj Such dynamics depend on the 
coherence between Zeeman populations, and confirmed 



the m agnetized nature of the initial state (Chang et al. 
2005b. 



(i) 



(F) 2 +p(F z )+q(F*). 



(34) 



The spin-dependent interaction term is rotationally 
symmetric. For > 0, as for the F = 1 spinor gas 
of 23 Na, this term favors the manifold of polar states, for 

which the (F) 



The transition between easy-axis and easy-plane ferro- 
magnetism at q = was measured by the Berkeley group, 
who applied in-situ magnetization-sensitive imaging to 
probe the spin structures formed in 87 Rb spinor conden- 
sates after several seconds of equilibration (Fig. |6|. The 
observed textures contained large commonly magnetized 
spin domains with the predicted spin-space anisotropy 



= 0. Such interactions are often denoted (Guzman et al, 2011) 



as "antiferromagnetic," though unlike solid-state antifer- 
romagnets, there is no Neel order. For < 0, pertinent 
to the F = 1 spinor gas of 8T Rb, this term favors the man- 
ifold of ferromagnetic states, for which (F) = 1. Both 
manifolds of states are inert because they depend only 
on the sign and not on the magnitude of c{\ 

The application of external fields, described by the ef- 
fective linear and quadratic Zeeman energies p and g, 
respectively, breaks rotational symmetry. The minimum- 
energy states of are represented in Fig. [5] 

Let us consider the case of p — 0, on which many 
experiments have focused. Three classes of inert states 
attain extremal values of the quadratic Zeeman term, 
and are thus identified as extrema of E^: the longitu- 
dinal ferromagnetic states ip = (1,0, 0) T and (0,0, 1) T ; 
the longitudinal polar state ip = (0, 1, 0) T , and the trans- 
verse polar states ip = (1, 0, e lc ^) T /\/2 for arbitrary (p. 
In the regime \q\ < qo = |2c^n|, we identify an addi- 
tional, non-inert energy extremum of partial transverse 
magnetization that arises from the competition between 
the spin-dependent interaction and quadratic Zeeman ef- 
fect. The energies of these extremal states are evaluated 
in Fig. [5] 

For ferromagnetic coupling, the non-inert transversely 
magnetized state emerges as the lowest energy state at 
q = qo, marking the point of a second-order phase tran- 
sition that breaks axial symmetry. At q — 0, we observe 
a transition between easy-plane (q > 0) and easy-axis 
(q < 0) ferromagnetism. 

Both transitions have been observed experimentally. 
The onset of transverse ferromagnetism at q = go was 
observed by the Chapman group ( Chang et al\ |2004| 
2005). In a first experiment, the Zeeman populations 



The antiferromagnetic interactions of F = 1 23 Na con- 
densates were identified by allowing such condensates 
to equilibrate in the presence of a magnetic field gra- 
dient. Stenger et al\ (1998) observed that, for suffi- 
ciently large q > 0, the spatial extent of the \m z = 0) 
population of the condensate was larger than predicted 
for spin-independent interactions, demonstrating that 
spin-dependent interactions make the longitudinal polar 
state lower in energy than the longitudinal ferromagnetic 
states. 

A second signature of antiferromagnetic interactions is 
the immiscibility o f the \m z = 0) and \ m z = —1) (or 



+1)) states ( |Miesner et al. 1999 Stenger et al. 



1998). This immiscibility (Fig. |22|) results from the rela- 
tion aij > y/auajj between the self- and cross-scattering 
s-wave scattering lengths for two equal-mass distinguish- 
able atomic states labeled i and j. For the F = 1 spinor 
system, a_i_i = ao-i = <22 and aoo = (2a2 + ao)/3. 
For antiferromagnetic interactions, ao > &2 and the im- 
miscibility condition holds; in the diagrams of Fig.[5j this 
immiscibility is reflected in the sharp transitions between 
the \m z = ±1) and \m z = 0) phases (for q > 0). 

Along the p — axis, the antiferromagnetic F = 1 
spinor condensate should always favor a polar ground 
state. Its alignment is predicted to vary from longitu- 
dinal, for q > 0, to tran sverse, for q < 0, as h as been 
observed experimentally QBookjans et al 2011b). 



2. Spin-2 case 

As in the spin-1 case, we assume that the system is 
spatially uniform and ignore the kinetic energy. Within 
the SMA and neglecting the MDDI, we can write the 
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FIG. 5 (Color online) Mean-field phase diagram for F — 1 spinor Bose condensates with (a,c) antiferromagnetic or (b,d) 
ferromagne tic s-wave in teractions. (a,c) Ground-states for variable linear (p) and quadratic Zeeman shifts (q), as derived 
by Stenger et al. (1998). For the antiferromagnetic case, the solid curve parameterized by p = ±^/2q (0 < q < 1/2) and 
p — ±(g + 1/2) (q > 1/2) marks the sudden (first-order) transition from the longitudinal polar state (\m z = 0)) to either 
one (\p\ > 1) or a superposition of both (\p\ < 1, shaded region) longitudinal ferromagnetic states. The dashed lines indicate 
smooth (second-order) transitions in state composition. For the ferromagnetic case, the dashed lines p = zLq (q > 0) and 
p — ±y/q(q — 2) (q > 2) bound the regions where the ground state transforms smoothly between the longitudinal polar and 
the longitudinal ferromagnetic states. (c,d) Extrema of the mean-field energy functional for p — 0, plotted for three inert states 
- the longitudinal ferromagnetic (red dashed line) , the longitudinal polar (blue dotted) , and the transverse polar states (green 
dot-dashed) - and the non-inert extremal state for \q\ < 2 (black). The antiferromagnetic ground state is either a longitudinal 
or transverse polar state for positive or negative q, respectively. The ferromagnetic ground state is either the longitudinal 
polar state (q > 2), the partially magnetized state (0 < q < 2) or the longitudinal ferromagnetic state (q < 0). Majorana 
representations are shown. In quench experiments, a condensate is prepared in the longitudinal polar state at low q, an extremal 
high energy state that is subject to dynamical instabilities. All energies scaled by |c^n|. 



spin-dependent part of the mean-field energy functional 



and KoasrnT2002) 



as (Ciobanu et a/., 2000 Koashi and Ueda, 2000 Ueda 



,( 2 ) 



n 



(Fy + ^—\(A)\^p(F z )+ q (Ft). (35) 



Compared with the spin-1 case, the energy 
functional involves an additional term A2 = 
Y^m=-2(~^) m ^rn4 ) -rn/V5 which is the annihilation 
operator of a spin-singlet pair of F = 2 atoms. 

The mean-field phase diagram of a spin-2 Bose- 



Einstein condensate is summarized in Table III and Fig. 



[7| Consider first the ground state phases in the ab- 

(2) 

sence of the quadratic Zeeman energy. When c\ ' < 

(2) (2) 

and c 2 > 4c\ , the interaction energy is minimized 
by maximizing the magnetization and minimizing the 

spin-singlet amplitude, so that the ground state is fer- 

(2) (2) (2) 

romagnetic. When c\ < and c 2 < 4c^ , the in- 
teraction energy is minimized by maximizing the spin- 
singlet amplitude and minimizing the magnetization, so 
that the ground state is antiferromagnetic with the spinor 
given by (0, 0, 1, 0, 0) T (and the family of states generated 
therefrom by rotations). It is not appropriate to call this 
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q/h= -5 Hz 



q/h=0 Hz 



q/h=5 Hz 




Time 0.25 0.50 0.75 1.0 1.5 2.0 s 



FIG. 6 (Color) Transverse (top) and longitudinal (bottom) magnetization of an F — 1 87 Rb spinor Bose gas that was prepared 
initially in a nondegenerate incoherent (1/3, 1/3, 1/3) population mixture of the Zeeman states and then cooled to quantum 
degeneracy. The degenerate spin texture evolves for variable time at a quadratic shift of q/h = -5 Hz (a), Hz (b), and 5 Hz 
(c) before the column density of the vector magnetization is measured. The transverse magnetization is represented using the 
color wheel shown, where the hue indicates the magnetization orientation and the brightness its magnitude. The longitudinal 
magnetization is represented by the color bar. The maximum brightness is set by a fully polarized condensate at 2 s evolution 
time. Textures form small domains of commonly oriented magnetization and then coarsen to where only a few domains span the 
condensate. The late-time images show a spin-space anisotropy; for positive q the transverse magnetization appears brighter 
while for negative q the opposite trend is o bserved. The condensates contain ~ 10 6 atoms and are two-dimensional with respect 
to spin dynamics. Figure reproduced from Guzman et al. (2011). 



state polar because it has no polarity. In fact, the or- 
der parameter of this phase is proportional to the m = 
component of the spherical harmonic function of rank 2: 
y 2 °(^7 <j>) 3 cos 2 — 1 which is even under — » + it. 
In contrast, the spin-1 polar state (0, 1,0) T does feature 
polarity, changing sign upon a rotation by tt about a 
transverse axis. 

The order parameter for the nematic phase is given 
by (lUeda and Koashil 120021) 



therefrom by rotations. This state, called either the 
cyclic or the tetrahedral state, has the symmetry of a 
tetrahedron under rotations and inversions. As in the 
F = 1 polar state, these discrete geometric operations 
must be accompanied by gauge transformations to leave 
the state unchanged; a fact that we discuss further in Sec. 



sin r] sin r\ 

, 0, cos 77, 0, 



V2 



V2 



(36) 



Depending on the value of 77, this order parameter de- 
scribes the uniaxial nematic (UN) phase (77 = 0), the 
biaxial nematic (BN ) phase (77 = 7r/6), and the dihedral- 
2 phase (otherwise) ( jUchino et al. 2010). Although these 



VII The cyclic/tetrahedral phase has neither spin-singlet 
amplitude nor magnetization, but it breaks time-reversal 
symmetry. 

In the presence of an external magnetic field, the 
quadratic Zeeman term lifts the degeneracy between the 
UN and BN phases in favor of the UN (BN) phase 
for q > (q < 0). The quadratic Zeeman effect 
also lifts the degeneracy of the two cyclic/tetrahedral 
phase described above, distorting these states into either 
C c = (sin 0/V%O,-zcos (9,0, sin 6/y/2) with cos 2 = 

l/2+5<?/(4 2) n) or C M = (cos (9, 0, 0, sin (9, 0) with cos 2 = 
(2) 

1/3 — g/(3c\ ; n), in different portions of the phase dia- 



gram (Saito and Ueda, 2005). These phases retain ei- 



states are degenerate at the mean-field level, zero-point 
fluctuations in the Bogoliubov modes lift the degeneracy, 
causing a first-order (order-by-disorder) quantum phase 

transition between the UN phase w i th > and the ther the original three-fold (( M ) or four- fold (£ c ) sym- 
metries, but they are distorted along the z direction by 
the quadratic Zeeman effect (see Fig. [7]). The phase ( c 



BN p hase with ) < (|Song et al 



2007 



Turner et al. 



2007). It is interesting to note that the parameter rj char- 



acterizes the geometric shape of the order parameter, as 
seen, for example, in the Majorana representations of the 
various states in Fig. [7[ However, it is not related to a 
symmetry of the Hamiltonian, and thus does not lead to 
a conserved current, even though the eigenspectrum is 

degenerate with respect to rj. 

(2) (2) 
When c\ > and c 2 > 0, neither ferromag- 
netic nor antiferromagnetic order is energetically favor- 
able. Rather, the zero-field ground state is the state 
(1/2, 0, — i/y/2, 0, 1/2) T and the set of states generated 



remains nonmagnetic, whereas £ M acquires magnetiza- 



tion proportional to 2 cos 2 — sin 



-g/(3c^n); the 



"mixed" phase is named after such magnetic hybridiza- 
tion. 

At present, the only experimentally realized spin 
2 spi nor Bose-Einste i n con d ensate is a 87 Rb c o nden- 
sate (Chang et al] 



Schmaljohann et al] 



2004 
2004^ . 



Kuwamoto et al. 2004 



Although this state is in 



the upper hyperfine manifold, the lifetime is as long as 
a few hundreds of milliseconds. A theoretical calcu- 
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Phase 


Order Parameter 


Parameter Region 


F ( n <r (\\ 


V 1 ? U 5 U 5 U 5 U J 


r (2) > 4 _(2) (2) |g| 


"Rivr ^ ^ ^^ 
olN < UJ 


f i n n n 1 \ 


(2) . A (2) (2) . 2|g| 
u 2 ^ 5 u 2 ^ n 


M < (J J 


(cos c/M, 0, U, sin c/m, UJ 


(2) A (2) (2) . |q| 


c 






BA (q > 0) 


(±a, 6, c, 6, ±a) 


+ (-) sign for cf } < (> 0) 


UN (q > 0) 


(0,0,1,0,0) 


numerically determined 



TABLE III Definitions of the various states and those of the parameter regions in which they are expected to be the ground 
state, where F, BN, M, C, BA, and UN stand for ferromagnetic, biaxial nematic, mixed, cyclic, broken-axisymmetric, and 

uniaxial nematic phases, respectively, q is the quadratic Zeeman coefficient, and cos 6m = y/l/3-4/(3(4 2) n)and cos#c — 



l/2 + 5<?/(4 2) n) (Saito and Ueda 



are determined numerically ( |Uchino et aT. 



2005). The phase boundary between the BA and UN phases and the components a, b, c 



2010). 



lation of the interaction parameters predicts that this 
condensate lies slightly on the antiferromagnetic side 
of the phase boundary between antiferromagnetic and 



cyclic/tetrahedral phases, i.e. that c 



(2) 



(2) 

> and c\ ' is 



only slightly positive (Klausen et al. 2001). 

TT] 

The coupling constant c\ characterizes the spin- 

(2) 

mixing processes with Amp = ±1, while c 2 character- 



izes the processes with Am^ = ±2. Schmaljohann et al 
( 2004 ) observed the time evolution of the F = 2 conden- 



sate by starting with the condensate with all atoms in 
the \m z = 0) state. The subsequent evolution shows that 
the m z = ±1 components first appear, followed by the 
appearance of the m z = ±2 components, indicating that 



is significantly larger than c^ ; . They also investi- 
gated the stability of various initial states against spin- 
mixing collisions, and found that the antiferromagnetic 
configuration \m z = 2) + \m z = —2) is stable and that 
the cyclic/tetrahedral configuration shows nearly no spin 
dynamics, suggesting that the ground state is slightly on 
the antiferromagnetic side of the phase boundary. How- 
ever, it should be noted that their experiments were per- 
formed in the presence of the quadratic Zeeman effect 
of the order of a few nK which may well stabilize the 
antiferromagnetic configuration. 

By placing a pair of atoms in each site of an optical 
lattice and letting them undergo spin-mixing oscillations, 
Widera et al ( 2006 ) were able to determine the values of 
scattering length difference within the 87 Rb F = 1 and 
F = 2 manifolds. They suggested that "the data points 
toward the antiferromagnetic ground state" for the F = 
2 state with the reservation that the cyclic/tetrahedral 
ground state cannot be excluded because of the error 
bars. 



(2) 



B. Exact ground states and fragmentation 

In the previous discussion, we assume a priori that 
the spinor Bose gas forms a Bose-Einstein condensate 



with a single macroscopically occupied single-particle 
ground state, described by the mean-field condensate 
spinor wavefunction. However, in the presence of spin- 
dependent interactions, is it necessarily the case that such 
a mean-field, coherent spin-state condensate correctly de- 
scribes the gas? 

In the theory of scalar Bose-Einstein condensates, 
th e validity of the me an-field assumption was quantified 
by Bogoliubov (1947). Due to interparticle scattering, 



the many-body ground state of the Bose gas includes a 
non-zero population in the excited single-particle states. 
However, this excited state population is kept small be- 
cause of the additional single-particle energy required to 
occupy such states. The fraction of atoms found outside 
the single-particle ground state, called the quantum de- 
pletion, remains small as long as the typical excited-state 
energy (scaling with density as n 2 / 3 ) is much larger than 
the typical interaction energy (scaling as no 2 ). 

However, for spinor Bose-Einstein condensates, the en- 
ergetic barrier against the interaction-induced mixing of 
different many-body spin states is absent. For exam- 
ple, the spin-1 many-body Hamiltonian contains the spin 
mixing term given in Eq. ( |20| ). The mean-field state of 
an TV-atom, polar F = 1 condensate aligned with the 
z axis, written as |0,7V, 0) in the Fock-state F z eigenba- 
sis, couples by spin mixing with the state |l,iV — 2, 1). 
The latter state overlaps with the polar mean-field state 
canted slightly from the z axis, a state that, for q = 0, is 
equal in energy to the initial state. Thus, spin mixing is 
seen to couple degenerate coherent spin states. Following 
the arguments of Bogoliubov regarding the scalar Bose 
gas, one might therefore expect the quantum depletion 
due to such spin mixing to be massive. 

Wi thout making the a priori mean-field approxima- 
tion, Law et al\ (1998) showed that, within the SMA, 



the many-body Hamiltonian of a spin-1 condensate can 
be diagonalized exactly. The spin-dependent interac- 
tions are symmetric both under rotation and under par- 
ticle exchange; thus, naturally, the many-body s-wave 
spin-dependent scattering Hamiltonian can be written in 
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FIG. 7 (Color online) Phase diagrams of spin-2 Bose-Einstein 
condensates in the presence of the quadratic Zeemen energy q. 
(a) Case of q < 0, where BN, C, F, and M stand for the biax- 
ial nematic, cyclic, ferromagnetic, and mixed phases, respec- 
tively. In the absence of an external magnetic field (q = 0), 
the C and M states can be transformed to each other by a 
rotation in space, (b) Case of q > 0, where the UN and 
BA stand for the uniaxial nematic and broken-axisymmetric 
phases, respectively. The broken-axisymmetric phase has a 
transverse magnetization and a spin-singlet pair amplitude, 
both of which depend on q, as in the case of the spin-1 broken- 
axisymmetric phase. The thick lines show the phase bound- 
aries. The Majorana representation is inserted in each phase. 
Figure adapted from|Uchino et al. (2010). 



terms of the total spin operator, F to t = 5^ F^, as 



H 



2 I N 



(37) 



The ground states of this Hamiltonian can be read off 
immediately. An Af-atom ferromagnetic F = 1 conden- 
sate has a ground-state manifold described by the states 
of maximum F to t = N. A constraint on the total longitu- 



dinal spin, e.g. F^tot = Mf? is accommodated by select- 
ing the corresponding magnetic sublevel of the F to t = N 
manifold. For example, in the case of strictly zero lon- 
gitudinal magnetization, the exact ferromagnetic ground 
state is |F to t = N,M F — 0). This state differs from the 
transversely oriented mean-field state by the exact equal- 
ity between the \m z = +1) and \m z = — 1) single-particle 
Zeeman state populations. 

The ground state of an antiferromagnetic F = 1 con- 
densate is more subtle. At zero magnetization, the 
ground state, in the case that N is even, is the unique 
state with F tot — 0. Unlike the mean- field ground 
states, this many-body state retains rotational symme- 
try. The character of this state is clarified by con- 
structing the many-body spin state using the two-particle 
spin-singlet creation operator (for F = 1 atoms) A\ = 

— 2$j_ 1 $L ] J /y/3. Generally, one finds the many- 
body spin state with atom number N = F tot + 2A2, spin 



Ftot 5 an d longitudinal spin Mp is written as ( Ho and Yip 



2000 Koashi and Ueda 2000) 



\N,F toU M F ) = Z-^ 2 A\"* F: 



N 2 



?F to 



-Mi 



^ tot I vac) 



(38) 

where here F_ is the collective spin lowering operator, Z 
is a normalization constant, and |vac) is the zero particle 

state. The state with |F tot = 0, M F = 0) ex A\ N 2 |vac) 
is thus seen as one in which all particles form spin-singlet 
pairs. From the form of A\, one sees that all single- 
atom magnetic sublevels in the many-body spin state 
are equally populated, a result also understood from the 
isotropy of the spin-singlet state. This distribution is dis- 
tinctly different from that expected for a mean-field polar 
state. The sublevel populations of the many-body ground 
state also exhibit large fluctuations, with the standard 
deviation of the order of N ( |Ho and Yip| [2000). 

Acc ording to the definition of Pen rose and On sager 
(1956), a system is Bose-Einstein condensed if one of the 
eigenvalues of the reduced single-particle density matrix 
is extensive, i.e., of the order of the number of particles. 



Nozieres and Saint James (1982) pointed out the possibil- 



ity that there is more than one such macroscopic eigen- 
value and called that state a fragmented Bose-Einstein 
condensate. The F = 1 spinor condensate with > 
at zero magnetic field is thus fragmented. Indeed, this 
state can be characterized as a superposition of mean- 
field Bose-Einstein condensate states (Mueller et al. 



2006), 



|F tot = 0,M F 



<m\e,t 



(39) 



where | #,(/>) is the polar coherent spin state aligned with 
the axis specified by the polar and azimuthal angles 
and (j). 

More generally, it has been shown that the spin state 
of an F = 1 condensate within the SMA can be expressed 
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in the form (Barnett et al., 2011a 2010) 



/ 



dn^(e,4>)\e,<p), 



(40) 



where \p(6,(p) is a properly normalized probability am- 
plitude function on the unit sphere. The state of the 
condensate is thus mapped onto that of a quantum ro- 
tor. When the quadratic Zeeman term q is absent, -0(0, <p) 
distributes uniformly over the entire unit sphere and the 
state (|40j reduces to Eq. p9|. When q > and the 
number of atoms is large, ip(&, (p) tends to be localized, 
and the mean-field polar solution is favored. It has been 
shown that the complete spectrum of the many-body 
Hamiltonian of Eq. (|3tJ) with c^p > agrees with the 
lowest set of eigenvalues of a quantum rotor Hamilto- 
nian. A similar representation has been developed for 
spin-2 condensates ( |Barnett et al , 2011a[ ). 

It has been argued that fragmented Bose-Einstein con- 
densates are fragile against symmetry-breaking perturba- 
tions ( |Mueller et a/. , 2006). In the present case, it can 
be shown that when the longitudinal spin of the system 
is M^, the population of the \m z = 0) Zeeman state is 



calcu lated to be (Ho and Yip 
2000) 



2000 Koashi and Ueda 



N -M F 
2M F + 3 



(41) 



This result suggests that if the isotropy of the system 
is broken and the system acquires magnetization Mp, 
the m z = population decreases rapidly; it decreases 
to y/N as M F increases to and the distribution 

of populations rapidly approaches a mean-field result 
(N/2, 0, N/2). Thus the mean-field theory breaks down 
when the system is exactly isotropic, but it is rapidly 
restored in the presence of symmetry-breaking perturba- 
tions. 



gases to accommodate a very rich variety of topological 
excitations such as vortices, monopoles, skyrmions, and 
knots. 



A. Order parameter manifolds 

W ithin the SMA phase diagrams developed above (Sec. 
V.A| ), each condensate phase is represented by a stan- 
dard spinor ip$. For example, the ferromagnetic and po- 
lar phases of a spin-1 BEC are represented by (1,0, 0) T 
and (0, 1, 0) T , respectively. In the absence of an external 
magnetic field and neglecting the MDDI, the free energy 
of the system is invariant under the £7(1) gauge trans- 
formations and the SO (3) rotations in spin space. The 
entire group of symmetry transformations that leave the 
free energy invariant is given by G = £7(1) x 50(3). 
Therefore, a general order parameter is given by per- 
forming a gauge transformation by the phase <p and 
rotating ip Q to an arbitrary direction specified by Eu- 
ler angles a, /?, and 7: ip = e^[7(a, /?, 7)^0, where 
[7(a,/3, 7) = e ~' lF z^ e -tF y f3 e -iF z ^^ we w -vj uge foe gu | > 

script S to denote such geometric (spin) rotations. The 
dependencies of 0, a : and 7 on space and time charac- 
terize spin textures. 

Each phase of the condensate has distinct symmetry 
characterized by the isotropy group whose elements 
h leave ip unaltered: hip = ip. Then, the order parame- 
ter space is given by the coset M = G/H. Acting every 
element m of M on the standard spinor generates a com- 
plete set of spinors, called an orbit, which represent the 
order parameter manifold ( jMermin , 1979). 

Let us consider the order parameter manifolds for sev- 
eral spinor condensate phases. The order parameter of 
the F = 1 ferromagnetic condensate transforms under 
arbitrary gauge transformations and spin rotations as 



VI. SPIN TEXTURES OF DEGENERATE SPINOR GASES 

A single coherent (i.e., unfragmented) Bose-Einstein 
condensate can be described by the order parameter 
which is given by the eigenfunction corresponding to the 
largest and extensive eigenvalue of the reduced single- 
particle density matrix. Here by extensive, we mean that 
the eigenvalue, which gives the number of condensate 
particles, is proportional to the total number of parti- 
cles. Apart from the particle-number density distribu- 
tion, a scalar Bose gas has only the phase degree of free- 
dom whose spatial variation gives rise to mass current. 
A degenerate spinor Bose gas has the additional degrees 
of freedom of spin whose spatial variation produces not 
only mass currents but also spin currents. 

Spatial variations of the spin configuration are called 
spin textures. Spin textures can vary over space and time 
in response to an external magnetic field and boundary 
conditions as well as due to spin-dependent contact inter- 
actions and the MDDI. Their rich structure allows spinor 
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(42) 

The crucial observation here is that the gauge angle (p 
and the rotation angle 7 about the spin direction appear 
as a linear combination <p— 7, so that the order parameter 
is invariant under the simultaneous gauge transformation 
and the spin rotation about its orientation by the same 
amount. Therefore, H = [7(1)^+5 , where the subscripts 
(p + S 1 is displayed to emphasize the combined nature of 



the group operations, 
given by (|Ho) 11998) 



The order parameter manifold is 



M 



17(1) x SO(3) s 



17(1) 



SO(3) 4 



(43) 



The polar phase of the F = 1 spinor condensate trans- 
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forms as follows: 
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V2 



sin f3 



V2 



cos /3 
sin P 



(44) 



This expression illustrates that the polar state has a U(l) 
symmetry with respect to rotation about the quantiza- 
tion axis because ( p does not depend on 7. Also, we 
see that rotating the \m z = 0) polar state by tt about 
the y axis generates a minus sign, as can be seen by 
making the change j3 — » j3 + tt in Eq. ( 44 ) . Combining 



this rotation with a gauge transformation by 7r returns 
the original state. Thus the polar state manifold has a 
discrete Z2 symmetry. Altogether, the isotropy group is 
given as H = U(l)s^ x (7i2)<j>,s B and the order parameter 
manifold is given by ( | Zhou , [2001 [ ) 



m p = U(l)* x 50(3) s 



x5§ 



C/(l)s 7 x (Z 2 ) , S/3 (Z 2 )^ 



(45) 



We note here that the spin and gauge remain coupled in a 
discrete manner. As discussed in Sec. |VI.D[ this discrete 
spin-gauge symmetry causes the polar-state manifold to 



support half-quantum vortices (Leonhardt and Volovik 



2000). 



Spin-2 condensates have even richer order parame- 
ter manifolds. The ferromagnetic phase with ipQ = 
(1,0,0,0,0) T has the same symmetry as that of the spin-1 
case except that the order parameter manifold (Eq. ([43]) ) 
is divided by ^2)0, s ( |Makela et al\ |2003[ ) because the 
phase winding is doubled compared with the spin-1 case. 
The uniaxial nematic phase with = (0,0, 1,0, 0) T 
looks similar in form to the polar phase of a spin-1 con- 
densate, but has one important difference. Both of these 
systems have the spin inversion symmetry about an ar- 
bitrary axis perpendicular to the quantization axis, but 
only the spin-1 order parameter changes sign upon in- 
version. Thus, only the spin-1 condensate can couple 
the spin inversion with the gauge and generate a half- 
quantum vortex accompanied by mass circulation. The 
biaxial nematic phase with ipo = (1,0,0,0,1) T has four- 
fold symmetry about the quantization axis and tt rotation 
symmetry about four perpendicular axes. Therefore the 
isotropy group is the dihedral-four group ( Da)^, arid 
the o rder parameter manifold is given by ( |Song et~al 
2007) 



M 



>iaxial _ U(lh X SO(3)< 



The cyclic/tetrahedral phase has the tetrahedral symme- 
try with the isotropy grou p H = T. The order parameter 
manifold is thus given by ( Makela et al , 2003 Semenoff 
and Zhoul 12007) 



M cyi 



clic _ UiXU x £Q(3)s 
(T) ,s 



(47) 



In both these cases, as in the polar case, geometric rota- 
tions and inversions are accompanied by gauge transfor- 
mations to return to the original quantum state. Because 
two elements in dihedral-four and tetrahedral groups do 
not commute in general, some vortices created in these 
phases are non-Abelian (Kobayashi et ah, 2009). A list 



of symmetry groups of spin-1 and spin-2 condensates as 
we ll as liquid crystals and s uperfluid 3 He systems is given 
by |Kobayashi et aZ.| ( [2012] ). 



B. Superflow and spin-gauge symmetry 

In a spatially extended spinor Bose-Einstein conden- 
sate, the spin order parameter can everywhere minimize 
the local spin-dependent interactions that arise from s- 
wave scattering, but still be characterized by spatial vari- 
ation in the order parameter within its degenerate man- 
ifold. In the absence of magnetic dipole-dipole interac- 
tions, the resulting spatial pattern, a spin texture, is a 
low-energy phenomenon, with excess energy only due to 
the kinetic energy. As exhibited in Eq. (42 ) for the F = 1 



ferromagnetic case, the combined spin-gauge symmetry 
or the order-parameter manifold implies that local spin 
rotations (spatially varying Euler angles) can be regarded 
equivalently as local variations in the condensate phase 
0, to which the co ndensate responds by supercurrent flow 
( |Ho and Shenoy 1996). Conversely, a supercurrent can 
decay by forming a spin texture. 

For the F = 1 ferromagnetic spinor gas, this connec- 
tion is exhibited in the expression for the superfluid ve- 
locity: 



vf = ^[V(<A- 7 )-cos/?Va], 



(48) 



where M is the atomic mass. Unlike for the scalar Bose- 
Einstein condensate, in which the superflow is irrota- 
tional (V x v s = 0), her e, vf satisfies the Mermin-Ho 



relation (Mermin and Ho 1976): 



V x vf = 4 sin/3V/3 x Va. 



(49) 



This result implies that the superflow in the ferromag- 
netic phase is not a potential flow, and consequently the 
circulation alone is not quantized. Rather, the difference 
between the circulation along an arbitrary contour C and 
the Berry phase enclosed by it is quantized: 



(46) j> v ^ r_ ^/ (l-cos/3)Vadr = x integer. (50) 



For the F = 1 polar phase, the superfluid velocity is 
calculated to give 



vf = — V0 

s M 



(51) 



which is irrotational, and hence the circulation alone is 
quantized. It should be noted that although Eq. (51) 



27 



assumes the same form as that of scalar condensates, the 
quantum of circulation is h/2M which is half of the latter. 
This is because If the phase <j) changes only tt around a 
loop, the single- valuedness condition can be met if the 
Euler angle f3 changes by tt simultaneously, as can be 
seen from Eq. (44). 



C. Spin current and spin nematicity 

A spin current can flow in the absence of mass current 
by carrying the same number of particles with different 
spin states in the opposite directions. In the absence of an 
external magnetic field and the MDDI, the magnetization 
satisfies the continuity equation: 



dt 



V-jf n = 0, 



(52) 



where 



2Mi 



E (^)mn [0*„ - (VOn] (53) 



is the spin current. However, additional terms that 
break spin conservation appear on the right-hand side 
of Eq. ( 52 ) in the presence of an external magnetic field 



or the MDDI flKudo and Kawaguchi[ |2010| ): 



^ + v . jr = ^ (bx ^ + E (zx ^ + | e „ 



tflZV 1 v ZV 1 

(54) 

where Cdd is the strength of the MDDI, b is the dipole 
field, p and q are the strengths of the linear and quadratic 
Zeeman effects, respectively, and e^ zv is the Levi-Civita 
tensor, and N^ v is the spin nematic tensor. This equation 
shows that the dipole field and the linear Zeeman term 
generate spin torques and cause spin precessions, while 
the quadratic Zeeman term combined with spin nematic- 
ity produces a spin-nonconserving term. This last term 
implies that even if the initial state of the system is non- 
magnetic, a system having spin nematicity can dynami- 
cally develop transverse magnetization. Once the trans- 
verse magnetization grows, the MDDI leads to longitudi- 
nal magnetization due to the first term on the right-hand 
side of Eq. ([54]). 
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FIG. 8 (Color) (a) Three-dimensional configuration of the 
n field (shown as arrows) of a knot with the unit linking 
number, where n = (n x ,n y ,n z ) T gives the normalized spinor 
order parameter as e l6 ((—n x + in y )/y/2, n z , (n x + in y )/y/2). 
In the figure, only n on the x, y, and z axes are shown. 
The dashed circle follows the direction n = (0,0, —1) T . The 
color represents the value of n z according to the scale shown, 
(b) The torus on the figure shows the isodensity surface of 
|V>-i| 2 = 0.47 with the color showing the phase of ac- 
cording to the color gauge on the left. The white and black 
tubes on the right are the cutouts on which n points in the 



directions (0, 0, — 1) 



< -0.95) and (1,0,0)' 



> 0.95), 



respectively, The two tubes cross once (linking number 1). (c) 
Cross sections of the density of the m = ±1 (left) and m = 
(right) components in the x-y plane. The torus shapes of the 
m = dzl components manifest themselves as double rings in 
the cross section, which can serve as the signature of a knot. 
Figure adapted from Kawaguchi et al. ( 2008 ) . 



D. Topological excitations 

Topological excitations are the defects that are sta- 
ble against weak perturbations. Such stability is en- 
sured by discrete quantum numbers that characterize 
the topology of the order parameter manifold. Natu- 
rally, the types of topological excitations depend cru- 
cially on the properties of the order parameter manifold 
which, in turn, is determined by the phase of the con- 
densate. A rich variety of order-parameter manifolds of 
spinor condensates allow many distinct types of topolog- 
ical excitations. Examples include integer and fractional 



vortices (Isoshima et al., 2001 Leonhardt and Volovik 



20001 |Makela et al] |2003| |Semenoff and Zhouj |2007[ 
Yip[ [19991 |Zhou| |2001[ ) , non- Abel ian vortices ( [Kobayashi 
et al] |2009| >7V Ho oft-Pol y akov flStoof et al] |200ip and 
Dirac monopoles (Blaha[ |1976 Ruostekoski and An- 



glin, 2003 fVolovikl 2003), skyrmions (particle-like soli- 
tons) QA1 Khawaja and Stoof] |2001 



Volovik and Mineyev[ | 1977| ), and knots (Kawaguchi et al. 
20081). 



Shankar 1977 



The underlying physics that supports fractional vor- 
tices is discrete spin-gauge symmetries. Consider, for 
example, the polar state of a spin-1 condensate whose 
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order parameter changes its sign upon inversion in spin 
space. Now, we consider a closed loop in real space 
and assume that the order parameter undergoes inver- 
sion in spin space as one circumnavigates the loop. The 
single-valuedness of the order parameter is met if the 
order parameter has an additional phase change of 7r, 
which is realized if the winding number of the state is 
a half integer. Therefore, under such circumstances, the 



system possesses half-quantum or Alice vortices (Leon- 



hardt and Volovik 2000). Let us next consider the 



cyclic / tetrahedral phase of a spin-2 condensate. As dis- 
cussed earlier, the order parameter has a three- fold sym- 
metry with respect to 2tt/3 and 4tt/3 geometric rotations 
about the symmetry axes of the tetrahedron. This geo- 
metric rotation adds a phase of 2tt/3 or 4tt/3 to the spinor 
wavefunction. Consider a loop in real space and assume 
that the order parameter rotates geometrically by 2tt/3 
as one makes a complete circuit of the loop. The single- 
valuedness of the order parameter is only met if the gauge 
phase changes by an integer multiple of 2tt/3. Thus, the 
cyclic phase can accommodate vort ices with circulations 
of one-third or two-thirds of h/M ( |Semenoff and Zhou 
2007). 



The biaxial nematic and cyclic/tetrahedral phases of a 
spin-2 condensate can host non-Abelian vortices because 
their isotropy groups are non-Abelian ( jKobayashi et al] 
2009). The unique feature of non-Abelian vortices lies in 
the fact that the generators of vortices do not commute 
with each other. Consequently, when two vortices col- 
lide, they cannot reconnect or pass through each other. 
Rather, they form a rung vortex that sticks the two vor- 
tices together. This implies that a collection of vortices 
will eventually form an interwoven network in sharp con- 
trast to the case of Abelian vortices which usually recon- 
nect upon collisions and tend to produce vortices with 
smaller scales. Such a distinction is expected to lead to 
a major difference in statistical properties of quantum 
turbulence. 

A spin-1 polar condensate is predicted to accommo- 
date knot excitations (Kawaguchi et al. , |2008| ). Here the 
objects that form knots are field lines along which local 
spins point in a given direction. Suppose that a polar con- 
densate is held in an optical trap and then a quadrupole 
magnetic field is applied on the condensate. Then, the 
magnetic moment of each atom undergoes Larmor pre- 
cession according to the local magnetic field. Numeri- 
cal simulations show that one after another knot, which 
consists of two loops linking together with unit linking 
number, enters the condensate from its periphery. Ex- 
perimentally, such a knot formation can be probed by 
performing a Stern- Gerlach experiment on a sliced con- 
densate; if there is a knot, one can observe a double-ring 
density distribution on every spin component. 

Experimental efforts to create and study these myriad 
topological objects are ongoing. Investigations thus far 
have revealed two such objects: the polar-core spin vor - 
tex ( |Sadler et al\ |2006| ), s kyrmions (|Choi et 



et a/.[|2009a" ). 

The polar-core vortex was observed in F = 1 87 Rb con- 
densates following a quench of an \m z = 0) condensate to 
a setting of the quadratic Zeeman shift q that favored fer- 
romagnetic ordering; the relevant phase diagram is shown 
in Fig. [5] and the dynamics of such a quantum quench are 
detailed in Sec. [IXJ Summarizing, the quench leads to a 
dynamical instability that introduces spontaneous spin 
currents within the initially uniform condensate. The 
spontaneous symmetry breaking occurs independently in 
different regions of the condensate. Altogether, this spon- 
taneously generated spin-current field is then expected 
to contain circulating spin currents. The polar-core spin 
vortex represents one quantum of spin-current circula- 
tion. An isolated polar-core spin vortex may be described 
by the following condensate wavefunction (in cylindrical 
coordinates (r, 0)): 



( f(r) x e ie 
g(r) x 1 
y h(r) x e~ l{ 



(55) 



The radial functions /(r), g(r) and h(r) have the lim- 
iting values lim^oo (/, g, h) = (1/2, 1/V% 1/2) and 
lim r ^o(/ 7 9-, h) = (0,1,0). The former limit yields a 
transversely magnetized texture away from the vortex 
core, a configuration that is locally in the ground state for 
the ferromagnetic F = 1 condensate with q c± 0. The lat- 
ter limit applies because the \m z — ±1) components have 
non-zero circulation, and thus they must vanish at the 
origin. However, the non-rotating \m z = 0) component 
remains within the vortex core, where the spinor con- 
densate remains in the higher-energy polar state. That 
is, such a topological defect captures the "false vacuum" 
represented by the pre-quench state of the condensate. 
A polar-core spin vortex of opposite circulation would be 
described by replacing 9 —> — in Eq. (55). 

The polar-core vortex was identified using in situ 
magnetization-sensitive imaging (Fig. [9]). A closed path 
within the condensate texture was identified as having 
non-zero transverse magnetization along the entire path, 
with the magnetization orientation winding by ±27r. 
Within that closed path, the magnetization was found 
to be consistent with zero, allowing this magnetization 
winding to be identified as a polar-core spin vortex, which 
has a non-magnetic core, rather than as a Mermin-Ho 
texture, which would be longitudinally magnetized at its 
core. The polar-core spin vortex represents the sponta- 
neous breaking of chiral symmetry and is reproduc ed in 
numerical si mulations of such quench experiments ( |Saito| 



et al, 2006). 



and Mermin-Ho textures (Jae-yoon et al> 2012 Leslie 



2012), 



Two-dimensional skyrmion spin textures were gener- 
ated in an antiferromagnetic F = 1 spinor condensate in 
the polar state. Such a skyrmion is defined as a defect 
in which the entire order parameter manifold is mapped 
continuously onto the plane, in such a manner that the 
order parameter is identical at the periphery of the two- 
dimensional system - or at infinity in an infinitely ex- 
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FIG. 9 (Color online) In- situ detection of a polar-core spin 
vortex. Spatial maps of the transverse magnetization (a) mag- 
nitude (IMjcj/I) and (b) orientation (0) are shown for a gas 
imaged 156 ms after being quenched from the polar to the 
ferromagnetic state. Data from a portion of the image are 
magnified, showing the transverse magnetization (c) magni- 
tude and (d) orientation and also (e) the magnitude of the 
longitudinal magnetization (|M 2 |). The transverse magneti- 
zation orientation along a closed path (indicated in yellow) 
shows a net winding of 27T, revealing the presence of a spin- 
vortex defect. The core (within the closed path) shows no 
significant longitudinal magnetization, identifying the defect 
as a polar-core spin vortex. G ray scale fo r images (a,c,e) is 
the same. Figure adapted from Sadler et al. ( |2006| ). 



tended sample. Such a skyrmion was created by prepar- 
ing the condensate in a spatially uniform polar state with 
the nematic director n aligned with the z axis, and then 
rotating the condensate spin in a coherent, spatially de- 
pendent manner to generate a texture with n = +z at 
the center and n = — z at the outer edge of the gas. To 
effect such a rotation, following earlier realizations of a 
similar technique with ferromagne t ic textures ( Leanhardt 
et al. H2002||Nakahara ef^|2000| ), |Choi et 



2012F 



pro- 
duced a spherical-quadrupole magnetic field, of the form 
B = B' [xx + yy — 2zz] where the Cartesian coordinates 
are defined with respect to the location of zero field. The 
location of the field zero was gradually translated per- 
pendicular to and across an optically trapped 23 Na con- 
densate (Fig. 10). Everywhere within the condensate, 
the magnetic field reverses its orientation about an axis 
that varies with the azimuthal angle in the plane of 
the condensate. Far from the trajectory of the field zero, 
the field is always strong and the atomic spin follows the 
rotating field adiabatically, rotating the director n by 7r 
radians. At the location of the field zero, the magnetic 
field changes its direction diabatically; the atoms cannot 
follow this field reversal, and thus their nematic director 
remains unchanged. At intermediate distances from the 
field zero, the spin texture smoothly interpolates between 
these settings, generating a continuous skyrmion texture. 

The texture was characterized by time-of-flight mea- 
surements of the Zeeman- state distributions, preceded by 
rf-pulses which allowed for a characterization of the phase 
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FIG. 10 (Color online) Creation of a skyrmion in a two- 
dimensional antiferromagnetic F = 1 spinor condensate, (a) 
A spherical quadrupole magnetic field, with field orientations 
shown by arrows, is translated from being centered above to 
below the Bose condensate (yellow horizontal region at the 
center). The nematic director n follows the rotating field at 
large distances from where the magnetic field zero crosses the 
condensate, and remains fixed near the crossing, generating 
the skyrmion texture, (b) Stern-Gerlach absorption imaging 
shows the \m z — 0) component is divided between an inner 
region (where n = +z) and an outer ring (where n = — z). 
Between these regions, the nematic director lies in the trans- 
verse plane; the state is thus a superposition of the \m z = dzl) 
sub levels as seen in the data. The skyrmion decays after 100 
ms. At long times, the entire gas returns to the \m z — 0) 
state. Figure reproduced from Choi et al. ( 2012| . 



winding in each of the spin components. The skyrmion 
was stable for about 100 ms. The mechanism for its de- 
cay has yet to be identified. 

Applying a similar procedure to an F = 1 condensate 
prepared initially in a uniform ferromagnetic state pro- 



duces a Mermin-Ho spin texture ( 


Jae-yoon et al. 2012 


Leanhardt et al 


2002| Nakahara et al t 2000). This tex- 



ture differs from a skyrmion because the edge of the con- 
densate, while magnetized in a common direction, car- 
ries a 47r phase winding. This difference between the 
polar and ferromagnetic states is reflective of the fact 
that the SO (3) group does not support two-dimensional 
skyrmions. 

A texture similar to the Mermin-Ho ferromagnetic tex- 
ture was produced in an untrapped F = 2 87 Rb spinor 



gas using optical Raman transitions (Leslie et al. , 2009a). 
Starting with a gas spin-polarized in the \m z = 2) state, 
circular polarized light was used to drive Ami? = 2 tran- 
sitions, populating the \m z = 0) and \m z = —2) states. 
One of the beams used for this transition was in the first- 
order Laguerre-Gaussian mode, so that the Raman tran- 
sition also imprinted a +2tt phase winding onto the prod- 
uct Zeeman state. The spin configuration produced with 
such Raman transitions resembled the Mermin-Ho tex- 
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tures of ferromagnetic F = 1 spinor condensates, which 
are composed of m z = (+1,0,-1) states with circula- 
tions of (0, 1, 2) quanta. 

Topological excitations are also expected to appear 
upon rotating spinor condensates. Ferromagnetic spin- 
1 condensates are expected to accommodate rotation by 
forming Mermin-Ho textures , which are shown to be 
stabilized by such rotation (Mizushima et al. 2002b 



Mueller, 2004). Under certain circumstances, such spinor 



condensates may also adopt non-axisymmetric vortices, 
ones in which all Zeeman components are circulating, but 
in which the vortex cores of the components do not over- 
lap. Such spin configurations are n o longer fully ferro- 
magnetic (Mizushima et al., 2002a). Antiferromagnetic 



spin-1 condensates are expected to store angular momen- 



tum in the form of half- vortices ( |Isoshima and Machida 



2002 Mueller, 2004). The landscape of spin condensates 



under rapid rotation is expected to be very complex (Kita 



et 



2002 Mizushima et al. 2004). 



E. Hydrodynamic description 

The Gross-Pitaevskii equation deals with the order pa- 
rameter, from which physical quantities such as magne- 
tization, spin nematicity, mass current and spin current 
can be calculated. Alternately, one may develop descrip- 
tions that are equivalent to the Gross-Pitaevskii equation 
and deal directly with the quantities of experimental in- 
terest. Such hydrodynamic theory has been instrumental 
in describing scalar condensates, and it is of interest to 
ask if it is possible to describe spinor condensates in terms 
of observable quantities only. 

To produce manageable equations of motion, one may 
confine the description to low-energy variations within 
the order-parameter manifold. Lamacraft derived such 
a set of low-energy hydrodynamic equations for the case 
of a spin-1 ferromagnetic condensate (Lamacraft 2008). 



When the particle density is constant and the system is 
locally fully magnetized, the normalized magnetization 
vector F obeys a modified Landau-Lifshitz equation that 
describes an advection of the magnetization by a super- 
fluid velocity v: 



(ft+v-V)F = Jfx V 2 F, 



(56) 



where the superfluid velocity satisfies the incompressibil- 
ity condition V • v = 0. A self-contained set of hydrody- 
namic equations that are equivalent to the spin-1 Gross- 
Pitaevskii equation involve the spin quadrupole tensor in 
addition to t he particle density, spin density and mass 
current (Yukawa and Ueda, 2012) 



The spin-gauge coupling of the spinor wavefunction 
permits some arbitrariness as to how the overall phase 
of the wavefunction is apportioned between variations of 
the spin texture orientation and of the superfluid phase. 
The choice of such apportionment defines a vector poten- 
tial, a, which enters into the definition of the superfluid 



velocity: v = V(j) — a, where <fi is the superfluid phase. As 
in the vector potential in electromagnetism, this vector 
potential has an arbitrariness of the gradient of a scalar 
potential, but its curl is fixed by the Mermin-Ho relation: 



V x a = - — tijkFiVFj x VFfc, 



(57) 



In analogy of electromagnetism, the right-hand side may 
be interpreted as an effective magnetic field which is gen- 
erated by spin textures. 

In addition to the vector potential, we may also in- 
troduce a scalar potential r] via the Berry connection 
r] = i(ip\dt\ip). The set (a, 77) then forms a four- vector 
a a , and by analogy with electromagnetism, we may in- 
troduce the electromagnetic tensor f a /3 = d a a$ — dpa^ 
whose time and space components give the electric and 
magnetic fields, respectively, and the right-hand side of 
the Landau-Lifshitz equation ([56]) now includes an ana- 



logue of the Lorentz force ( jBarnett et al. 2009 Volovik 
19871. 



Including the effects of the MDDI changes the hydro- 
dynamic descripti on because it permits n on-local spin- 
spin interactions. iTakahashi et al. (2007) included the 



MDDI in a hydrodynamic description of classical spins 
(ignoring spin-gauge coupling). A hydrodynamic treat- 
ment for the ferromagnetic condensate including the 
MDDI and the quadra tic Zeeman energy is discussed by 
Kudo and Kawaguchi (2010). 



VII. EXAMINING THE CONNECTION BETWEEN 
MAGNETIC ORDER AND BOSE CONDENSATION 

Magnetism in solid-state systems may be ascribed 
broadly to two different origins: the interaction among 
local magnetic moments, or the interaction among itin- 
erant fermions. In spinor Bose gases, we encounter mag- 
netic ordering that occurs in a system of itinerant bosons. 
We argue that this magnetization occurs due to Bose- 
Einstein condensation, by which the influence of weak 
interactions is amplified through bosonic enhancement. 

Here, we explore the question of how this statistically 
enhanced magnetic order responds when the effects of 
bosonic enhancement are weakened. We consider three 
scenarios: increasing the temperature toward or above 
the Bose-Einstein condensation temperature, restricting 
the dimensionality of the system to one or two dimen- 
sions so as to enhance the role of quantum and thermal 
fluctuations, and placing the spinor Bose gas within a pe- 
riodic potential so as to increase the low-energy density 
of states and enhance the role of interactions. 



A. Spinor Bose gases at non-zero temperature 

The magnetic properties of a gas of non-interacting 
bosons can be derived from basic statistical mechanics 
( |Yamada] |1982[ ) . A phase diagram evaluated as a func- 
tion of the applied magnetic field B and the temperature 
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FIG. 11 (Color) Phase diagram for a uniform non-interacting 
F = 1 spinor Bose gas as a function of applied magnetic 
field B and temperature T. The normal density for the 
m z Zeeman sublevel is determined by its chemical poten- 
tial /io — \iBvn z /ksT ', where /x is the magnetic moment and 
/jlo is determined by number conservation. The unmagne- 
tized gas at B = condenses at the critical temperature 
T c (0) = 2<Kh 2 /mk B x (nC(3/2) /3) 2/3 determined by the den- 
sity n/3 of each Zeeman sublevel. For \jiB /ksTc^ | ^> 1, the 
critical temperature (solid black line) increases by a factor 
(2F + l) 2 / 3 (indicated by the dashed line). Contours show 
the longitudinal magnetization, with color bar shown. The 
superfluid state is ferromagnetic, as indicated by the discon- 



(o) 



tinuous magnetization at B — and T < T c 



T, given in Fig. [TT] specifically for a spin-1 gas, shows two 
notable features. First, the quantum degeneracy temper- 
ature varies with the magnetic field. This occurs because 
the Zeeman energy differentiates the chemical potential 
of the different Zeeman sublevels, increasing the frac- 
tional population, and hence the phase-space density, in 
the lowest-energy Zeeman state. Thus, the Bose-Einstein 
condensation temperature T c increases by a maximum 
factor of (2F + l) 2 / 3 in the limit that the Zeeman split- 
ting greatly exceeds the thermal energy. Second, ferro- 
magnetism accompanies Bose-Einstein condensation, as 
evident in the diverging zero-field susceptibility of the gas 
below the Bose condensation temperature. It may be said 
that ferromagnetism is a "parasitic" phenomenon, taking 
advantage of the zero entropy of the condensed fraction 
to establish magnetic order. Refined treatments of this 
system, e.g. utilizing renormalization group theory and 
evaluating critical exponents, support the picture that 
the two phenomena of condensation and magnetization 



occur through the same transition (Caramico D'Auria 
[19961 [FYota et al\ [19841 [Guj |2QQ3| |Simkin and 



et 



Cohen 1999). 



In the case considered here, the ferromagnetic phase 
- as opposed to nematic, cyclic, or other phases - is 
selected due to the application of a symmetry-breaking 
magnetic field and the absence of spin-dependent inter- 
actions. This ferromagnetic tendency persists even under 
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FIG. 12 (Color) Thermodynamic phase diagram of a 
chromium spinor Bose gas. The large dipolar relaxation 
rate of chromium allows the longitudinal magnetization (ab- 
scissa, defined as the per-atom spin divided by K) to achieve 
its thermodynamic equilibrium value as a function of the 
applied magnetic field and the gas temperature (ordinate, 
scaled by the transition temperature T c o of a fully magne- 
tized gas). The Bose-Einstein condensate transition temper- 
ature increases with magnetization, as predicted (see Fig.fTTj). 
Below the condensation temperature, the Bose-Einstein con- 
densate is either fully magnetized (region B) or partially mag- 
netized (region C), according to whether the Zeeman energy is 
larger or smaller than the non-ferromagnetic spin-dependent 
interaction energy. Figure from Pasquiou et al. ( 2012| . 



arbitrary spin-independen t interactions, as can be shown 
both for pseudo-spi n 1/2 (Sig gia and Ruckenstein 1 1980 



and Liebl 2002). 



Yang and Li, 2003]) and higher-spin systems (Eisenberg 



This basic thermodynamic picture was confirmed in ex- 
periments with chromium gases (|Pasquiou et al 2012). 



Due to the large dipolar relaxation rate in this gas, the 
longitudinal magnetization is a free parameter, deter- 
mined by thermodynamics, and also by spin-dependent 
interactions. The Bose-Einstein condensation tempera- 
ture was observed to vary with magnetic field, and hence 
with the magnetization of the £ as, as predicted (Fig.[l2|). 
Below the Bose-Einstein condensation temperature, the 
gas adopts a magnetization that is determined by both 
the linear Zeeman effect and also the spin-dependent in- 
teractions in the gas (Pasquiou et al, 2011[): at higher 



magnetic field, where the Zeeman energy dominates, the 
condensate is fully magnetized, while at lower magnetic 
field, where the non-ferromagnetic spin-dependent inter- 
action dominates, the condensate occupies multiple Zee- 
man states. 

In the absence of dipolar relaxation, such as for the 
alkali spinor Bose gases, the equilibrium states of spinor 
gases need be evaluated under the constraint of constant 
longitudinal magnetization, rather than as occurring at 
constant magnetic field. To the extent that the mag- 
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netization in Fig. 11 spans all values for T > T c (0) , the 
constrained equilibrium state can be taken as that which 
occurs at the value of B (now a Lagrange multiplier) that 
gives the correct longitudinal magnetization. 

However, for T < , the equilibrium diagram of Fig. 
11 fails to describe spinor Bose gases within a range of 
longitudinal magnetization around zero. In this case, a 
gas of non-zero longitudinal magnetization is predicted to 
undergo a sequence of Bose-Einstein condensation transi- 



tions ( 


Isoshima et al. t 2000 Kis-Szabo et al. 


2007 


Zhang 


et al. 


2004 


). A first condensation, at T > T 


produces 



a purely longitudinally polarized condensate. Below this 
temperature, the non-condensed populations in all Zee- 
man states approach a common value, whereupon, at 
T < Tc°\ Bose-Einstein condensation occurs into mul- 
tiple Zeeman sublevels. The magnetic order of that low- 
temperature condensate depends on the character of the 
spin-dependent interactions. 

Experimental evidence of this two-step condensation 
phenomenon has been obtained by the Paris group work- 
ing with sodium condensates of non-zero magnetization 5 . 
Two-step condensation can also occur due the quadratic 
Zeeman shift, which, similarly, can imbalance the distri- 
bution of atoms within the non-condensed gas so as to 
initiate Bose-Einstein condensation into a single Zeeman 
sublevel (the \m z = 0) state if q > 0) before condensation 



2011) 



can occur into other states as well (Phuc et al 

The role of spin-dependent interactions on the Bose- 
Einstein condensation transition has been debated. In 
a mean-field description, magnetic order in the non- 
condensed or the condensed gas produces a Weiss field 
through the spin-dependent s-wave interaction. Calcu- 
lating thermodynamic functions in the presence of this 
molecular field indicates that, in the case of a ferromag- 
netic F = 1 spinor gas (c^ < 0), a ferromagnetic state 
has a lower free energy than the paramagnetic state at 
the Bose-Einstein condensation temperature, i.e. ferro- 
magnetism occurs separately and before Bose-Einstein 
condensa tion, even for a weak ferromagnetic s-wave inter- 
action flGu et al.||2004||Gu and Klemm||2003HKis-Szabo 



et 



2005 Tao et 



2008) 



However, 



Natu and Mueller 



(2011) point out that such 



a calculation neglects the dominant effect of exchange 
interactions. As discussed in Sees. III. A and III.D[ zero- 
range interactions take place between identical bosons 
or fermions if their total spin F pa i r is even. For exam- 
ple, for two electrons in an atom, the strong Coulomb 
repulsion at short range is thereby minimized for spin- 
triplet states. For two bosons, the repulsion due to a 
strong spin-independent s-wave interaction is minimized 
by pairing the bosons in an odd spin state. This direct 
bosonic exchange effect disfavors ferromagnetism in the 
nondegenerate gas, and thus a weak spin-dependent in- 



teraction cannot yield ferromagnetism before condensa- 
tion. For condensed atoms, the direct exchange is absent, 
and weak spin-dependent interactions can indeed dictate 
the magnetic order of the gas. 

The magnetic properties of spinor Bose gases close 
to the Bose-Einstein condensation transition were ex- 



amined experimentally by Vengalattore et al. (2010). 
A non-magnetized gas was prepared well above the 
Bose-Einstein condensation transition, and then grad- 
ually cooled by evaporation from an optical trap. A 
time-of-flight analysis indicated the formation of a low- 
momentum, quantum degenerate gas at the temperature 
expected given the atom number per Zeeman state 
and the presence of a harmonic confining potential. In- 
situ imaging revealed local magnetic ordering at all tem- 
peratures below Tc°\ However, the experimental sensi- 
tivity to magnetization and the level of control over the 
gas temperature were insufficient to settle the question 
of whether magnetization appeared exactly at, or else 



(o) 



slightly above, T c 



B. Two- and one-dimensional systems 

Another interesting arena in which to examine the 
connection between magnetic order and superfluidity is 
in low-dimensional systems. For example, in a scalar 
Bose fluid in two dimensions, the link between superflu- 
idity and condensation is severed. Below the Berezinskii- 
Kosterlitz-Thouless (BKT) phase transition, the fluid be- 
comes suddenly superfluid, via a jump in the superfluid 
stiffness, but the gas eschews off-diagonal long-range or- 
der and is thus not Bose-Einstein condensed. The transi- 
tion can be understood from the behavior of topological 
defects in the scalar fluid - vortices of single positive or 
negative quantum of circulation - that are bound below 
the BKT transition and unbound above it. Signatures of 
the BKT transition have been observed in atomic gases, 
where the presence of a trap transforms the transition 
into a gradual crossover, and where the low-temperature 
phase is indeed Bose-Einstein condensed in that the co- 
here nce length becomes as large as t he finite size of the 
gas (Hadzibabic and Dalibard, 2011). 



Topological structures in spinor Bose-Einstein conden- 
sates differ from those of scalar condensates, and thus 
one expects different behavior in two dimensions. Iden- 
tifying the correct symmetry group of the F = 1 polar 
spinor Bose-Einstein condensate, Mukerjee et al. (2006) 



F. Gerbier, personal communication, 2012 



discussed the role of half-quantum vortices in the behav- 
ior of such a gas in two dimensions. From analytical ar- 
guments and supported by numerical studies, they iden- 
tified a BKT-like transition from a normal to a superfluid 
state. Given the lower kinetic energy of the half-quantum 
vortex, they predict such a transition is characterized by 
a larger jump in the superfluid stiffness than is observed 
in the scalar case, a tell-tale sign to be sought in future 
experiments. It is important to note that this larger jump 
was predicted much earlier in considerations of classical 
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two-dimensional nematic films flde Gennes , 1971) and of 
p-wave superfluids (|Kor shunovl 119851). 



Remarkably, this topological phase transition leads to 
a gas that is superfluid, but still lacking in long-range 
spin order. This phenomenon can be understood by con- 
sidering the Hamiltonian of the gas in the form (James 



and Lamacraft 2011) 



H 



t n (Vn) 2 +t (V0f -qn\ 



(58) 



Here, the gas is restricted to the ground-state manifold, 
and thus defined by the locally varying director n(r) of 
the polar state and the superfluid phase 0(r). The quan- 
tities t n and t<f) denote the magnetic and phase stiffness 
(kinetic energy) of the gas, derived for example from a 
lattice model of the fluid. I n the q = case consid- 
ered by |Mukerjee et al (2006), the dynamics of n and <\> 
are decoupled, owing to the lack of spin-gauge coupling 
for the polar state, except for the topological linking of 
nematic disclinations and vorticity of the two fields, re- 
spectively, by the half-quantum vortices. Below the BKT 
transition, the half-quantum vortices of opposite charge 
become bound to each other, and the two fields become 
unlinked from one another. Algebraic order develops in 
the quantity e 2 *^, yielding superfluidity without conden- 
sation. However, the director n, which explores the larger 
space of rotations, cannot order according to the Mermin- 
Wagner-Hohenberg theorem, and thus no nematic spin 
order is obtained. 

The application of a quadratic Zeeman energy breaks 
the rotational symmetry of the polar-state manifold. For 
q > 0, one expects the spins to order at sufficiently low 
temperature at an Ising-type phase transition (Fig. 13). 



A polar superfluid is expected also at q < 0, though such 
a scenario has been considered only und er the additional 
influe nce of tuning t n and separately ( |Podolsky et al 
2009). Numerical simulations have been performed tak- 



ing into account the finite size of a trapped polar gas, 
revealing similar spatial variations between coherent and 
incoherent portions of the gas as observed for trapped 
scalar gases, though the precision of such calculations is 
difficult to assess without a reliable non-zero temperature 



theory for the spinor Bose gas (Pietila et al 2010). 



In one dimension, the scalar Bose gas is described by 
Luttinger liquid theory, and, for delta-function interac- 
tions, can also be treated exactly by the coordinate and 
thermodynamic Bet he Ansatze. At high density, the low 
temperature phase is a quasi-condensate, essentially de- 
scribed by mean-field theory but lacking true long-range 
order. At low density, in the Tonks- Girardeau regime, the 
gas properties resemble those of non-interacting fermions 
as the repulsive interactions suppress the short-range 
two-body correlation function (Cazalilla et al , |2011 ). 

A one-dimensional Bose gas with spin still admits 
an exact solution in the case of fully symmetric delta- 
function interactions. For example, for the F = 1 spinor 
gas with SU(3) symmetric s-wave interactions (corre- 
sponding to = 0), the ground state is a "color fer- 




FIG. 13 (Color online) Phase diagram for a two-dimensional 
F = 1 polar spinor Bose gas vs. temperature and quadratic 
Zeeman shift. The dashed red line indicates the BKT transi- 
tion from the normal to the spin-disordered superfluid state 
(PS), marked by t he binding of half-quantum vortices ( Muk 



erjee et al 2006). A non-zero quadratic Zeeman shift fa- 
vors the longitudinal polar state. A spin texture containing 
a bound half- quantum- vortex pair necessitates regions where 
the polar- state director deviates from the z axis; these regions 
acquire a domain- wall width and tension defined by q. At high 
temperature, in the PS state, these domain walls proliferate, 
causing the spin texture of the superfluid to be disordered. 
Below an Ising-type phase transition (blue line), the domain 
walls are confined, allowing for long-range spin ordering and a 
conventional 2D superfluid (S). Here, t is the hopping energy 
in a discrete lattice model of the 2D gas. Figure reproduced 
from James and Lamacraft (2011). 



romagnetic" liquid: the phonon excitations and phase 
coherence properties are those of the scalar Luttinger liq- 
uid, while the spin sector shows loca l magnetic orde r and 
free-particle-like (oc k 2 ) excitations ( |Li et al\ |2003[ ). 

Significantly different behavior is seen for antiferro- 
magnetic spin-dependent interactions. As discussed for 
the low-temperature 2D polar gas, the dynamics of 
the gas separate into the charge-phase dynamics and 
the nematicity-director dynamics. The zero-temperature 
state shows a power-law decay in the spin-singlet ampli- 
tude Ai, a state which can be regarded as a spin liquid 



that is gapped against spin-triplet excitations (Cao et al 



2007 Essler et al 2009) 



The separated dynamics of the charge and spin de- 
grees of freedom hark to the phenomenon of spin-charge 
separation expected for one-dimensional Fermi systems. 
Here, the mass and spin degrees of freedom of the spinor 
Bose gas evolve sep arately, show i ng, fo r example, differ- 
ing sound velocities (Essler et al, 2009). However, in the 
Bose system, this separation is not unique to one dimen- 
sion, appearing also in three-dimensional gases. 

Adding a non-zero magnetization to the gas transforms 
the p olar spin-liquid s tate into a ferromagnetic liquid 
state. |Lee et al\ ( |2009| ) suggest that, between these two 
limits, the two liquids coexist as a two-component Lut- 
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tinger liquid. In this regime, the spin excitations of the 
liquid are no longer gapped. 



C. Spinor Bose gases in optical lattices 

An optical lattice makes it possible to control the filling 
fraction - the number of atoms per lattice site - realizing 
an artificial carrier-doped many-body quantum system 
that features, among other things, the superfluid-Mott- 



insula tor (SF-MI) transition ( |Greiner et al. , 2002 Jaksch 
et al. |1998[ ). The fluctuations in the number of atoms 
per lattice site are controlled by changing the ratio of the 
tunneling amplitude t to the on-site Hartree interaction 
(charge gap) Uo. At zero temperature, if the ratio t/Uo 
is below a critical value, the system is in the MI phase; 
otherwise, the system is in the superfluid phase. With a 
spinor gas of spin-F atoms loaded into the lattice, new 
energy scales concerning spin-dependent interactions, U{ 
(i = 1, • " come into play, making the system an 
ideal playground to study quantum magnetism. The 
dynamics of spinor gases in optical lattices also offers 
unique possibilities because spin-mixing dynamics can 
be coupled with density mod ulations such as Bloch os- 
cillations and self-trapping (Li et al., 2010). An optical 



lattice loaded with spinor gases can also be utilized to 
determine various scattering lengths of spinor gases pre- 



cisely ( Widera et al.\ 2006) and for quantum information 
processing ( |Widera et al.\ 2005 ) . 

A spinor gas in a spin-independent optical lattice 
can be described with the Bose-Hubbard Hamilto- 



nian ( Jaksch et al. , 1998| ) supplemented by the spinor 
interaction at each lattice site. For the spin-1 case, the 
Hamiltonian is given by 



(ij),m 



(59) 



where ^ m is the annihilation operator of a boson in 
the magnetic sublevel mp — 1,0,-1 at site z, hi = 
5^ m / 0j m ^im? an d /i is the chemical potential that con- 
trols the filling fraction of the system. The spin op- 
erator f\ has three components F ia = '$l rn (F a ) rnn 'ipi n 
(a = x, y, z), where F a is the a component of the spin-1 
operator. The ratio Ui/Uq is about 0.04 for 23 Na and - 
0.005 for 87 Rb (see Table[nj). The control parameters that 
determine fundamental properties of the system are t/Ui, 
filling fraction, dimensionality, and size of the atomic 
spin. 

When the number of atoms per lattice site is large, a 
ferromagnetic mini-condensate on each lattice has a large 
magnetic moment and interacts with other condensates 
through the MDDI. Such an array of large spins is shown 
to yield ferromagnetism (|Pu et al.\ |2001[ ) and magnetic 
solitons (Li et al. , 2005 Xie et all 2004) in one dimension, 



and d omain structures in higher dimensions ( Gross et al. 
[2002] ). 

When the number of atoms per lattice site is small, 
the MDDI can be ignored. Then, ferromagnetic conden- 
sates are expected to show similar behavior as scalar con- 
densates with respect to the MI-SF transition. For the 
antiferromagnetic case, the state in each Mott lobe is ex- 
pected to show properties such as a nematic insulator and 



a spin-singlet insulator (Demler and Zhou, 2002 Imam 



bekov et al. 2003), where the nematic state corresponds 



to the polar state of a bulk condensate. The nematic 
state breaks spin-rotation symmetry but preserves time- 
reversal symmetry and has gapless spin- wave excitations, 
while the spin-singlet state has spin-rotation symmetry 
and gapped excitations (spin gap). A remarkable feature 
of the antiferromagnetic MI is an even-odd parity effect 
in which even fillings stabilize the MI state because all 
particles can then form spin-singlet pairs with the spin 
gap in addition to the cha rge gap ( |Demler and Zhou 



2002| |Tsuchiya et a/.|[2004| ). 



In one dimension, antiferromagnetic spin-1 systems 
with odd fillings and small t/Uo are predicted to be 
dimerized in such a manner that adjacent atoms tend to 
form spin-singlet dimers with an excitation gap (|Demler| 



and Zhou[ [2002] [Yipl [2003] |Zhou and Snoek[ |2003| ). The 



underlying physics is that virtual hopping to neighboring 
sites gives rise to an effective antiferromagnetic coupling 
between adjacent sites. Since each atom has two nearest 
neighbors to be dimerized, the dimer phase has two-fold 
degeneracy and breaks lattice translational symmetry. 
A renormalization-group method confirms that fo r odd 
fillings the insulating phase is always dimerized ( |Rizzi| 
et al.\ [2005 ), although the obtained dimer amplitude is 



very small and the state might be quantum nematic near 
the SO(3)-symmetric point (Lauchli et al. 2006). On 



the other hand, quantum Monte Carlo calculations show 
that the entire first Mott lobe has the dimerized ground 
state up to a relatively large t/Uo which may be probed 
by the Bragg spectroscopy with peaks locating at half- 



integer multiples of the inverse lattice vector (Apaja and 



Syljuasen, 2006). A Monte Carlo study also shows that 



the SF-MI transition is first order for even fillings due to 
the spin gap of spin-singlet pairs, whereas it is s econd or- 
der for odd fillings as in the ferromagnetic case ( Ba trouni| 
[^rZ|[2009l ). 

The spin-1 system with small t/U{ can be de- 
scribed with an effective Hamiltonian called the bilinear- 
biquadratic Hamiltonian ( De Chiara et a/.^2011| |Imam-| 
bekov al] |2003| |Rodriguez et a/.| |2011| |Yip[|2003| ) 



(ij) 



(60) 



where the cos term arises from the virtual second-order 
hopping to neighboring sites, and D characterizes the 
strength of the quadratic Zeeman effect. In one dimen- 
sion, this model describes many different phases such as 
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the Haldane phase including the Affleck-Kennedy-Lieb- 
Tasaki point (De Chiara et aZ. , 2011). 



In higher dimensions, insulating states with an odd 
number of atoms per site are predicted to be nematic, 
whereas those with an even number are either nematic or 



singlet with a first-order transition between them (Imam 



|bekov et al. [ [2003). Based on the mean- field or Gutzwiller 



approximation, it has been shown (Kimura et al. 2005 



Krutitsky et al. 2005 ) that the SF-MI transition may be 



first order around the tip of each Mott lobe. This is in 
sharp contrast with the ferromagnetic case in which the 
phase bound ary is alwa ys second order, as in the scalar 
case (Fisher et al. , 1989). The physical origin of the first- 



order transition is that in the MI phase the spin state 
is fixed, whereas in the superfluid phase it is not; thus 
the spin symmetry can have a discontinuity at the phase 
boundary. A possible experimental signature of the first- 
order transition is a discontinuous jump in the transverse 
magnetization from zer o in the MI phase to a nonzero 
value in the SF phase (Kimura et al.\ 2005). The spin 



gap in the MI phase for even fillings leads to quantized 
magnetization plateaus, while the second-order virtual 
hopping of spinor atoms induces an effective quadrupolar 



interaction between spins (Imambekov et al. 2003 Snoek 
et al. 2009| ) (see the second term in Eq. ( 60 f) which leads 



to a canted nematic order in which magnetization is not 
quantized but changes contin uously as a functi on of an 
external magnetic field (Imambekov et al., 2004). As the 



ratio Ui/Uo increases, even Mott lobes expand and odd 
Mott lobes shrink until the latter disappear at a critical 
value around Ui/Uq = 0.5 (Lacki et al.\ 2011|), presum- 



ably because the spin gap due to spin-singlet pair for- 
mation works in favor of the even Mott lobes. At unit 
filling, phase diagrams in the presence of the quadratic 
Zeeman effect were investigated based on the density- 
matrix renormalization method and Lanczos diagonaliza- 
tion ( [Rodriguez et a/.[|2011 ), showing the field-induced 
transition between the transverse and longitudinal mag- 
netization (or nematic) for the ferromagnetic (or antifer- 
romagnetic) case. 

For the case of spin-2 atoms, we may expect even richer 
phase diagrams flZhou and Semenoff [ [2006] ) . The phase 
diagram was st udied by treating the intersite hopping 
perturbatively QEckert et al] |2007a[ |Hou and Ge[ |2003) 



and by a variational approach ( Snoek et al. t 2009 ) . 



n one 



dimen sion, a trimerized state was predicted ( |Chen et al] 
[2012] ). The SF-MI transition in spin-3 condensates with 



dipolar interactions was also discussed (Bernier et al. 



2007) 



We note that these theoretical predictions for quantum 
magnetism in optical lattices are extremely appealing, 
but also very challenging to achieve. The spin-dependent 
energies between atoms trapped at low filling within op- 
tical lattices are typically much smaller than the on-site 
interaction energies that have produced the Mott insu- 
lating states studied to date. The task of lowering the 
temperature of lattice-trapped gases to the point where 
such energy scales dominate, without the aid of the Bose 



enhancement which makes spinor Bose-Einstein conden- 
sates accessible, is a major hurdle. Moreover, open- 
ing the spin degree of freedom for lattice-trapped gases 
may add a large amount of entropy per particle that 
may be difficult to reduce thereafter by evaporation or 
other techniques. In this regard, it is worthwhile to men- 
tion two major experimental advances: the demagnetiza- 
tion gradient cooling method demonstrated by the MIT 
group ( |Medley et al] 2011) and the Pomerantchuk-type 
cooling observed with spinful fermions of Yb by the Ky- 
oto group (Sugawa et al. 2011). 



VIII. DIPOLAR SPINOR GASES 

The magnetic dipole-dipole interaction (MDDI) cou- 
ples the spin degrees of freedom with the orbital degrees 
of freedom. For magnetic dipoles \i x and i± 2 separated 
by the relative coordinate r, the interaction energy is 

K, d (r) = c dd ^-^- 3( ^ 3 r/r)( ^- r/r) . (61) 

The MDDI is invariant under simultaneous rotation of 
spins and coordinates, but not under the separate rota- 
tion of each. Therefore, in the presence of the MDDI, 
the constraint of global spin conservation, assumed fre- 
quently above, no longer strictly holds. The MDDI con- 
serves the total angular momentum, and, thus, a dipole- 
induced relaxation of the spin angular momentum should 
lead to a change in the orbital angular momentum, caus- 
i ng a previou sly stationary spinor gas to begin rotating 
(|Kawaguchi et al. 2006a). This phenomenon is analo- 



gous to the Einstein-de Haas effect in a ferromagnet ([Ein- 
stein and de Haas, 1916), where the spin angular momen- 



tum is converted to a mechanical rotation. Conversely, 
a rotation of the system can lead to its magnetization 



- a phenomenon known as the Barnett effect (Barnett 
19151). 



In many experiments on cold atomic gases, the MDDI 
is negligible. To assess its importance, one may com- 
pare the MDDI to the spin-independent contact interac- 
tion, which determines the compressibility and excitation 
spectrum of a spin-polarized Bose-Einstein condensate. 
To make such a comparison, recall that the magnetic 
field inside a sphere with uniform magnetization M is 
B = (2/3)/ioM, where /iq is the vacuum permeability. 
Setting M = fin for a gas of atoms with magnetic mo- 
ment fi and density n, the magnetic-dipole energy per 
atom, = (l/3)/ioM 2n 7 shows the same scaling with 
density as the mean-field interaction energy per atom 
E = Anh 2 an/M with a being the scattering length and 
M the atomic mass. The ratio of their energies is thus 
independent of density: 



Cdd,a 



127rh 2 a 



(62) 



For alkali atoms, this ratio is quite small away from 
Feshbach resonances where the interaction energy can be 
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tuned near zero ( jFattori et al. , 2008 Pollack et al. 2009); 
for example, for 87 Rb in the F = 1 state, edd,a = 0.002. 
Thus, the long-range dipole interaction can typically be 
neglected in experiments on scalar alkali gases. In con- 
trast, for Cr (edd,a = 0.15), and even more so for Dy 
(^dd,a — 1 for a typical scattering length), the dipolar 
interaction is an important factor in the structure and 
stabil i ty of a spin-polarized gas (see reviews by Baranov 
et a/.|(|2002l);|Baranov|(|2008|);|Lahaye et aZ.| ([2QQ9|)). 



However, for spinor Bose gases, the MDDI can have 
significant effects, both because it is spin-dependent and 
because it is long-ranged. As a spin-dependent interac- 
tion, the MDDI may compete with the spin-dependent 
contact interaction to determine the local spin ordering 
of the degenerate spinor gas. Here, the importance of 
the MDDI may be assessed by examining the ratio of 
the spin-dependent contact to the dipolar interaction en- 
ergy, given as edd,Aa i n which we now use the scatter- 
ing length difference Aa that governs the spin-dependent 
interactions. For example, in the F = 1 manifold of 
87 Rb, considering Aa = (<22 — ao)/3 = 0.25 A, we find 
e dd,Aa = 0.4, so that we may expect dipolar interactions 
to pla y a significant role, as firs t pointe d out by |Yi et 
(120041) and as discussed in Sec. IVIII.C 



Even if the MDDI does not disrupt the local spin or- 
der selected by the spin-dependent contact interaction, 
it may nevertheless dictate the long-range spin ordering 
of the gas. Here, the proper comparison is between the 
MDDI and the quantum confinement energy, h 2 /(2ml 2 ). 
Equivale ntly, one ma y define a dipolar healing length, 
^dd — \/h 2 /(2M^dd) and compare it to the length I of 
the condensed gas. As a result, for spatially extended 
samples, the ground state spin configurations are pre- 
dicted to be inhomogeneous spin textures that spon- 



taneously rotate (Kawaguchi et al. 2006b Yi and Pu 



2006), as discussed in Sec. VIII.D 



A. Dipolar interactions in the cold-collision regime 

The s-wave approximation cannot be used to treat 
the MDDI, owing to its peculiar scattering properties in 
three dimensions. The scattering phase shift 5t(k) for 
the lih. partial wave is proportional to the wave vector k 
in the low-energy limit k — »• 0, regardless of £. This im- 
plies that all partial waves contribute to the low-energy 
scattering (Landau and Lifshitz 1977). Moreover, the 



anisotropy of the MDDI induces coupling between differ- 
ent partial waves. 

Rather, the MDDI is retained explicitly in the inter- 
action Hamiltonian. At low magnetic fields, where the 
hyperfine spin F is a good quantum number, the matrix 
elements of the magnetic dipole operator are given by 
{F,m\n\F 1 m , ) = gFl^B^mm' ? where gp is the Lande g- 
factor. The MDDI is then written in second quantization 



as 



Kid = ( ~y II d,r 1 dr 2 : 



F-F-3(F-r/r)(F-r/r) 



(63) 

where in the integral r = ri— r 2 , and Cdd = Ho(9f^b) 2 
The coupling between different hyperfine manifolds can 
be ignored because the hyperfine splitting ~ 100 mK is 
much larger than the ~ 1 nK energy scale of the MDDI. 



B. Interactions between rapidly precessing dipoles 

The rapid precession of the gas magnetization un- 
der a strong applied magnetic field alters the effective 
form of the MDDI. For conditions where the strength 
of the dipolar field, on the order of B ~ gpfiofiBn, 
is much smaller than the applied field, we account for 
the rapid precession in the interaction between two spa- 
tially fixed dipoles by identifying the dipole moments 



^( rot ) i n the Larmor precessing frame as fi£ 0t ^ = Hz an d 
^(rot) ^ ^ rot ) ± ifiy 0t ^ = e ±lUJLt fi±, where we have spec- 
ified z as the magnetic field axis b. Taking the time 
average, we obtain the effective interaction as 



(rot) 

Mo Ml 



Q(r 



(rot) 



47T 



where 



1 - 3((r/r) • b) 2 3b^b u - 5^ 



(64) 



(65) 



The effective MDDI, written now for arbitrary orien- 
tation of the magnetic field, distinguishes the longitu- 
dinal from the transverse magnetization, and expresses 
a non-local XXZ-type interaction between atomic spins. 
This interaction can effect variations in the magnetiza- 
tion even when the dipolar field is small compared to the 
applied field, for the reason that the dipolar field pro- 
duced by a transversely magnetized region of the gas is 
modulated at precisely the Larmor precession frequency. 
This time- varying field can cause resonant Rabi nutation 
of the magnetization in another portion of the gas. 



C. Dipolar interactions in the single spatial-mode regime 

Let us revisit the mean-field description of spinor Bose- 
Einstein condensates in the single-mode approximation, 
presented above in Sec. |V.A| wherein the condensate spin 
character is described simply by a spatially independent 
spinor wavefunction. 

Accounting for the MDDI introduces a new spin- 
dependent energy, potentially changing the local order 
preferred by the gas. Within the SMA, this additional 
energy is determined by taking the spatial integral of Vdd 
- or, in the presence of rapid Larmor precession, of VJ<f - 
over the condensate density distribution. Considering a 
density distribution that is cylindrically symmetric about 
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the z axis, the spatial integral of V dd yields the fo llowing 
addition to the energy functional ( |Yi et al 2004): 



^dd = c dd n (2(F Z ) 2 - 
dri dr 2 



c 'dd 



Cdd 

n 



(F X ) 2 ~ (Fyf 

n(ri)n(r 2 ) 2 



ri -r 2 



(1-3 cos 2 e 12 ) 



where n(r) is the condensate density, and #12 is the polar 
angle of the vector ri — y 2 . 

We see that the dipolar interaction adds an energy sim- 
ilar to that of the spin-dependent contact interactions, so 
that local effects of the MDDI can be regarded, roughly, 
as changing the value of c\. However, the MDDI energy 
term is anisotropic, i.e. the coupling between spin and or- 
bital degrees of freedom causes the spatial anisotropy in 
the density of the condensate, determined by the trap- 
ping geometry, to translate into spin-space anisotropy. 
By varying the trapping geometry, one can change the 
sign of c' dd and thus tune the effective sign of the spin- 
dependent dipolar energy. 

|Yi et q/.| p004) examined the impact of dipolar interac- 
tions on the exact many-body groun d-state pr edicte d for 
the F = 1 case by |Law et oZ] ( |1998D (see Sec. |V.B[ ). For 
ferromagnetic interactions, the dipolar interaction breaks 
spin-rotational symmetry, so that the ground state be- 
comes the state of maximum spin, \F tot = N, F Zj tot — 
Mp), where now Mp = N for c' dd < (prolate) and 
Mp = for c' dd > (oblate). The dipolar interaction 
also stabilizes these maximal-spin states as the ground 
state even in the presence of weak antiferromagnetic in- 
teractions. A similar stabilization of the ferromagnetic 



ground states was noted by Kjall et al (2009), who con- 



sidered also the impact of the quadratic Zeeman shift on 
the F = 1 mean-field phase diagram. Dipolar interac- 
tions, treated within the single-mode approximation, are 
also predicted to affect the dynamics of the c ondensate 
magnetization und er applied magnetic fields ( jYasunaga 
and Tsubota 2011 ), an effect which will ultimately affect 



the accuracy of spino r-gas based magnetometers (Ven- 
galattore et al, 2007) (discussed in Sec. XI). 

The dipolar interaction introduces additional spin- 
mixing terms to the Hamiltonian, e.g., terms of the form 
(^q^q^+i^-i + ^^l^o^o) f° r the case of the F = 1 
spinor gas. Thus, the MDDI modifies the spin-mixing 
dynamics, such as coherent oscillations and dynamical 
instabilities, discussed in Sec. IX For example, dipolar 



interactions can induce a transition across the separatrix 
that divides two types of spin-mixing trajectories, so that 
the spin-mixing dynamics from a given initial spin state 
will depend strongly on the trapping geometry (Rong 



et al, 2005). Dipolar interactions are also predicted to 



break the rotational symmetry of dynamically unstable 
modes in single- mode and spatially extended \m z — 0) 
condensates quenched to low values of the quadratic Zee- 
man shift. This dipole- induced anisotropy makes the 
critical quadratic Zeeman shift qo governing the instabil- 
ity depend on the orientation of the z axis with respect 
to the trap shape, and should also produce anisotropic 



patterns of spin domains after the quench (Sau et al 
[2009] ). 

D. Flux-closure relation 

The MDDI favors spin waves and magnetic domain 
structures. This can be seen by expressing it in momen- 
tum space: 

(Vdd) oc J d 3 k [3|F(k) • k/k\ 2 - |F(k)| 2 ] , (66) 

where F(k) is the Fourier transform of the magnetization 
vector F(r). The dipole interaction energy is minimized 
by minimizing the first term and maximizing the second 
term on the right-hand side. The latter is achieved by 
placing a large weight on a particular spin-wave compo- 
nent F(k), while the former can be achieved by requiring 
that k be perpendicular to F(k). In real space, this con- 
dition is written as V-F(r) = 0, which implies that no net 
magnetic flux emanates from any point in real space, that 
is, the magnetic flux closes upon itself. This is known as 



the flux-closure relation (Landau et al 



1984]), which fa- 
■state ferromag- 



vors magnetic domain structures in solic - 
nets to minimize the magnetic-field energy. 

In addition to the flux-closure relation, a ferromagnetic 
Bose-Einstein condensate features spin-gauge symmetry 
(see Sec. VLB ) that induces supercurrent when the spin 
texture of the system is nonuniform. As discussed in 
Sec. VI.E[ the spin textures also give rise to vorticity 
via the Mermin-Ho relation (Eq. (57)). Since trapped 



systems in general produce nonuniform spin textures via 
the MDDI, we may expect a supercurrent to flow even in 
the ground state. In fact, the F = 1 ferromagnetic Bose- 
Einstein condensates may undergo symmetry breaking 
into at least three different phases that feature inhomo- 
geneous spin textures and circulating currents (Fig. [l4| . 
These phases are selected by changing the ato m number 
and trap freq uencies ( |Huhtamaki et al 2010| [Kawaguchi 



et al, 2006b). 



The phase diagram of dipolar spinor condensates at 
zero magnetic field has also been examined as a func- 
tion of the dipole-dipole interaction and trap aspect ratio 
for both ferrom agnetic and antiferromagnetic cases by Yi 
and Pu (2006), and spin textures are predicted to emerge 



in a pancake trap and relatively strong dipole-dipole in- 
teractions. Spin t extures in the MD DI-dominated regime 
are discussed by Takahashi et al (|2007[) under the as- 



sumption that the spin behaves as a classical vector. 



E. Experimental signature 

Effects of the MDDI may have been observed in the 
dynamics of spin textures studied by the Berkeley group 
( |Vengalattore et al\^00^ . A ferromagnetic F = 1 87 Rb 
condensate in a quasi-two-dimensional trap was prepared 
in a uniform transversely magnetized state. Thereupon, 
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FIG. 14 (Color) (a) Phase diagram of a ferromagnetic spinor- 
dipolar condensate, where PCV, CSV, and FL stand for polar- 
core vortex, chiral spin-vortex, and flower phases, respec- 
tively, and £ sp , £dd, and Rtf denote the spin healing length, 
the dipole healing length, and the Thomas-Fermi radius, re- 
spectively. The solid curve and dashed line show the first- 
order and second-order phase boundaries, respectively. The 
total angular momentum J per particle is for PCV, 
CSV and FL. The red and black circles with arrows 
mass and spin circulations, respectively. The F = 1 
spinor gas traces the dotted line indicated as "Rb line, 
(d) Typical spin texture in each phase. The top and bottom 
panels show the unit spin vector projected onto the x-z and 
x-y planes, respectively. The color of the arrows represents 
the magnitude of the normalized sp in density according to th e 
bottom scale. Figure adapted from ||Kawaguchi et a/.| (|2006b). 
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a gradient of the magnitude of the magnetic field was 
briefly applied. Larmor precession in the presence of this 
gradient gave rise to a helical spin texture, in which the 
orientation angle <fi of the magnetization in the transverse 
spin plane increased linearly with distance along the di- 
rection of the field gradient. The pitch of the helix was 
determined by the strength and duration of the pulsed 
gradient. Thereafter, the texture was allowed to evolve 
under a uniform applied magnetic field before the gas was 
imaged. 



Considering just the s-wave interactions, a spin helix 
with a long pitch A is a low-energy texture, with ex- 
cess energy per particle given just by the kinetic energy 
h 2 /2m\ 2 . It was therefore surprising that such a helix 
underwent a rapid evolution toward a finely modulated 
spin texture (Fig. 15) in which, given the short typical 
length scale of spin domains, the kinetic energy is much 
larger than its initial value. 

The MDDI provides an additional source of energy to 
the spinor Bose-Einstein condensate, and may account 
for the observed dynamics. To test this possibility, a 
rapid sequence of rf pulses with random spacing was ap- 
plied to randomize the orientation of the atomic spins 
with respect to spatial coordinates. Such rf pulses had 
the effect of diminishing the strength of the short length- 
scale modulation, supporting the picture that the MDDI 
causes such modulation to emerge. We note that rf-pulse 
techniques have been suggested also as a way to isolate 
and emphasize, rather than to suppress, the effe cts of the 
MDDI in spinor gases (Yasunaga and Tsubota, 2008). 



The stability of the helical texture was explored 
theoretically, both from the perspective of a low- 
energy hydrodynamic theory of the ferromagnetic state 
( Lamacraftj 2008 ), a nd als o using a Bogoliubov treat- 
ment ( |Cherng et al\ |2008[ ). The spectrum of modula- 
tional instability was found to reflect the anisotropy of 
dipolar interactions, suggesting a mode- softening toward 



a spatially pe riodic ground-state texture (Cherng and 
2009[ ). The possible emergence of modulated 



Demler 



ground-state spin textures due to dipolar interactions was 
explored, focusing on the balan ce between the quadratic 
Zeeman effect and the MDDI (|Kawaguchi et al. , 2010 
Kjall et all 120091 iZhang and Hoi 120101). When the 



quadratic Zeeman energy is positive and stronger than 
the MDDI, the longitudinal magnetization is suppressed, 
and the transverse magnetization forms a helix. In 
the opposite case in which the MDDI dominates the 
quadratic Zeeman energy, the longitudinal magnetiza- 
tion forms a domain structure to minimize the MDDI. 
However, these works have not yet provided a consistent 
explanation for the experimental findings ( |Kawaguchi| 



et 



2010) 



IX. NON-EQUILIBRIUM QUANTUM DYNAMICS 

Cold atomic gases, and specifically spinor Bose gases, 
allow one to examine non-equilibrium dynamics in a man- 
ner that would be difficult with solids or liquids. In con- 
trast with condensed-matter systems, which are studied 
typically near equilibrium, cold-atom materials can be 
prepared far from equilibrium by applying state prepa- 
ration techniques or rapidly changing the system Hamil- 
tonian. Atomic gases equilibrate slowly. The kinetic de- 
grees of freedom relax on millisecond timescales defined 
by the rate of elastic collisions. The internal degrees of 
freedom in spinor gases equilibrate even more slowly. The 
decay of the total longitudinal magnetization is either 
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FIG. 15 (Color) Dissolution of helical spin textures in a ferromagnetic F 



7 Rb spinor Bose-Einstein condensate. A transient 



magnetic field gradient is used to prepare transversely magnetized (b) uniform or (a, c) helical magnetization textures. The 
transverse magnetization column density after a variable time T of free evolution, in a homogeneous magnetic field, is shown in 
the imaged plane, with orientation indicated by hue and amplitude by brightness (color wheel shown), (b) A uniform texture 
remains homogeneous for long evolution times, while (c) a helical t exture with pitch A = 60 /xm dissolves over ~200 ms, evolving 
into a sharply spatially modulated texture. Figure from Vengalattore et al. (2008). 



slow (as for 52 Cr) or nearly imperceptible (as for many 
alkali gases). The weak spin-dependent contact and dipo- 
lar interactions also produce dynamics over time scales - 
many milliseconds - that are long compared to the time 
needed to vary the trapping potential, the applied fields, 
or to probe the system. 

In addition, unlike condensed-matter materials, cold- 
atom systems generally lack true reservoirs with which 
to exchange energy, angular momentum, or magnetiza- 
tion. This raises the question of whether cold-atom sys- 
tems equilibrate at all, whether equilibration is achieved 
within a given time frame for certain properties and 
not for others, and, if an equilibrium state is achieved, 
whether that state differs from a straightforward thermal 
equilibrium state. 

Here, we highlight studies on spin mixing dynamics in 
degenerate spinor Bose gases. Such dynamics are asso- 
ciated with redistribution of atomic populations among 
spin states, the formation of spin domains and spin tex- 
tures, symmetry breaking and phase transitions, and 
nonlinear matter- wave optics. 



A. Coherent spin mixing in single mode experiments 

1. Microscopic spin dynamics 

The spin-dependent contact interactions in a spinor 
gas can cause temporal oscillations of spin-state popu- 
lations. To exhibit this effect, consider a system with 
two F = 1 atoms in the ground state of a tightly con- 
fining trap for which we make the single-mode approxi- 
mation. Consider an initial state ^(0) = |0,2,0) written 
as a Fock state in the eigenbasis of a projection of the 
single-particle angular momentum, with magnetic quan- 



tum numbers mp = {1,0,-1}, respectively. Owing to 
its being a superposition of F tot = and F tot = 2 states, 
the initial state evolves as 



-iojot 



V2(l 



-iuj\t\ 



|0, 2, 0) - v 3 qi,0,l> 

(67) 

where ujq = (cq 1 ^ + c^)n/2H and uj\ = 3c^n/2h with n 
being the average density. In a spin-mixing collision, a 
pair of atoms in the \m z = 0) spin state scatter inelas- 
tically into the \m z = +1) and \m z = —1) states. This 
interaction is coherent, leading to periodic variation of 
the spin composition of the system. 



Such coherent microscopic spin mixing oscillations 



were observed by [Widera et al. (2005) (Fig. 16). Pairs 
of atoms were held in sites of an optical lattice, pre- 
pared with rf and microwave pulses, allowed to interact 
for a fixed time, and then probed by a Stern- Gerlach 
analysis. Each atom pair serves as a clean microscopic 
labora tory to measure the strength of the spin-mixing 
terms (IWidera et al. 2006), potentially also sensitive to 



the magnetic dipole interaction between individual cold 
atoms ( |Sun et al. , 2006| ). Spin mixing oscillations were 
observed within both the F = 1 and F = 2 manifolds 
of states of 87 Rb. The experiments explored also the 



impact of quadratic Zeeman shifts ( jGerbier et al. 2006) 



adding a controlled energy difference between the |0, 2, 0) 
and 1 1,0, 1) two-atom spin states. When this difference 
exceeds the spin-mixing interaction strength, the oscilla- 
tions increase in frequency and decrease in amplitude, as 
expected for off-resonant Rabi oscillations. 
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( 34 ) is then expressed as 



Hold time t (ms] 



Magnetic field 6(G) 



FIG. 16 (Color online) Spin-mixing oscillations of two 
trapped F — 1 87 Rb atoms, (a) An atom pair is prepared in 
the lowest motional state of an optical lattice site, in the ini- 
tial state 1 0, 2, 0) . Spin mixing collisions convert the atom pair 
to the state |1, 0, 1), which is offset in energy from the initial 
state by the quadratic Zeeman energy 2q oc B 2 . Spin mix- 
ing occurs in a large ensemble of such two-atom systems, (b) 
Population oscillations are measured by Stern-Gerlach/time- 
of-flight imaging, (c) Several coherent spin-mixing oscilla- 
tions are observed, (d) the frequency of which varies with 
the quadratic Zeeman shift. The frequency differs from that 
expected based on th e scattering lengths calculated by |van| 
Kempen et al (2002) (theory curved (dashed), with system- 
atic uncertainty indicated by shading); indeed, 20% discrep- 
ancies remain among reported value s of a<i — ap for the F = 1 
87 Rb system. Figures adapted from |Widera et al (2006). 



2. Single-mode mean-field dynamics 

Remarkably, the coherence of microscopic collisional 
dynamics of atom pairs is retained also in the spin-mixing 
dynamics of macroscopic collections of atoms. Such dy- 
namics in tightly confined Bose-Einstein condensates can 
be treated in a single-spatial-mode mean-field theory, i.e. 
specifying the condensate wavefunction, \I/ = $(r)^;(t), 
as a product of a common, fixed spatial wavefunction 
and a time- varying spin state. We consider the ensu- 
ing spin-mixing dynamics in the presence of the con- 
tact interaction and a quadratic Zeeman shift q along 
z. For an F = 1 spinor gas, while the spinor wavefunc- 
tion ip is specified by six real parameters, the dynamics 
is constrained by the overall U(l) gauge symmetry of 
the spinor wavefunction (associated with conservation of 
the norm of ip) and rotational symmetry about a given 
direction, say, z (associated with conservation of the z- 



projection of the spin, defined as m 
per atom), leaving just two dynamical variables (Kro- 



de Passos, 2004 Zhang et a/., 2005). Following Zhang 



njager et al\ [20051 [2008) |Pu et al\ Jl999j |Romano and 



et al, we define po = as the fractional population 

of the \m z =0) Zeeman state, and use the phase differ- 
ence — arg ( , + i'0_i'0o'0o) which controls the degree of 
magnetization or nematicity. The spin-dependent energy 
functional for the uniform spin-1 spinor gas given in Eq. 



£7(1) 



„(1) 



npo 



(1 - Po) + V(l - Po) 2 ~ m 2 cosO 



+9(1 - Po)- 



(68) 



This energy functional implicitly describes coherent 
spin dynamics in a single spin mode. The dynamical 
variables po and satisfy canonical relations, yielding 
equations of motion as = (2/K)dE /dpo and po = 
— (2/K)dE/d0. The energy E is conserved under such dy- 
namics; thus, the trajectory of po and follows a constant 
energy contour (Fig. 17). The energy extrema are thus 
identified as stationary states. Below, we will discuss 
extensively the stability of the higher-energy stationary 
states. 

Spin mixing occurs in two distinct dynamical regimes. 
Consider the case of zero longitudinal magnetization 
(m = 0). F or \q\ > 2|c ( j 1 ^|n, in the so-called "Zee- 
man regime" (|Kronjager et al 2006), all trajectories 



are running solutions in 0, and the spin undergoes oscil- 
latory orientation-to-alignment conversion - alternating 
between magnetized and nematic states - as expected 
for high-spin objects u nder quadratic spin Hamiltonians 



( |Fano and Macek| |1973| ). For \q\ < 2\c ( ± ) \n, the so-called 
"interaction regime," phase space is divided in two dy- 
namical regions: one region contains oscillations that are 
bounded in (akin to oscillatory pendular motion) , and 
the other contains running solutions in (similar to over- 
the-top pendular motion). Indeed, several works have 
drawn analogies between such spin-mixing dynamics and 



cillations (Chang et al 


, 20071 2005[ Jing et al\ 2009[ Ya- 


sunaga and Tsubota| 2 


010a|b[|Zhang et al\ 2005). How- 



ever, one should note the difference in topology between 
the phase space of a pendulum and that of the spin-1 
spinor gas under this representation. The phase space 
of a pendulum has the topology of an infinite cylinder: 
periodic in the angle and unbounded in the angular veloc- 
ity. In contrast, the po — phase space has the topology 
of a cone (for non-zero longitudinal magnetization) or a 
sphere (for zero longitudinal magnetization). 

The transition between the Zeeman and interaction 
regimes is signified dramatically in the frequencies of 
small-amplitude spin oscillations. Consider the case with 
m = 0, and a state initialized near the transverse mag- 
netized state (po = 1/2, = 0). Deep in the Zeeman 
regime, E^ is dominated by the term q(l — po), whence 
we obtain ~ —2q/h as expected for single-atom dynam- 
ics. In the interaction regime, for q = 0, the deviations 
5po — Po — 1/2 and from the transverse magnetized 
state remain small. To lowest order in these deviations, 
we find = — (2c^n/h) 2 0, so that spin mixing oscil- 
lations occur at the frequency 2\c^\n/h determined by 
the spin-dependent contact interaction. Between these 
regimes, the spin dynamics are slowed, reaching a "spin- 
mixing resonance" at \q\ = 2|c^ \n where the spin-mixing 
period diverges (Black et al 2007 Chang et al> 2005; 
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FIG. 17 (Color) Energy functional E (1) of an F = 1 spmor 
Bose-Einstein condensate with zero longitudinal magnetiza- 
tion (m — 0). Here, po = ^o^o and = arg (^i^-i^^o) for 
the spinor (^+1, ^o, ip-i) T written in the F z eigenbasis. For 
antiferromagnetic (ferromagnetic) interactions, minimum en- 
ergies are shaded in blue (red) . Spin mixing causes the spinor 
gas to evolve along constant-energy contours (gray lines). In 
the interaction regime (\q/c^n\ < 2), a non-inert energy ex- 
tremum emerges (marked by black dot) and also a separatrix 
(dashed curve), lying at — q (E^ = 0) for positive 

(negative) q/c^n, separates the bounded and running spin 
mixing trajectories. 



Kronjager et al. , 2006). The resonance value of q varies 
with magnetization m dBlack et all 120071). 

Spin mixing oscillations in both ferromagnetic ( 87 Rb) 
and antiferromagnetic ( 23 Na) F = 1 spinor Bose-Einstein 
condensates have been observed. The Chapman group 
initiated spin-mixing oscillations of 87 Rb in the inter- 
action regime by preparing condensates in spin states 
that were neither ferromagnetic nor polar, making use 
of microwave transitions between hyperfine levels. From 
time- of- flight based measurements of po? the consequent 
trajectory in spin space was revealed to undergo sev- 



eral complete revolutions through pQ-6 phase space be- 
fore the system rela xed to the predicted ground state 
(Chang et al. 2005). The Hamburg group focused on 



oscillations in the Zeeman regime, using it/ 2 rf pulses 
to rotate the atomic spin prior to Stern-Gerlach anal- 
ysis, thereby permitting t he observation of orientation- 
to-alignment oscillations ( [Kronjager et al 2006). Fi- 
nally, working with 23 Na condensates, the NIST group 
employed both Stern-Gerla ch anal ysis and also disper- 
sive optical detecti on (Sec. IV.C.2[) to measure both p p 



and (9, respectively ( |Black et al. H2007||Liu et al. [|2009bp . 



Altogether, these measurements enabled a determination 
of the spin-mixing energy from dynamics in the interac- 
tion regime (Black et al] 2007 Chang et al. 2005), and 



the observation of the spi n-mixing resonance at nonzero 
quadratic Zeeman shifts (|Black et al. , 2007 Kronjager 



et al, 2006 Liu et al, 2009b). 



Spin mixing dynamics in F = 2 spinor gases are richer 
owing to the larg er phase space spanned by such dynam- 
ics. As noted by Kronjager et al. (2008), in addition to 



the aforementioned gauge and rotational symmetries, the 
equations of motion for the spin state are invariant under 
rotations by tt about a transverse spin axis, with the sub- 
stitution m —> —m. For the F = 2 case, this additional 
non-trivial symmetry reduces the number of dynamical 
variables to four. Deep in the interaction and Zeeman 
regimes - the limits of \q\ being either much smaller or 
much larger than the spin-dependent interactions - the 
spin-mixing again shows simple oscillatory behavior, as 
confirmed both by experiments and by theoretical ap- 



proximations (Kronjager et al, 2006 



Kuwamoto et al. 



20041) . In the interaction regime (small |g|), the specific 
trajectory of spin- mixing dynamics depends on the rela- 
tive strengths of the spin-dependent interaction parame- 



ters and c 2 



(2) 



providing an experimental means to de- 
termine such relative strengths clearly (|Saito and Ueda| 



2005 ) . Intriguingly, between the interaction and Zeeman 



regimes, where one obtains a spin-mixing resonance in 
the F = 1 case, here the mean-field theory predicts the 
onset of chaotic motion at the mean- field level (Kronjager 
2008]). 



et 



Finally, we mention several studies of spin-mixing dy- 
namics in dr i ven spinor g as systems. Pu et al. (2000) and 
Zhang et al. (2002 2001 ) considered theoretically the sta- 



bility of spin mixing oscillations in the presence of exter- 
nal field noise to assess whether experimental limitations 
on such noise would bar the observation of coherent spin 
mixing dynamics. It will be interesting to revisit such cal- 
culations in light of the prediction for chaotic dynamics in 
high-spin spinor gases. They also assessed the response 
of the spinor gas to a periodic magnetic fields, finding 
regimes where the conventional Rabi oscillations of non- 
interacting spins give way to dynamical spin localization, 
akin to the physics of Jose phson junctions (see also Ya- 
sunaga and Tsubota ( 2010b[ )). Dynamical localization is 



also predicted in the case that the spin- d ependent scat- 
tering length is periodically modulated (Zhang et al.\ 



2010). Pulsed magnetic excitation is predicted to lead to 
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chaotic spin dynamics, even in the F = 1 case (Cheng 
2010), similar to quantum ch aotic phenomena observed 
in nthpr nnlspH snin svstpms fNakamnra 1QQ4K 



3. Many-body single-mode dynamics 

For F = 1 spinor gases, a many-body single- spatial 
mode treatment of spin-mixing dynamics at q = is ob- 
tai ned directly from the exact energy s pectru m derived 
by Law et al (1998), as explained in Sec. V.B Indeed, an 



analytic Heisenberg-picture solution for spin-mixing dy- 
namics base d on this spectrum has been obtained ( |Chen 
et al 2008). To compare with mean- field theory, let us 
consider the specific case of spin mixing dynamics for a 
state prepared initial ly nea r a ferromagnetic state. As in 
experiments ( Chang et al\ 2005), we consider the evolu- 
tion of the Zeeman populations in the \m z = 0) state. 
Such populations, considered above at the mean-field 
level, are determined from the many-body spin states as 
the expectation value of a sum of single-atom projection 
operators. This sum is an operator of rank 2, so the tem- 
poral oscillations of the \m z = 0) state population occur 
at the frequency determined by the energy difference be- 
tween many-body states of total angular momentum L 
and L — 2: 



n N - 



(1) 

c\ J n 
2HN 



(4L-2) 



(69) 



This re sult reproduces both the two-atom result (Sec. 
IX. A. 1) and also, for large atom number N and with 



L = N near the ferromagnetic initial state, the mean- 
field result discussed above (Sec. IX.A.2). 

At longer evolution times, the anharmonicity of the en- 
ergy spectrum becomes evident (|Law et al 1998). Con- 



sider a state initialized in a coherent spin state that, if 
substantially different from the ferromagnetic state, rep- 
resents a superposition of angular momentum eigenstates 
spanning a range AL ~ y/~N. Spin mixing dynamics now 
involve a range of Bohr frequencies spanning a band- 
width of Auj ~ Sl^TV -1 / 2 . After an evolution time 
t > (Acj) _1 , the coherent spin- mixing oscillations pre- 
dicted by mean-field theory are lost to "quantum phase 
diffusion," similar to that discussed in the physics of 



scalar Bose-Eins t ein condensates ( Lewenstein and You 
1996} |Pitaevskii[ |1997| [Wright e£ al.| |1996|). However, 



given that all Bohr frequencies remain integer multiples 
of a common factor, even beyond the collapse of mean- 
field dynamics, the exact many-body dynamics neverthe- 
less exhibit remarkable structure as seen in the predicted 
probability distributi on for measurem ents of Zeeman- 
state populations (Diener and Hoi 2006). The experimen- 



tal observation of such many-body effects awaits the cre- 
ation of well-controlled mesoscopic spinor Bose-Einstein 
condensates. 



B. Spin mixing instability 



In the interaction regime, the energy functional of Eq. 



(68) contains also energy maxima. As a specific example, 



.a) i 



for a ferromagnetic F = 1 system, with < q < 2\c\ 
the high-energy state with p = 1 (inert state) is, at 
the mean-field level, stationary. Whereas the mean-field 
ground state has a non-zero magnetization oriented in 
the transverse spin plane, the p = 1 stationary state 
does not select a transverse spin orientation. 

Experiments and theory have shown that this mean- 
field stationary state is, in fact, dynamically unstable. 
Akin to a marble placed at the apex of a convex potential 
landscape, the collective spin state will diverge exponen- 
tially from the unstable stationary state on a trajectory 
initiated by even microscopic fluctuations. Starting from 
an unbroken symmetric state, identified by zero- valued 
expectations of the relevant magnetic moments, the sys- 
tem evolves toward states of strong non-zero transverse 
magnetization and broken symmetry. The relevance of 
such dynamics to topics such as quantum noise, quantum 
amplification, symmetry breaking and phase transitions 
has led to many studies of the spin mixing instability 
in both single-mode and spatially extended spinor Bose 
gases. 

To discuss this full range of topics, let us consider the 
dynamics of the inert \m z = 0) condensate under the 
Bogoliubov approximation, retaining terms up to second 
order in the fluctuations atop the condensate wavefunc- 
tion. We will find it convenient to work in the Carte- 
sian or polar-state basis, with ijj a being the \m a = 0) 
eigenstate for the F a spin operator (a G {x,y,z}). The 
spin-dependent Hamiltonian is approximated as follows: 



H 



P={x,y} 



d 6 r$l (#o + 4 1} n(r))^ 



n(r) 



($ 2 (r)- 



- 2 

(r) 



(70) 



Here we account for spatial variation of the spinor con- 
densate: n(r) is the density of the initial-state conden- 
sate, and x/jp is the field operator for particles in the ini- 
tially unoccupied polar states labeled by j3. Such par- 
ticles experience the single-particle Hamiltonian of the 
form Ho = — h 2 V 2 /2m + V(r) + q where, for simplic- 
ity, we assume all particles experience the same potential 
V(r). 

Now, let us simplify the problem by considering the dy- 
namics of a single spin-excitation (magnon) mode, with 
two polarizations as labeled by j3. Such a mode can 
be identified either as a discrete magnon mode in the 
single-spatial-mode regime discussed above, as a linear 
combination of momentum eigenstates in the case of a 
homogeneous condensate, or by careful consideration of 
the appropriate Bogoliubov-de Gennes equations for this 



problem (Baraban et al 2008 Sau et al 2009 2010). 



Letting ipp now be the creation operator for particles 
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in this magnon mode, the dynamics obey the following 
Hamiltonian: 



H 



P 



where the energy e gives the kinetic and potential energy 
of the magnon mode. 

Following the analogy of a marble rolling downhill, let 
us define Hermitian operators Zp = + $p )/2 an d 
Pp = — $p)/2, and now obtain (Lamacraft 



2007) 



^ = E[( 6 +^) i l + ( 6 +^ +2c i 1)n )^ 



(72) 



In the case that the coefficients of Zi and Pj have the 
same sign, as obtained in the Zeeman regime when \q\ 
dominates the other energy scales, each polarization of 
the magnon mode evolves stably, precessing in the phase 
space spanned by the Zp and Pp operators at a frequency 

given as u = VEP/h with E 2 = (e + q)(6 + q + 2c{ ] n) > 0. 
In contrast, when the coefficients have opposite sign, 
which occurs when e + q lies between and — 2c^n, 
the magnon mode represents a dynamical instability of 
the initial-state condensate. While (Zp) = (Pp) = 
in the initial state, the fluctuations of these operators 
become amplified parametrically: one quadrature of the 
Zp — Pp phase space grows exponentially at the rate 
uj = +V—E 2 /h, while the other quadrature is attenu- 
ated. 

This dynamical instability was noted first by a high 
sensitivity to noise in numerical simulations of spin mix- 
ing dynamics (Pu et al., 1999), prompting investiga- 



tions of the effects of field fluctuations on such dynam- 



ics (Pu et al. 2000). A more focused study by |Robins] 



et al. (2001 ) was informed by the analogy between coher- 
ent matter-wave dynamics and optical four-wave mixing 
(Deng et aZ.[ |1999|) . Robins et al. examined the evolution 



of spatially uniform inert states of an F = 1 condensate - 
either the polar state (Case 1 in their work) or the mag- 
netized state (Cases 2-4), in the absence of a quadratic 
Zeeman shift. A modulational instability is observed only 
for ferromagnetic interactions and for the polar state, at 
modulation wavevectors k for which e = e& = h 2 k 2 /2M 
gives E 2 < as defined previously. The instability was 
predicted to cause an exponential rise in the population 
of previously unoccupied Zeeman states, and strong spa- 
tial modulation of the condensate spinor wavefunction. 

The predicted transfer of Zeeman populations was ob- 
served in the first experiments on dyn amics of ferromag- 
netic F = 1 condensates QChang et al. , 2004). A conden- 
sate prepared in the \m z = 0) state at low q evolved 
slowly for several hundred milliseconds before rapidly 
tending toward a mixture of Zeeman states. Similar dy- 
namics were observed also for antiferromagnetically in- 
teracting spinor gases (c^ > 0) at negative quadratic 
shifts. For example, the stability analysis for a 87 Rb 



F = 2 condensate prepared in the \m z = 0) initial 
state is similar to that of the F = 1 antiferromagnetic 
gas. Such condensates were found to undergo collision- 
induced spin mixing at low dc magnetic fields, for which 
q < in the F = 2 87 Rb spin manifold. The dynamics 
are substantially faster than that for the F = 1 manifold, 
owing to the larger spin-dependent interaction strength 



(Kuwamoto et a/.[|2004| Schmaljohann et al. 2004b. The 



spin-mixing instability was also demonstrated recently 
for the antiferromagnetic F = 1 23 Na condensate, for 
which off-resonant microwaves were used to induce a neg- 



ative quadratic Zeeman shift ( |Bookjans et al.\ 2011b). 

Spinor Bose condensates display other instabilities as 
well. Returning to the single-mode energy functional for 
the F = 1 system (Eq. (68)), the energy extremum de- 
fined by = (dubbed the "phase-matched state" in 

q/(2c[ 1) n)] /2, 



Matuszewski et al. 



(|2008)) and p = 



which exists for \q\ < \2c\ ~ y n|, is an energy maximum un- 
der antiferromagnetic interactions (denoted by the black 
line in Fig.[5jc)). A wave-mixing analysis shows that this 
stationar y state is also subject to modulational instabil- 
ity ( |Matuszewski et a~ 2008) that can produ ce intricate 
spin domain structures (Matuszewski, 2010). Modula- 



tional instability was observed experimentally in a sim- 
ilar system - the transversely mag netized F = 2 87 Rb 
spinor condensate - generating spin domain structures 
with a periodicity controlled by the quadratic Zeeman 



shift (Kronjager et al. 2010) (Fig. 18). 



C. Quantum quench dynamics and the Kibble-Zurek 
mechanism 

Considering the \m z = 0) energy extremum for the fer- 
romagnetic F = 1 single-mode spin energy, we observe a 
transition from dynamical stability to instability coinci- 
dent with the state's transition from an energy minimum 
to a maximum, i.e. with a second-order zero-temperature 
phase transition. The evolution of systems prepared ini- 
tially at the extremal state thereby exemplifies the dy- 
namical response of a system quenched rapidly across a 
phase transition, evolving from a state of unbroken sym- 
metry (here, SO (2) axisymmetry) toward a manifold of 
broken- symmetry ground states. 

In classical systems rapidly quenched across a finite- 
temperature phase transition, symmetry-breaking occurs 
inhomogeneously, with different portions of the system 
adopting different symmetry-broken phases. Such a pro- 
cess was considered by |Kibble| ( |1976[ ), in the context of 
symmetry breaking in the early universe, and by |Zurek 
(1985), in the context of low-energy laboratory exper- 
iments on material systems. Their treatment, known 
as the "Kibble-Zurek mechanism," translates the phe- 
nomenology of critical scaling of equilibrated systems 
near phase transitions into predictions regarding the typ- 
ical size of the symmetry-broken phases. Their theory 
also discusses the types of topological defects produced 
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FIG. 18 (Color) Spontaneous pattern formation from modu- FIG. 19 (Color) Spontaneous symmetry breaking in a fer- 
lational instability of a transversely mag netized F — 2 87 lib 
spinor condensate. The condensate is prepared in a quasi one- 
dimensional optical trap. The initial state (a) and saturated 
spin patterns arising for increasing magnetic field (b-f) are 
determined from Stern- Gerlach/time-of- flight images that re- 
veal the in-trap axial density distributions of the five Zeeman- 
state populations. Cases (b) and (c) represent the interaction 
regime (small q), while cases (e) and (f) represent the Zeeman 
regime (large q). The characteristic length scale of the spin 
modulation varies stro ngly with q. Figure reproduced from 
|Kronjager et a/.| ( |2010| . 



in such a quench, their initial density (related to the 
aforementioned typical size), and their subsequent evo- 
lution. Laboratory experiments have explored aspects 
of the Kibble- Zurek mechanism, e.g. using liquid crys- 



tals ( Bowick et a/. , 1994) Chuang et a/. , 1991), pressure- 
quenched 4 He (Hendry et al. 1994 ), neutro n-bombarded 
3 He (|Bauerle et a/.[|1996||Ruutu et al. jl996[), and Joseph- 



son junctions (Carmi et al.\ 2000 Monaco et al. 2002). 



Given theoretical progress in understanding quantum 
phase transitions, and experimental progress in preparing 
isolable low-temperature quantum systems, the Kibble- 
Zurek idea was naturally extended to the quench of a 



system across a quantum phase transition (Dziarmaga 



2005 Polkovnikov, 2005 



Zurek et . 



2005). For an ideal 



experiment, the quench dynamics are now quantum me- 
chanical, in that quantum fluctuations (rather than ther- 
mal or technical ones) are amplified into coherent super- 
positions of macroscopically distinct, broken symmetry 
states. Ultracold atomic gases are unique materials with 
which to study such ideal quantum quenches. 

The evolution of a quenched spinor Bose-Einstein con- 



romag netic F — 1 Rb spinor Bose-Einstein condensate 
quenched across a phase transition. A condensate is prepared 
in the \m z — 0) state at high q, after which q is suddenly re- 
duced to q < \c^\no where spin-dependent interactions favor 
a transverse ferromagnetic state (no is the central column- 
averaged condensate density), (a) The transverse column- 
integrated magnetization, measured for condensates at vari- 
able times after the quench (indicated on top) , is shown with 
the magnetization orientation (0) indicated by hue, and mag- 
nitude (A) indicated by brightness. The maximum brightness, 
indicated by the color wheel at left, corresponds to full mag- 
netization of the condensate center. For the data at 216 ms, 
the magnetization density (b) and orientation (c) are shown 
separately. The spin texture at position I is characterized by 
an extended region with magnetization of non-zero magnitude 
and slowly varying orientation. Smaller magnetic domains at 
position II are divided by domain walls with zero magnetiza- 
tion. The gray scale in (c) covers the range to 2ty. Regions 
outside the condensate are indicated in black. Figure repro- 
duced from |Sadler et a/.| ( |2006| ). 



densate was studied experimentally by Sad ler et al 
(120061). A large F = 1 87 Rb spinor condensate, spa- 
tially extended in two dimensions, was prepared in the 
\m z =0) initial state at a high quadratic Zeeman shift. 
Following a rapid reduction of the quadratic Zeeman 
shift, which initiated the spin- mixing instability, the gas 
evolved for a variable time before it was probed by disper- 
sive magnetization-sensitive imaging. The gas showed lit- 
tle magnetization for tens of milliseconds, after which an 
inhomogeneous magnetization landscape spontaneously 
emerged (Fig. [l9|). 

For the sudden quench studied in this experiment, the 
spatial spin correlations of the inhomogeneously mag- 
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netized condensate are determined by, and are, there- 
fore, reflective of, the spin-mixing amplification spectrum 
(Lamacraft 2007). Considering the simplest case of a 



homogeneous initial state, at early times after the on- 
set of the spin-mixing instability, the evolution of the 
initial spin fluctuations at each wavevector can be con- 
sidered independently according to Eq. (72), where the 



temporal evolution is defined through the magnon dis- 
persion relation E 2 (k, q) = (e& + q)(ek + q + 2c^n). For 
a quench to q c± 0, this highest temporal gain occurs at 
the wavevector /c max = l/£ sp , leading to domains with 

dominant size 7r£ sp where £ sp = h/ \J 2mc^n is the spin 
healing length for this system, and with time constant 
of r = H/(c^ri). These predictions were in rough agree- 
ment with experiment. However, details of the spatial 
spin correlation function were une xplained by the pu rely- 
quantum, homogeneous theory of Lamacraft (2007). Re- 



fined theories were constructed that account explicitly 
for the anisotropy of the trapped conden sate by defining 
the p roper magnon modes for that case ( |Baraban et al] 
2008). Numerical simulations obtained agreement with 



measurements b y varying t he spatial Fourier spectrum of 
the initial noise (Saito et al. 2007b), perhaps indicating 



technical imperfections in the initial state preparation. 

Following this work, several theoretical papers con- 
sidered the realization of the Kibble-Zurek mechanism 
in spin-ordering transitions of the ferromagnetic F = 1 
spinor Bose-Einstein condensate 



|Lamacraft 



(|2007J) dis- 
n) in which 



tinguished between deep quenches (q < \c' m 
spin domains will emerge at a length scale determined by 
the depth of the quench, and shallow quenches (|c^n| < 

q < 2\c^\n) for which spin correlations will show "light- 
cone" dynamics and expand linearly in time. The latter 
behavior is reminiscent of the Lieb-Robinson expansion 
of quasiparticle correlations and entanglement predicted 
for quantum quenches that do not necessarily cross phase 



transitions (Bravyi et al. 



2006| Calabrese and Cardy 



2006 Cramer et al, 2008). Gradual linear quenches 
across the phase transition were considered, essentially 
by considering evolution under the magnon dispersion 
E 2 {k 1 q(t)) with q(t) varying linearly in time, and scal- 
ing laws for the characteristic domain length and spin- 



vortex density were obtained ( 


Damski and Zurek, 2007 


Lamacraft[ |2007| Saito et al. 


2007a| Uhlmann et al. 


2007). Saito et al. (2007a) noted that the spin- mixing 



instability has the special feature of local spin conser- 
vation; while this conservation requires that regions of 
spin anti-cor r elatio n exist at short range (as observed by 
|Sadler et al.\ ( |2006| )), the standard Kibble-Zurek mecha- 
nism still applies over longer length scales. 



plored this question by testing whether the macroscopic 
magnetization patterns observed after the quench could 
indeed be ascribed to the quantum-limited parametric 
amplification of quantum noise. In the spatially extended 
systems which they studied, the spin-mixing instability 
is represented as a spatially broadband parametric am- 
plifier, i.e. one with a wavelength dependent amplifica- 
tion, the spatial spectrum of which is defined by the 
magnon dispersion relation E 2 (k,q). That amplification 
spectrum is detected experimentally by quenching to dif- 
ferent values of g, feeding spatially white noise at the in- 
put of the amplifier and detecting the post-amplification 
output via the spatial Fourier spectrum of the magne- 
tization. The amplification spectrum was indeed ob- 
served to vary with q, with shallow quenches amplify- 
ing longer wavelength features, while deeper quenches 



produced small length scale spin domains (Leslie et al. 



2009b Sau et al., 2009). However, the experiments did 
not verify the quantum nature of the seeding and ampli- 
fication due to uncertainties o f the spin- dependent inter- 
action energies for 87 Rb (Chang et al] 2005 van Kemp en 



et 



2002; Klausen et all 2001 Widera et all 2006). 



Stronger evidence for fully quantum mechanical evo- 
lution has been obtained in single-spatial-mode experi- 
ments, performed with spinor Bose gases in tightly con- 
fining optical traps. The Hannover group studied the 
spin-mixing instability of the \m z = 0) state for the 
F = 2 87 Rb spinor condensate, in a situation where 
discrete spin-excitation modes were separated energeti- 
cally by gaps, on the order of the transverse trapping 
energy, that were made larger than the spin-mixing en- 
ergy. By tuning the quadratic Zeeman shift to match 
the excitation energies e, the spin-mixing instability led 
to the amplification of discrete spin-excitation modes for 
which the amplification gain was largest. Comparing the 
Zeeman populations generated after a fixed duration of 
spin mixing for the various magnon modes, they con- 
cluded that higher-order spatial modes were less contam- 
inated by technical noise, such as that due to imperfect 
purification of the initial spin state, potentially to the 
level where the initial fluctuations of those higher-order 
modes was due to vacuum fluctuations alone (Klempt 



et al. 



2009, 2010). In their cylindrically symmetric trap, 
spin-excitation modes can be labeled by radial and angu- 
lar quantum numbers, with the states (n, ±1) being de- 
generate. Remarkably, they observed the amplification 
of spin fluctuations selected a coherent superposition of 
such states which varied between runs of the experiment 
(Fig. [20]), representing a spontan eous breaking of spi n 
and spatial rotational symmetries (Scherer et al, 2010). 



D. Parametric spin amplification 

In describing the experiments discussed above as quan- 
tum quenches, what is meant s pecifically by the w ord 
"quantum?" Subsequent work by Leslie et al. ( 2009b ) ex- 



E. Quantum spin-nematicity squeezing 

Distinctly quantum-mechanical evolution was observed 
in recent experiments on spin squeezing generated by the 
spin mixing instability. The concept of spin squeezing 
has been derived most extensively with reference to the 
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FIG. 20 For a spinor condensate trapped in the single spatial- 
mode regime, the spin mixing instability can be tuned to dis- 
crete spin excitation modes by varying the quadratic Zee- 
man shift, (a) The fraction of atoms transferred into the 
\m z = ±1) sublevels from the initial \m z = 0) state, within 
18.5 ms of the onset of dynamical instability (determined 
from images such as (b)), shows a mult i- resonant structure; 
here, experiments are performed with F — 2 87 Rb atoms 
for which q is negative for dc magnetic fields. The la- 
bels (n, I) give the radial and angular quantum numbers of 
transverse spatial output modes. Resonances occur when 
q matches the output mode energy (theoretical predictions 
shown by gray bars). The spatial structure of the output 
modes matches the predicted patterns for averages over many 
realizations (c:experiment, d:theory), and for single realiza- 
tions (eiexperiment, f:theory). The single realizations show 
spontaneous breaking of both spin and spatial symmetry as 
particular superpositions of degenerate output modes are se- 
lected at random. Figures reproduced from Sch erer et al.\ 

( poTo| . 



collective spin state of an ensemble of spin- 1/2 particles. 
In the absence of interparticle correlations, such an en- 
semble is described by the single-body density matrix 
and, thereby, defined by its net vector spin F to t = (Ftot) 
where F to t is the symmetric sum of vector spin opera- 
tors over the ensemble's particles. Allowing the orien- 
tation of F t ot to define the z axis, the uncertainty rela- 



tion ((AF x ) 2 )((AF y ) 2 ) > (F z ) 2 /4 is satisfied by parti- 
tioning the uncertainties equally between the two trans- 
verse spin components. Correlations among the particles 
allow one to apportion the uncertainty unevenly, allow- 
ing the measurement uncertainty of one transverse spin 
component to be lower than the standard quantum limit 



/SQL 



Wineland et al. 



\{F z ) \/2 ( [Kitagawa and Ueda[ |1993 
1992). Such squeezed states can pro- 
vide sub-quantum-limited performance for atom-based 



sensors such as atomic clocks ( Appel et al. , 2009 Schleier- 



Smith et al. 2010), magnetometers (|Koschorreck et al] 



2010), and interferometers (Gross et all 2010) 



Spin squeezed states can be generated naturally by 
the evolution of an ensemble under a system Hamilto- 
nian that contains terms quadratic in the spin operators 
(Kitagawa and Ueda, 1993). Such terms are generated 
by binary collisional interactions, e.g. in Bose-Einstein 
condensates, as pointed out in Sorensen et a/.| ( |2001[ ) and 



recently realized for a pseudo-spin- 1/2 condensate in a 



state-dependent potential (Riedel et al. 2010). Natu- 



rally, researchers suggested that spin squeezing would be 
generated by the spin-dep endent interaction s in spinor 
Bose-Einstein condensates (Duan et al., 2002). However, 



the spin- 1/2 theory of vector spin squeezing is clearly in- 
adequate for describing squeezing in higher-spin spinor 
Bose gases, since the single-particle density matrix is 
characterized by more observables than just the vector 
spin. 

Several frameworks have been su ggested for describ- 
ing spin squeezing in this context. |Duan et al.\ ([2002 ) 



suggest that the definition of spin squeezing for spin- 1/2 
systems can be adapted to the higher-spin case by select- 
ing any two orthogonal single-particle states, or, equiva- 
lently, any two orthogonal bosonic modes based on those 
single-particle states, and denoting them as pseudo-spin 
+1/2 and -1/2 states, respectively. 



A similar approach is taken by Sau et al. (2010), who 
point out that for any coherent spin state of a spin-F en- 
semble, in which all particles occupy the single-particle 
spinor |^), one can select 2F spinors that form an or- 
thonormal set with Based on these, one defines spin 

fluctuation observables m£ ] = V. (|f a )(^| + h/>)(fa|) 7 - 



J 3 VIW/ \ t I ' IT/ \^ U J 

and Mi = J2j i (IfaX^I _ where the sum is 

taken over all particles in the ensemble. These observ- 
ables, each of which has an expectation value of zero and 
variance equal to the standard quantum limit for the co- 
herent spin state |^), represent 2F independent polariza- 
tions of spin fluctuations, and play the role of transverse- 
spin operators for independent pseudo- spin- 1/2 systems. 

For a spinor Bose condensate prepared initially in the 
|^) = \m z = 0) spin state, choosing the orthogonal 
states {|£q;)} = { \ip x ), l^y}} from the polar- state ba- 
sis, one finds the two polarizations of spin-fluctuation 



observables to 



r(i) 



~N yz ,My 



be (m< 
(2) --^). 



-N XZ ,M£ 



(2) 



and 



F x ) . That is, in contrast to 



the vector spin squeezing familiar for spin-1/2 particles, 
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here, spin squeezing implies the reapportionment of fluc- 
tuations between conjugate spin- vector (F x and F y ) and 



spin-quadrupole components (N yz and N xz ) (Leslie et al. 



2QQ9b| Sau et al. , 2Q1Q[ ) . More recently, a rigorous set of 
criteria for squee zing and entanglement in high- spin gases 
was provided by Vitagliano et al. (2011). 



Such spin-nematicity squeezing has been recently ob- 
served (Fig. 21). An \m z = 0) condensate prepared in 
a tight optical trap was evolved under the spin mix- 
ing instability for several spin- mixing times. To observe 
the squeezed quadrature, the gas was first exposed to a 
pulsed quadratic Zeeman shift, which rotates spin fluctu- 
ations in the spin-nematicity phase space, and then the 
transverse spin was measured by applying a tt/2 spin ro- 
tation and then measuring F z via Stern-Gerlach analysis. 
Significant squeezing, at the level of around 8 dB below 
the standard quantum li mit, was observed w hen probing 
the squeezed quadrature (Gross et al. 2011). 



Sim ilar experiments were conducted by |Gross et 



(2011), working with a sample of several hundred con- 
densed atoms in the F = 2 manifold of 87 Rb. They inter- 
pret the measurement of coherences between the initial 
\m z = 0) state and the \m z = ±1) states produced by 
the dynamical instability as a form of homodyne mea- 
surement akin to those performed in quantum optics. 
However, sub-shot-noise squeezing was not detected in 
this work due to technical limitations. 

Spin-nematicity squeezing records the coherence gener- 
ated between the particles emitted by spin mixing (in the 
\m z ± 1) states) and the condensate from which they are 
produced. Alternately, one may focus on the correlations 
generated between the emitted particles themselves. By 
spin conservation, the populations generated by spin mix- 
ing in the \m z = ±1) states should be strictly identical; 
indeed, experimental measurements of such populations 
have shown sub-Poissonian, resolution-limited differences 



in their populations (Bookjans et al. , 2011a Chang et al. 



2004 Liicke et al. 2011). However, such number corre- 



lations do not, on their own, imply any quantum cor- 
relations or entanglement. Rather, to detect such cor- 
relations, one probes whether the particles emitted by 
spin mixing inhabit two (one for each spin state) common 
bosonic modes. Under such conditions, the \m z = ±1) 
atomic populations are represented as a superposition of 



twin Fock states (Duan et al, 2000 Mias et al. 2008 



Pu and Meystre| |2000fl , |N + ,lV_77with N+ = N_, simi- 
lar to the photonic states produced by parametric down- 
conversion. 

One feature of this twin Fock state is that the variance 
in the population difference between the two modes is 
initially minimal, but grows rapidly if one first applies a 
coherent rotation by an angle <p between the two modes 
(similar to the action of an optical beam splitter on light). 
This feature was observed in recent experiments by the 
Hannover group, who point out that the rapid increase 
in measurement variance gives a quantum-enhanced sen- 
sitivity to the rotation angle (j) dLiicke et a/.[|2011[). An 
additional feature, pointed out by|Mias et al. (2008), is 



that, in the twin-Fock-state superposition, the popula- 
tions N+ = N- should show thermal-like fluctuations 
even though the emitted particles inhabit a zero-entropy 
quantum state. Such fluctuations, analogous to Hawking 
radiation, have not yet been studied experimentally. 



X. COARSENING AND EQUILIBRATION 

The equilibration of a quenched system, such as a fer- 
romagnet quenched suddenly below its Curie point or a 
binary alloy undergoing phase separation, consists gen- 
erally of two stages of evolution: an initial stage dom- 
inated by the fastest timescales for establishing locally 
ordered regions, and an ensuing coarsening stage during 
which the order parameter becomes ordered on length 
and time scales that grow ever longer. In degenerate 
spinor Bose gases, the former evolution is described well 
in terms of the coherent quantum dynamics of collective 
spin mixing following a quantum quench. The nature of 
the latter stage of evolution, the long-time relaxation of 
spin-mixing oscillations and the phase-ordering kinetics 
of spin textures, remains largely unexplored. Here, we 
describe several aspects of relaxation and equilibration. 



A. Relaxation of spin-mixing oscillations 

Spinor Bose-Einstein condensates prepared in non- 
equilibrium spin states undergo coherent spin-mixing dy- 
namics, marked by the deterministic initial evolution 
both of the mean-field spinor wavefunction and its many- 
body quantum fluctuations. At longer times, however, 
the signatures of this coherent evolution are lost. This 
behavior is exhibited most clearly in two experimental 
reports. Studying spin- mixing dynamics of an F = 1 
87 Rb condensate, the Chapman group observed coherent 
oscillations of the Zeeman populations that were damped 



out after a few oscillations (|Chang et al. 2005 ) . The ob- 



served damping coincided with the appearance of spatial 
structure along the long axis of their prolate gas samples. 
In those experiments, the comparison between the con- 
densate radii of (3.2, 3.6, 36) /im and the expected size 
of spin domains 7r£ sp = 9 fim supports the notion that 
the observed damping can be explained as resulting from 
inhomogeneous broadening of the spin-mixing frequency 
along the long axis of the gas ( |Mur- Petit et al. , 2006). 

A second experiment, by the NIST group, probed the 
spin-mixing oscillations of F = 1 23 Na condensates ( |Liu| 
\et al.\ [2009a). For that work, the condensate was pre- 
pared with spatial dimensions (radii ranging from 6 to 8 
/im) roughly equal to the expected size of spin domains 
(^Csp = 8/im), so that the single-mode approximation 
should apply. Experiments were conducted with a small 
quadratic Zeeman shift, for which the mean- field spin- 
mixing dynamics divided between the closed-orbit tra- 
jectories (in pQ-0 space; see Sec. IX. A. 2) of the higher- 
energy spin states, and the running trajectories of the 
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FIG. 21 (Color online) Observation of spin-nematicity squeezing. An F = 1 87 Rb condensate in the single-mode regime is 
prepared in the \m z = 0) initial state and allowed to evolve under a spin-mixing dynamical instability. The gas is then probed 
by first using a brief pulse with variable q (controlled by a microwave drive with detuning shown on the top axis) to rotate the 
spin state in the transverse spin-nematicity plane (by the quadrature phase shift A0), and then using a tt/2 spin rotation and 
Stern-Gerlach analysis to measure one component of the transverse spin. (Left) Choosing AO to measure the quadrature that is 
squeezed by the dynamical instability, the measurement variance is found to be below the atomic shot-noise level by as much as 
8.6 dB. (Right) From measurements at all quadrature angles one reconstructs the Wigner function in the spin-nematicity phase 
space produced following (a) 15, (b) 30, (c) 45, and (d) 65 ms of parametric amplification by the spin-mixing instability. The 
squeezed and anti-squeezed quadratur es are clearly observed. Black ellipses indicate the 1/y/e measurement width calculated 
from theory. Figures reproduced from Hamley et al. (2012). 



lower-energy spin states. These two dynamical regions 
are divided by a separatrix marked by larg e, slow varia- 
tions in the \m z = 0) Zeeman population (Black et al. 



2007). 



The spinor gases were prepared initially in the high- 
energy, transverse ferromagnetic state. The ensuing spin- 
mixing dynamics, monitored by observing the Zeeman 
populations, indicated that the spin system underwent 
"orbital decay," evolving from high- to low-energy spin- 
mixing trajectories. The authors modeled this decay 
phenomenologically by adding a damping term to the 
mean-field dynamical equations. However, the true na- 
ture of this damping remains unclear. En do and N ikuni 
(2011) have provided a framework for describing dissipa- 



tion of spin mixing dynamics through damped collisional 
dynamics of the non-condensate fraction. However, a 
quantitative prediction has not yet been derived from 
this framework to test its applicability to the experimen- 
tal situation. 



B. Elements of coarsening dynamics 

Far outside the single- mode regime, the equilibration 
dynamics of a spinor gas that is quenched across either a 
quantum or thermal phase transition may be explained 
using the theories of phase-ordering kinetics and dynam- 
ical critical phenomena that hav e been developed to de- 
scribe a wide variety of systems (Bray 2002 Hohenberg 
and Halperin 1977[ ). Elements of such theories include 
the hypothesis that domain structures and correlation 
functions at different times in the equilibration process 



differ only in the overall length scale, and that this length 
scale grows in time by the power law l(t) ex t x l z , with z 
being a dynamical critical exponent. Such scaling has 
been examined numerically, and, to a limited extent, ex- 
perimentally in several models of relaxation dynamics 
that differ in the symmetry of the order parameter and 
in the presence of conserved quantities. Scaling dynam- 
ics can often be related to the behavior of topological 
structures, such as vortices and domain walls. In this 
sense, studies of topological structures in spinor Bose 
gases should shed light on their coarsening dynamics. 

The mechanisms of equilibration in spinor Bose gases 
were explored in experiments on the appearance of sharp 
domain walls in F = 1 spinor condensates of sodium, 
and on their role in impeding phase ordering of the gas 
( |Miesner et~al\ |1999| |Stamper-Kurn et al] |1999|). Rf 



pulses were used to prepare a condensate in a superpo- 
sition of the \m z = 0) and \m z = 1) Zeeman sublevels, 
at a high value of the quadratic Zeeman shift where spin 
mixing into the \m z = —1) state is energetically sup- 
pressed so that the number of atoms in each Zeeman 
state is conserved. In accordance with the phase dia- 
gram of Fig. |5j with the relation aoi > y^oo a n be- 
tween s-wave scattering lengths in the case of antiferro- 
magnetic F = 1 interactions, the initially overlapping 
states underwent rapid phase separation (Fig. 22). Such 
dynamics can be described as due to a dynamical insta- 
bility that is most dominant at the wavevector equal to 



the inverse spin healing length (Goldstein and Meystre 
19971 IGraham and Walls] [1998 



Isoshima et all 11999 



Timmermans 1998). This experiment represents another 



: 'quantum quench," in which a system is prepared in the 
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non-equilibrium admixed state under conditions that fa- 
vor the breaking of an Ising symmetry by the process of 
phase separation. The observed dynamics resemble the 



(Sadler et al. (2006); see Sec. IX. C) and antiferromag- 



classical spinodal decomposition of binary alloys ( Ao and 



Chui||1998||2000| ). 



The condensates in these experiments were extremely 
prolate, so as to be effectively one-dimensional with re- 
spect to spin dynamics. Accordingly, the textures formed 
after the quench showed only axial variations in spin 
composition, comprised of single-component domains of 
length ~ 40 fim separated by sharp domain walls. 

In the case of a conserved order parameter, classical 
phase ordering proceeds by diffusion, which drives parti- 
cles from the surface-tension-induced high chemical po- 
tential of small domains to the lower chemical potential 
of larger domains ( |Bray| [2002 ). In the phase-separated, 
quantum-degenerate Bose gas, such diffusion describes 
the evaporation of condensed atoms from one portion of 
the gas, and the recondensation of the normal fluid into 
another portion of the gas (where the chemical potential 
is lower). The quantum-degenerate gas also features a 
second means of transport: quantum tunneling of the su- 
perfluid across the potential barriers erected b y domains 
and domain walls (|Stamper-Kurn et al. 1999). 



Both these modes of coarsening were observed by 
preparing gases initially in the simpler starting condi- 
tion of having just one large domain of each component. 
A magnetic field gradient was then applied, which fa- 
vored the spatially reversed rearrangement of spin do- 
mains along the long-axis of the trap. In this case, the 
chemical potential difference that drives particle motion 
is controlled by the strength of the field gradient, rather 
than by surface tension. 

The quantum tunneling current across the spin domain 
wall is regulated by the height of the potential barrier, 
which is proportional to the spin-dependent mean-field 
energy (and hence to the condensate density), and by its 
width, which is inversely proportional to the magnetic 
field gradient. At high condensate density and weak mag- 
netic field gradients, particle transport was dominated by 
the diffusion and recondensation of the thermal gas. The 
particle current due to this recondensation was slow and 
was relatively insensitive to the field gradient. At lower 
condensate densities and stronger field gradients, trans- 
port was dominated by the direct quantum tunneling of 
the superfluid. The rapid (exponential) variation of the 
tunneling rate led to a rapid rearrangement of the gas 
once the gradually decaying condensate density fell be- 
low a threshold value. 

Upon allowing for spin-mixing among all substates of 
the spinor gas, one expects transport and coarsening 
mechanisms to differ from those studied in the aforemen- 
tioned experiments. For example, in the F = 1 spinor 
gas at low quadratic Zeeman shift, the \m z — 0) popu- 
lation is no longer a locally conserved quantity. Thus, 
domains of \m z = 0) atoms can be transformed by spin 
mixing into the \m z = ±1) states. This transformation 
is seen in quench experiments with both ferromagnetic 



netic ( |Bookjans et al.\ ( |2011b[ ); see Fig. [2| F = 1 conden- 
sates when the \m z = 0) gas is made locally unstable. 
The MIT group studied the fate of \m z = 0) domains 
under conditions where the domains are locally stable 
(due to a positive quadratic Zeeman shift), but globally 
energetically unstable (due to a magnetic field gradient) 
against spin mixing (|Miesner et al. 1999). The \m z = 0) 



domain evaporated via thermally activated spin mixing, 
with the spin mixing products eventually recondensing 
into the \m z = ±1) spin domains that flank the \m z — 0) 
domain on either side. 



C. Long-time equilibration dynamics 

While some of the elemental processes in coarsening 
have been studied, our understanding of the mechanisms 
and timescales for equilibration remains incomplete and 
poses vexing problems for future research. Indeed, re- 
cent experiments on relaxation and coa rsening in spa- 
tially extended F = 1 87 Rb condensates (jGuzman et al. 



2011 Vengalattore et al. 2010) raise the question of 



whether spinor Bose gases can be reliably equilibrated 
within experimentally accessible timescales. These ex- 
periments were motivated by observations made on de- 
generate spinor gases that had undergone dynamical in- 
stabilities, brought about either by rapid sweeps of the 
quadratic Zeeman shift (Leslie et al. 2009b Sadler et al. 
2006) or by the imposition of a helical spin texture (Ven- 



galattore et al. t 2008 ) , and allowed to evolve for several 



times the dynamical instability timescale. Under a vari- 
ety of initial conditions, such gases developed long-lasting 
spatial modulations of the spin texture, comprised typi- 
cally of ~ 10 fim long domains of alternating magnetiza- 
tion. Following the realization that the MDDI should sig- 
nificantly influence the structure of low-energy spin tex- 
tures in rubidium spinor condensates, the question arose 
whether these spatial modulations could be characteristic 
of the thermal equilibrium state. 

This question was explored in two works in which un- 
magnetized, nondegenerate spinor gases were produced, 
cooled gradually by evaporative cooling to below the de- 
generacy temperature, and then held at constant tem- 
perature for a variable equilibration time before being 
probed either by time-of-flight or in situ imaging. The 
first study (Vengalattore et al. 2010), which was lim- 



ited to moderate equilibration times of several 100 ms, 
explored the textures produced at variable temperature 
and for compositions of the nondegenerate gas that dif- 
fered in the fractional population of the \m z =0) state. 
For atoms prepared in an initial (1/3,1/3,1/3) popula- 
tion mixture of the m z = (+1,0, —1) Zeeman substates, 
the fine-scale magnetization modulations were observed 
at all temperatures below the degeneracy temperature. 
In contrast, for atoms prepared in the (1/4, 1/2, 1/4) ini- 
tial state, the modulated spin textures were observed at 
the lowest temperatures, while at higher temperatures 
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FIG. 22 A two component 23 Na spinor Bose-Einstein condensate is quenched by preparing it in a uniform superposition of the 
\F — l,m z = 0, +1) Zeeman sublevels in a narrow, quasi-one-dimensional optical trap, and at a high quadratic Zeeman shift 
for which spin mixing into the unoccupied \m z = — 1) state is strongly suppressed. Stern-Gerlach imaging is used to determine 
the spin-state distribution at a variable time after the state preparation. The gas undergoes phase separation, initially in 
the narrow radial direction (10 ms) but soon settling into a distribution of pure-spin domains along the long axis of the gas 
(vertical in image), with domain lengths of 40 /im or less. At the low temperatures of this experiment (no non-condensed gas is 
discernible) and without a magnetic-field gradient to induce directed diffusion or quantum tunneling, no significant coarsening 
occurs over the few-second lifetime of the trapped gas. Figure reproduced from |Miesner et al. (1999). 



the degenerate gases were nearly uniformly transversely 
magnetized. 

The persistence of the observed spin textures was ex- 
amined theoretically, primarily by computing either the 
ground-state or the temporal evolution determined by 
the zero-temperature Gross-Pitaevskii equation, account- 
ing for bot h the spin-dependent s -wave interactions and 
the MDDI. |Cherng and Demler| ( |2009D found that uni- 
formly magnetized rubidium condensates were subject to 
a roton-like instability with particular wavevectors and 
spatial anisotropy, suggesting an instability toward the 
formation of a spatially periodic ground-state texture. 
The possibility that the observed structures were long- 
lived spin- vortex lattices was suggested by an imagi nary- 
time calculatio n for the ground-state spin texture ( Zhang 
and Ho 2010), but belied by computations of the real- 
time dynamics ( |Kawaguchi et al , 2010). Equilibrium 
phases with striped magnetization patterns were pre- 
dicted to be generated by the competition of the kinetic- 
energy stiffness of the condensate magnetic order param- 
eter and the MDDI ( |Kjall et a/.||2009| ); however, as with 
the other theory works, the predicted length scale of the 
ground-state modulation was several times larger than 
that observed experimentally. 

The mismatch between experiment and theory was 
partly explained by the fact that these finely modulated 
textures were non-equilibrium states, which coarsen over 
several seconds to form larger-scale magnetic domains 
(Fig. [6|. Studying sp inor condensates prod uced at the 



lowest temperatures, Guzman et al. ( |2011| ) found that 



the equilibration time for condensates with different ini- 
tial spin compositions to reach the same final steady-state 
spin compositions depended strongly on the quadratic 
Zeeman energy: for small |g| < hx 10 Hz, condensates ar- 
rived at a common steady-state spin distribution within 
several seconds, while for larger \q\, their equilibration 
was dramatically slowed (Fig. 23). The equilibration dy- 
namics at small \q\ was marked by four trends: First, 
the Zeeman populations were observed to evolve from 



their initial values to a common final steady state. Sec- 
ond, the degenerate gas became increasingly magnetized 
as measured by the mean value of |F 2 (r)|. Third, the 
spin texture developed a spin-space anisotropy, marked 
by the difference between the longitudinal (|F^(r)|) and 
transverse ((|F^(r)| + |F^(r)|)/2) magnetization where z 
denotes the axis defined by the quadratic Zeeman shift. 
As predicted by mean-field theory, the textures showed 
predominantly transverse magnetization (easy-plane) for 
q > and longitudinal magnetization (easy-axis) for 
q < 0. Fourth, the size of the commonly magnetized 
regions coarsened over time from the I ~ 10 /im domain 
length seen at short equilibration times to a final length 
comparable to the overall size of the trapped samples. All 
four trends occurred on a similar time scale of a few sec- 
onds, although the mechanisms that determine the time 
scales for each of these trends are unknown. 

An ambitious target for future experiments is to re- 
late these observed trends quantitatively to the predic- 
tions for phase-ordering kinetics. The presence of a 
Bose-Einstein condensate and the condition of locally 
conserved magnetization from the rotationally symmet- 
ric s-wave interactions (ignoring the MDDI) qualify the 
F = 1 spinor Bose gas as obeying the "Model F" descrip- 
tion^ coarsening, as defined by |Hohenberg and Halperin 
(1977). This assignment is supported by numerical com- 



parisons between the coarsening dynamics predicted by 
the stochastic time-dependent Landau- Ginzburg and the 
energy- cons erving Gross-Pitaevskii equation ( |Mukerjee 



et 



2007). Such a comparison is valid in the q < re- 



gion, where a ferromagnetic F = 1 condensate should de- 
velop longitudinal magnetization with long-range Ising- 
type order. The growth of the charact e ristic domain 
size l(t) was reported by [Guzman et al. (2011), but it 
is doubtful whether a meaningful comparison to the pre- 
dicted power-law growth (l(t) ex t 1 / 3 ) can be made over 
the limited time- and length-scales accessed in that work. 

For ferromagnetic condensates evolving in two dimen- 
sions with < q < 2|c^n|, one does not expect long- 
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FIG. 23 (Color) Equilibration of condensate Zeeman popula- 
tions following the preparation of an unmagnetized thermal 
gas, and evaporative cooling to below the quantum degener- 
acy temperature. The F = 1 87 Rb gas is held in a quasi two- 
dimensional optical potential, well outside the single spatial- 
mode regime. The average population in the \m z = ±1) 
states, (P+i + P-i)/2, is determined at various evolution 
times t = 200 ms (purple circles), 1 s (red squares), and 3 
s (black diamonds) for initial thermal spin populations of (a) 
(1/3,1/3,1/3) or (b) (1/4,1/2,1/4) in the m z = (1,0,-1) 
states. In (c), the long-time (2 seconds) Zeeman populations 
are shown for the initial spin mixtures of figure (a) (green 
triangles), figure (b) (blue circles), and for the populations 
(0,1,0) (yellow squares). For a narrow range \q\/h < 10 
Hz (shaded), Zeeman populations evolve toward a common 
steady state that is qualitatively in accord with mean-field 
predictions (solid line in c) based on the spin-dependent con- 
tact interaction. Outside that range, initial population differ- 
ences persist throughout the experimentally accessed equili- 
bration times. Figure reproduced from |Guzman et a77| ( |2Q1 1| ) . 



range order to develop at thermal equilibrium. Numerical 
studies suggest that the dynamics of such gases following 
a quantum (or thermal) quench may show a long-lived 
"prethermalized" state with well-characterized correla- 
tions that differ fr om those at thermal equilibrium ( |Bar 
|nett et a/.||2011b| ). 



XI. APPLICATION OF SPINOR GASES TO 
MAGNETOMETRY 

Aside from being objects of fundamental interest, 
spinor Bose gases also lend themselves to the applica- 
tion of magnetic field sensing with high sensitivity and 
high spatial resolution. Atomic Bohr frequencies uj may 
generally be measured using an atomic ensemble by the 
following procedure: First, the ensemble is prepared in a 
common polarized state, e.g. through optical pumping or 
Stern- Gerlach separation. Next, to initiate the measure- 
ment, all atoms are placed in a coherent superposition of 
internal states, with a relative phase <\> between the inter- 
nal states. During the measurement time r, <j) advances 
by an amount <j> = ujt. To conclude the measurement, 
the phase <p is estimated, by a single projective measure- 
ment or several weak measurements. 

Subdividing the ensemble spatially allows one to mea- 
sure the spatial inhomogeneity of the Bohr frequency, 
which may come about from an inhomogeneous environ- 
mental variable, such as the magnetic field. For an un- 
corrected initial state, the measurement within a single 
measurement pixel is limited by atomic shot noise to the 
imprecision 6u(l) = [N(l)Dt}~ 1/2 where N(l) oc l d n is 
the number of atoms within a pixel of length I and di- 
mension n is the relevant (volume, column or line) 
density, D is the measurement duty cycle and t the total 
time allotted to the measurement. 

Several features of spinor Bose-Einstein condensates 
qualify them for use in such spatially resolving sensors. 
The Bose-Einstein condensate is a high-density, low- 
energy medium. The high density reduces the atomic 
shot noise as defined above. The low-energy permits the 
use of shallow optical traps, which are benign to the 
atomic spin and from which spontaneous scattering is 
minimal. The absence of thermal diffusion implies that 
atoms remain well localized for long times. Finally, the 
rotational symmetry of the s-wave interactions ensures 
that Larmor precession of a fully magnetized gas un- 
der a linear Zeeman shift is unaffected by interactions. 
Altogether, these qualities allow for small measurement 
imprecision and systematic error. 

High-resolution magnetometry with two-dimensional 
spatial resolution was demonstrated using an F = 1 87 Rb 



spinor Bose-Einstein condensate ( Vengalattore et al. 
2007). The gas was prepared in a transversely mag- 



netized state and allowed to Larmor precess for r = 
250 ms before the precession phase was read out us- 
ing magnetization-sensitive phase contrast imaging. The 
sensor was tested by measuring the measurement noise 
both in the absence and the presence of a short length- 
scale magnetic signal, which was generated via the vec- 
tor ac Stark shift of a focused, circularly polarized, off- 
resonant light beam (Fig.[24|). The Allan deviation, mea- 
sured for measurement areas A = I 2 in the range 6 - 
400 /im 2 , was limited by photon shot noise at a level 
slightly higher than the atomic shot-noise limit. Even in 
this proof-of-principle experiment, the achieved measure- 
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FIG. 24 (Color) Spatially resolved magnetometry with a 
spinor Bose-Einstein condensate. A transversely magnetized 
condensate undergoes Larmor precession before the preces- 
sion phase is measured by magnetization sensitive imaging, 
(a) Repeated images of one component of the transverse mag- 
netization, from which (b) the Larmor precession phase (indi- 
cated by color) and (c) the magnetic field is spatially resolved 
across the gas. (d) The field sensitivity determined from 
imaging noise (grey line, assuming unity duty cycle) agrees 
well with the photon shot noise limit (PSN) on magnetiza- 
tion imaging, and remains above the atomic shot-noise limit 
(ASN). The sensitivity determined also from measurements 
of the light-induced magnetic field (shown in (a-c)) is plot- 
ted for the actual (filled circle) or for unity duty cycle (open 
circle). Also shown are the quantum noise limit for SQUID 
magnetometers (dot dashed), and the sensitivities achieved in 
dc SQUID magnetometers (triangles) of Kirtley et al. ( 1995 ) 
and |Lee et al] ( |1996| . Dotted line indicates the worsening 
of sensitivity due to atomic motion, an effect that can be 
avoided by use of an op tical lattice. Figure reproduced from 
Vengalattore et aU (12007). 



merit sensitivity was comparable to that achieved with 
low-noise SQUID-based magnetic microscopes at similar 
spatial resolution 6 . 



Hot-atom magnetometers achieve much better field sensitivity 
( |Budker and Rom alis 2007), but are limited in spatial resolution 
at the millimeter scale 



Magnetic field sensing has also been demonstrated us- 
ing scalar Bose-Einstein condensates, based on the obser- 
vation that spatial variations of the magnetic field on the 
order of 1 nT produce variations in the magnetic poten- 
tial energy, which then lead to visible density variations 



in the compressible condensed gas (Wildermuth et al. 



2005). However, the detection of Larmor precession al- 
lowed for measurements of far smaller field variations 
that were in the 1 pT range. Spatially resolved Larmor 
precession in a cold, optically trapped, non-condensed 
gas was recently demonstrated, but was limited to poorer 
spatial resolution (I 2 > 1000 /mi 2 ) and short m easure- 



ment times (r < 1ms) ( |Koschorreck et al > 2011). 

Several measurement targets for such a cold-atom sen- 
sor have been suggested. The aforementioned scalar con- 
densate magnetometer has been used to image the mag- 
netic field corrugations generated by current flow through 
a disordered thin conductor ( |Aigner et al] |2008[ ) . This 
measurement established that cold atoms can be used as 
a direct probe of solid-state science. Cold atomic gases 
have also been used as a "scanning-probe microscope" to 
visualize three-dimensional micro- and nano-structures 
(at micron-s cale resolution) ( |Gierling et al. 2011 Hunger 



et al. 



2010). Mechanically vibrating structures that are 
either ferromagnetic or carry current will produce time- 
varying magnetic fields that can induce spin flip tran- 
sitions in nearby spinor gases, allowing for detection of 
small scale motion and, perhaps, for coherent coupling 
between the two systems ( Kalma n" et al\ |2011| |Steinke| 

|2007| ). Alternately, sens- 



ed 



2011 Treutlein et al 



ing temporal fluctuations in the Larmor precession fre- 
quency may provide a means to sense, and, indeed, to 
image vortices migrating through thin superconducting 
films (IScheel et al] [20071. 



Several extant or foreseeable technological advances 
may make such spinor-based magnetometry more prac- 
tical and powerful. The duty cycle of the sensor, which 
was just D = 0.003 in the work of Vengalattore et al. 



( 2007 ) , can be increased by shortening the time required 
to form the spinor condensate (to around 1 s) and by 
increasing the measurement time; indeed, the Berkeley 
group has recently observed coherent Larmor precession 
in a spinor condensate lasting over 3 s. The appara- 
tus required for producing large spinor condensates can 
be shrunk and simplified significantly, e.g. as demon- 
strated by |Salim et al] ( |2011[ ). Thin membranes can 
be employed to separate the spinor gas from its mea- 
surement target, similar to methods used to probe room- 
temperature objects using cryogenic-temperature SQUID 
microscopes ( |Lee et al] [1996), so that the measurement 
target need not be UHV compatible and can be rapidly 
exchanged. Additional improvements may stem from uti- 
lizing the physical properties of the spinor condensates 
themselves. For example, spin-mixing instabilities can 



be used to am plify weak spatially varying signals (Leslie 
et al, 2009b), allowing higher sensitivity measurement 



even at low optical-detection efficiency. Spatially multi- 
mode spin-nematic squeezing, generated either by atomic 
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interactions or by measurement, may provide measure- 
ment resolution below the atomic shot-noise limits (Leslie 



et 



2009b Sau et 



2010) 



XII. CONCLUSIONS 

Several concepts regarding spinor Bose gases were 
emphasized in this Review. First, we have empha- 
sized the properties of spinor gases that are necessitated 
by rotational symmetry. In contrast to generic multi- 
component quantum fluids, which can be represented 
formally in terms of a pseudo-spin internal-state wave- 
function, spinor gases are comprised of particles whose 
internal-state wavefunction responds under geometric ro- 
tations as a spherical tensor of rank 2F. If it were not for 
experimental conditions that break symmetry, the Hamil- 
tonian of a self-interacting spinor gas would necessarily 
be invariant under rotations. 

In Sec. |III[ we have discussed how this rotational in- 
variance constrains the interactions among spinor gas 
atoms. If we first neglect magnetic dipolar interactions, 
the ultracold nature of spinor gases indicates the use of 
partial waves to describe interparticle scattering, with 
only the lowest partial wave contributing to low-energy 
collisions. By rotational symmetry, such s-wave binary 
collisions are distinguished only by the total spin of the 
colliding pair, not by its orientation. By Bose symmetry, 
s-wave collisions occur only when this total spin is even. 
Altogether, the interaction Hamiltonian is enormously 
simplified by these symmetries. 

The symmetries of this interaction Hamiltonian are 
reflected also in its mean-field stationary solutions and 
its exact many-body ground states. In Sec. [TTJ we pre- 
sented a scheme to visualize these solutions graphically. 
Rotational invariance of the Hamiltonian implies that 
these stationary states reside in a degenerate manifold 
of states, and we outlined consequences of this fact in 
describing some of the unusual textures and topological 
defec ts tha t occur in spinor Bose-Einstein condensates 
(Sec. VI. D[ ) and in treating their low-energy dynamics 
Finally, the rotational symmetry of inter- 



VI.E 



(Sec. _ 

actions implies that, for certain internal states, the dy- 
namics of interacting spinor Bose gases under applied 
magnetic fields are the same as those of free particles. 
On the one hand, in the absence of dipolar relaxation, 
this fact allows one to gauge away the effects of uniform 
magnetic fields and study low-energy dynamics of certain 
spinor gases even in the presence of substa ntial and un- 
controlled magnetic fields (discussed in Sec. IV.B.l). On 
the other hand, the immunity of Larmor precession to 
interaction effects makes th e sp inor Bose gas suitable for 
magnetic field sensing (Sec. XI). 

A second concept central to this Review is the inter- 
play between magnetic order and superfluidity. Both 
phenomena signify phase coherence among atoms in the 
fluid: magnetic order stems from the phases controlling 
internal-state coherences, while superfluidity, at least in 



three dimensions, stems from coherence of the global 
phase of the gas-particle wavefunctions. The long-range 
phase coherence of a superfluid implies that a spinor 
Bose gas display also magnetic order, though the con- 
verse implication does not necessarily hold. The question 
of whether Bose-Einstein condensation and magnetic or- 
dering must occur simultaneously was taken up in Sec. 



VII where we considered three scenarios in which long- 



range coherence and the influence of Bose statistics are 
weakened, namely non-zero temperature, restriction of 
the spinor gas to lower dimensions, and the effects of 
optical lattice potentials. 

The connection between superfluid and magnetic order 
was discussed also in the context of spin textures in Sec. 



VI We have highlighted the relation between superflow 



and the spin-rotation symmetry of the ground-state man- 
ifolds of a spinor Bose-Einstein condensate (Sec. VLB), 
a relation familiar from studies of superfluid 3 He that is 
manifest in the properties of the varied textures and de- 
fects predicted for spinor condensates. This relation also 
arose in the description of low-energy hydrodynamics of 



spinor Bose-condensates, as discussed in Sec. VI.E 



A third focus of this Review was non-equilibrium quan- 
tum dynamics. Quantum degenerate atomic gases are a 
new class of materials, and many of the scientific concepts 
explored in research on such gases are similar to those 
explored with condensed matter. On the one hand, the 
short lifetime of the gas, its isolation from a thermal envi- 
ronment, and the presence of conserved quantities make 
it challenging to associate the observed behavior of the 
gas with a state of thermal equilibrium. On the other 
hand, these same properties make the gaseous systems 
particularly appealing for studies of non-equilibrium dy- 
namical effects that would be hard to simulate in denser 
solids and liquids. Focusing on non-equilibrium spin dy- 
namics is particularly appealing because the spin state of 
the gas can be initialized with very high fidelity, at noise 
levels much smaller than might be expected at thermal 
equilibrium given the nK-scale kinetic temperatures of 
the gas. 



In Sec. IX we discussed the non-equilibrium phenom- 
ena of coherent spin-mixing oscillations and instabilities. 
We highlighted the connection between few-body, many- 
body, and mean-field treatments of such dynamics, il- 
lustrated by experimental examples. In discussing the 
rich area of spin-mixing instabilities, we drew connec- 
tions to four- wave mixing of optical and matter waves, 
the dynamics of systems near quantum phase transitions, 
symmetry-breaking in closed quantum systems, and the 
dynamical generation of correlated quantum spin states 
of practical use for metrology. 

With these far-reaching developments, the authors of 
this review consider the field of spinor Bose gases still 
in its infancy. We conclude this review by enumerating 
some outstanding issues. One important goal is to make 
quantitative test of mean-field and Bogoliubov theories 
on which so many theoretical studies are based. This 
can be done most straightforwardly by measuring the 
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collective modes. An equally pressing issue is to un- 
derstand why non-equilibrium states persist so long in 
spinor gases. This persistence marks a sharp contrast 
with scalar gases which relax to their equilibrium states 
quickly, and also raises the question of whether it is ex- 
perimentally feasible to explore equilibrium properties of 
spinor Bose gases, and, if so, under what circumstances. 
The question of equilibration also hinges on the proper- 
ties of spinor Bose gases at non-zero temperature, which 
remain poorly understood and unexplored. 

Another inviting research target is the characterization 
of correlated spin states produced either dynamically or 
at equilibrium in spinor Bose-Einstein condensates. Re 
cent studies on spin-nematicity squeezing (|Hamley et al. 



2012 Liicke et al., 2011 ) have opened this line of investi- 



gation. The beautiful prediction of a fragmented ground 
state for antiferromagnetic spinor condensates by Law 
et al.\ {1998) may be confirmed by many of the same 
techniques used in recent experiments. 

Topological aspects are also worthwhile to investigate 
because spinor condensates are among the rare systems 
for which the needed manipulations to prepare topolog- 
ical objects are experimentally available. The evolution 
of topological structures in spin textures, and their role 
in determining the character of a spinor Bose gas in re- 
stricted dimensions, will offer important examples and 
fresh perspectives on the dynamics of quantum fields. 
New topological structures may also be found at inter- 
faces between different phases of the spinor Bose gas; for 
example, it would be interesting to ask whether such in- 
terfaces could accommodate topological defects such as 
a boojum, which is a point singularity on a s urface in 
superfluid helium-3 (iBhattacharyya et al, 1977). 



The topic of spinor gases in optical lattices is espe- 
cially rich. Such lattices add new twists to the physics of 
spinor gases because the dimensions, filling fraction, and 
geometry of the lattices can introduce strong correlations, 
phase transitions, and spin frustration in the same sys- 
tem. Together with newly developed techniques for syn- 
thesizing gauge fields and for generating more complex, 
spin-dependent, lattice potentials, such lattice-trapped 
spinor gases offer a new paradigm for exploring quantum 
phenomena and simulations. 



Acknowledgments 

The authors are indebted to our close colleagues within 
our research groups with whom we have explored the 
physics of spinor Bose gases. At Berkeley, we acknowl- 
edge the experimental contributions of A. Chikkatur, J. 
Guzman, J. Higbie, S. Inouye, G.-B. Jo, S. Leslie, K. 
Murch, L. Sadler, V. Savalli, F. Serwane, C. Thomas, M. 
Vengalattore, A. Wenz, and also the theoretical inputs of 
J. Sau and M. Cohen. M.U. acknowledges his scientific 
collaborators including M. Kobayashi, Y. Kawaguchi, M. 
Koashi, M. Nitta, H. Saito, and S. Uchino. We are also 
grateful to A. Lamacraft, W.V. Liu, J. Moore and E. 



Mueller for valuable discussions during the preparation 
of this work, and to C. Raman and Y. Liu for assistance 
with several figures. D.S.-K. thanks the NSF, the ARO 
with funds from the DARPA OLE program, the Division 
of Materials Sciences and Engineering at the Office of Ba- 
sic Energy Sciences (DOE), the Sloan and Packard Foun- 
dations, and the Hellman Faculty Fund and the Miller In- 
stitute at the University of California, for their support. 
M.U. acknowledge support by Grants-in-Aid for Scien- 
tific Research (Grant Nos. 22103005 and 22340114), a 
Global COE Program "the Physical Sciences Frontier," 
and the Photon Frontier Network Program, from MEXT 
of Japan. 



References 

Aigner, S., L. D. Pietra, Y. Japha, O. Entin-Wohlman, 

T. David, R. Salem, R. Folman, and J. Schmiedmayer 

(2008), Science 319, 1226. 
Aikawa, K., A. Frisch, M. Mark, S. Baier, A. Rietzler, 

R. Grimm, and F. Ferlaino (2012), Phys. Rev. Lett. 108, 

210401. 

Al Khawaja, U., and H. Stoof (2001), Nature 411, 918. 
Andrews, M., M.-O. Mewes, N. van Druten, D. Durfee, 

D. Kurn, and W. Ketterle (1996), Science 273, 84. 
Ao, P., and S. Chui (1998), Phys. Rev. A 58, 4836. 
Ao, P., and S. Chui (2000), Journal of Physics B-Atomic, 

Molecular and Optical Physics 33, 535. 
Apaja, V., and O. F. Syljuasen (2006), Phys. Rev. A 74, 

035601. 

Appel, J., P. J. Windpassinger, D. Oblak, U. B. HofT, 
N. Kjaergaard, and E. S. Polzik (2009), Proc. Natl. Acad. 
Sci. USA 106, 10960. 

Arndt, M., M. Ben Dahan, D. Guery-Odelin, M. Reynolds, 
and J. Dalibard (1997), Phys. Rev. Lett. 79, 625. 

Arodz, H., J. Dziarmaga, and W. H. Zurek (2003), Patterns 
of symmetry breaking, NATO science series Series II, Math- 
ematics, physics, and chemistry (Kluwer Academic Pub- 
lishers, Dordrecht ; Boston). 

Baraban, M., H. F. Song, S. M. Girvin, and L. I. Glazman 
(2008), Phys. Rev. A 78, 033609. 

Baranov, M., L. Dobrek, K. Goral, L. Santos, and M. Lewen- 
stein (2002), Physica Scripta T102, 74. 

Baranov, M. A. (2008), Physics Reports- Review Section of 
Physics Letters 464, 71. 

Barnett, R., H.-Y. Hui, C.-H. Lin, J. D. Sau, and 
S. Das Sarma (2011a), Phys. Rev. A 83, 023613. 

Barnett, R., S. Mukerjee, and J. E. Moore (2008), Phys. Rev. 
Lett. 100, 240405. 

Barnett, R., D. Podolsky, and G. Refael (2009), Phys. Rev. 
B 80, 024420. 

Barnett, R., A. Polkovnikov, and M. Vengalattore (2011b), 

Phys. Rev. A 84, 023606. 
Barnett, R., J. D. Sau, and S. Das Sarma (2010), Phys. Rev. 

A 82, 031602. 

Barnett, R., A. Turner, and E. Dernier (2006), Phys. Rev. 

Lett. 97, 180412. 
Barnett, R., A. Turner, and E. Dernier (2007), Phys. Rev. A 

76, 013605. 
Barnett, S. J. (1915), Phys. Rev. 6, 239. 



55 



Barrett, M. D., J. A. Sauer, and M. S. Chapman (2001), 

Phys. Rev. Lett. 87, 010404. 
Barton, G., and M. A. Moore (1974), J Phys C Solid State 

7, 4220. 

Barton, G., and M. A. Moore (1975), J Phys C Solid State 

8, 970. 

Batrouni, G. G., V. G. Rousseau, and R. T. Scalettar (2009), 

Phys. Rev. Lett. 102, 140402. 
Bauerle, C., Y. M. Bunkov, S. N. Fisher, H. Godfrin, and 

G. R. Pickett (1996), Nature 382, 332. 
Berglund, A. J., S. A. Lee, and J. J. McClelland (2007), Phys. 

Rev. A 76, 053418. 
Bernier, J. S., K. Sengupta, and Y. B. Kim (2007), Phys. 

Rev. B 76, 014502. 
Bhattacharyya, P., T. L. Ho, and N. D. Mermin (1977), Phys. 

Rev. Lett. 39, 1290. 
Bigelow, N. P., J. H. Freed, and D. M. Lee (1989), Phys. Rev. 

Lett. 63, 1609. 

Black, A. T., E. Gomez, L. D. Turner, S. Jung, and P. D. 

Lett (2007), Phys. Rev. Lett. 99, 070403. 
Blaha, S. (1976), Phys. Rev. Lett. 36, 874. 
Boesten, H. M. J. M., A. J. Moerdijk, and B. J. Verhaar 

(1996), Phys. Rev. A 54, R29. 
Bogoliubov, N. (1947), J. Phys. (USSR) 11, 23. 
Bookjans, E. M., C. D. Hamley, and M. S. Chapman (2011a), 

Phys. Rev. Lett. 107, 210406. 
Bookjans, E. M., A. Vinit, and C. Raman (2011b), Phys. 

Rev. Lett. 107, 195306. 
Bowick, M. J., L. Chandar, E. A. SchifT, and A. M. Srivastava 

(1994), Science 263, 943. 
Brahms, N., T. P. Purdy, D. W. C. Brooks, T. Botter, and 

D. M. Stamper-Kurn (2011), Nature Physics 7, 604. 
Bravyi, S., M. B. Hastings, and F. Verstraete (2006), Phys. 

Rev. Lett. 97, 050401. 
Bray, A. J. (2002), Adv Phys 51, 481 . 
Breit, G., and I. I. Rabi (1931), Phys. Rev. 38, 2082. 
Budker, D., W. Gawlik, D. F. Kimball, S. M. Rochester, V. V. 

Yashchuk, and A. Weis (2002), Rev. Mod. Phys. 74, 1153. 
Budker, D., and M. Romalis (2007), Nature Physics 3, 227. 
Burke, J., Jr., J. Bohn, B. Esry, and C. Greene (1997), Phys. 

Rev. A 55, R2511. 
Calabrese, P., and J. Cardy (2006), Phys. Rev. Lett. 96, 

136801. 

Cao, J. P., Y. Z. Jiang, and Y. P. Wang (2007), Epl 79. 
Caramico D'Auria, A., L. De Cesare, and I. Rabuffo (1996), 

Physica A: Statistical Mechanics and its Applications 225, 

363. 

Carmi, R., E. Polturak, and G. Koren (2000), Phys. Rev. 
Lett. 84, 4966. 

Carusotto, I., and E. J. Mueller (2004), J. Phys. B 37, S115. 
Cazalilla, M. A., R. Citro, T. Giamarchi, E. Orignac, and 

M. Rigol (2011), Rev. Mod. Phys. 83, 1405. 
Chang, L., Q. Zhai, R. Lu, and L. You (2007), Phys. Rev. 

Lett. 99, 080402. 
Chang, M.-S., C. D. Hamley, M. D. Barrett, J. A. Sauer, 

K. M. Fortier, W. Zhang, L. You, and M. S. Chapman 

(2004) , Phys. Rev. Lett. 92, 140403. 

Chang, M.-S., Q. Qin, W. Zhang, L. You, and M. Chapman 

(2005) , Nature Physics 1, 111. 

Chen, A. X., W. Y. Qiu, and Z. P. Wang (2008), Chinese 

Physics B 17, 4204. 
Chen, P. C, Z. L. Xue, I. P. McCulloch, M. C. Chung, and 

S. K. Yip (2012), Phys. Rev. A 85, 011601. 
Cheng, J. (2010), Phys. Rev. A 81, 023619. 



Cherng, R. W., and E. Dernier (2009), Phys. Rev. Lett. 103, 
185301. 

Cherng, R. W., V. Gritsev, D. M. Stamper-Kurn, and 
E. Dernier (2008), Phys. Rev. Lett. 100, 180404. 

Choi, J. Y., W. J. Kwon, and Y. I. Shin (2012), Phys. Rev. 
Lett. 108, 035301. 

Chuang, I., R. Durrer, N. Turok, and B. Yurke (1991), Sci- 
ence 251, 1336. 

Ciobanu, C, S.-K. Yip, and T.-L. Ho (2000), Phys. Rev. A 
61, 033607. 

Cohen-Tannoudji, C, and J. Dupont-Roc (1972), Phys. Rev. 
A 5, 968. 

Cramer, M., A. Flesch, I. P. McCulloch, U. Schollwck, and 

J. Eisert (2008), Phys. Rev. Lett. 101, 063001. 
Crubellier, A., O. Dulieu, F. Masnou-Seeuws, M. Elbs, 

H. Knockel, and E. Tiemann (1999), European Physical 

Journal D 6, 211. 
Cui, X. L., Y. P. Wang, and F. Zhou (2008), Phys. Rev. A 

78, 050701. 

Damski, B., and W. H. Zurek (2007), Phys. Rev. Lett. 99, 
130402. 

De Chiara, G., M. Lewenstein, and A. Sanpera (2011), Phys. 

Rev. B 84, 054451. 
Dernier, E., and F. Zhou (2002), Phys. Rev. Lett. 88, 163001. 
Deng, L., E. W. Hagley, J. Wen, M. Trippenbach, Y. Band, 

P. S. Julienne, J. E. Simsarian, K. Helmerson, S. L. Rolston, 

and W. D. Phillips (1999), Nature 398, 218. 
Deutsch, C, F. Ramirez-Martinez, C. Lacroute, F. Reinhard, 

T. Schneider, J. N. Fuchs, F. Piechon, F. Laloe, J. Reichel, 

and P. Rosenbusch (2010), Phys. Rev. Lett. 105, 020401. 
Diener, R. B., and T.-L. Ho (2006), "Quantum spin dynamics 

of spin-1 bose gas," Preprint, arXiv:cond-mat/0608732. 
Donley, E., N. Claussen, S. Cornish, J. Roberts, E. Cornell, 

and C. Wieman (2001), Nature 412, 295. 
Du, X., Y. Zhang, J. Petricka, and J. E. Thomas (2009), 

Phys. Rev. Lett. 103, 010401. 
Duan, L. M., J. I. Cirac, and P. Zoller (2002), Phys. Rev. A 

65. 

Duan, L.-M., A. Sorensen, J. Cirac, and P. Zoller (2000), 

Phys. Rev. Lett. 85, 3991 . 
Dziarmaga, J. (2005), Phys. Rev. Lett. 95, 245701. 
Dziarmaga, J. (2010), Adv Phys 59, 1063. 
Eckert, K., L. Zawitkowski, M. J. Leskinen, A. Sanpera, and 

M. Lewenstein (2007a), New Journal of Physics 9, 133. 
Eckert, K., L. Zawitkowski, A. Sanpera, M. Lewenstein, and 

E. S. Polzik (2007b), Phys. Rev. Lett. 98, 100404. 
Einstein, A., and W. J. de Haas (1916), P K Akad Wet- 

Amsterd 18, 696. 
Eisenberg, E., and E. H. Lieb (2002), Phys. Rev. Lett. 89, 

220403. 

Endo, Y., and T. Nikuni (2011), Journal of Low Temperature 

Physics 163, 92. 
Essler, F. H. L., G. V. Shlyapnikov, and A. M. Tsvelik (2009), 

Journal of Statistical Mechanics- Theory and Experiment . 
Est eve, J., C. Gross, A. Weller, S. Giovanazzi, and M. K. 

Oberthaler (2008), Nature 455, 1216. 
Falke, S., H. Knoeckel, J. Friebe, M. Riedmann, and E. Tie- 
mann (2008), Phys. Rev. A 78, 012503. 
Fano, U., and J. H. Macek (1973), Rev. Mod. Phys. 45, 553. 
Fattori, M., G. Roati, B. Deissler, C. D'Errico, M. Zaccanti, 

M. Jona-Lasinio, L. Santos, M. Inguscio, and G. Modugno 

(2008), Phys. Rev. Lett. 101, 190405. 
Fisher, M. P. A., P. B. Weichman, G. Grinstein, and D. S. 

Fisher (1989), Phys. Rev. B 40, 546. 



56 



Frota, H. O. d., M. S. Silva, and J. S. G. Rosa (1984), J. 

Phys. C 17, 1669. 
Fuchs, J. N., D. M. Gangardt, and F. Laloe (2002), Phys. 

Rev. Lett. 88, 230404. 
Gawryluk, K., M. Brewczyk, K. Bongs, and M. Gajda (2007), 

Phys. Rev. Lett. 99, 130401. 
de Gennes, P. G. (1971), Symposia of the Faraday Society 5, 

16. 

Gerbier, F., A. Widera, S. Foiling, O. Mandel, and I. Bloch 

(2006), Phys. Rev. A 73, 041602. 
Gerlach, W., and O. Stern (1924), Ann. Phys.-Berlin 74, 673. 
Gerton, J. M., C. A. Sackett, B. J. Frew, and R. G. Hulet 

(1999), Phys. Rev. A 59, 1514. 
Gierling, M., P. Schneeweiss, G. Visanescu, P. Federsel, 

M. Haeffner, D. P. Kern, T. E. Judd, A. Guenther, and 

J. Fortagh (2011), Nature Nanotechnology 6, 446. 
Giorda, P., P. Zanardi, and S. Lloyd (2003), Phys. Rev. A 

68, 062320. 

Goldstein, E., and P. Meystre (1997), Phys. Rev. A 55, 2935. 
Gorlitz, A., T. L. Gustavson, A. E. Leanhardt, R. Low, A. P. 

Chikkatur, S. Gupta, S. Inouye, D. E. Pritchard, and 

W. Ketterle (2003), Phys. Rev. Lett. 90, 090401. 
Graham, R., and D. Walls (1998), Phys. Rev. A 57, 484. 
Greiner, M., O. Mandel, T. Esslinger, T. W. Hansen, and 

I. Bloch (2002), Nature 415, 39. 
Gross, C., H. Strobel, E. Nicklas, T. Zibold, N. Bar-Gill, 

G. Kurizki, and M. K. Oberthaler (2011), Nature 480, 

219. 

Gross, C., T. Zibold, E. Nicklas, J. Esteve, and M. K. 

Oberthaler (2010), Nature 464, 1165. 
Gross, K., C. P. Search, H. Pu, W. P. Zhang, and P. Meystre 

(2002), Phys. Rev. A 66, 033603. 
Gu, Q. (2003), Phys. Rev. A 68, 025601. 
Gu, Q., K. Bongs, and K. Sengstock (2004), Phys. Rev. A 

70. 

Gu, Q., and R. Klemm (2003), Phys. Rev. A 68, 031604. 
Guery-Odelin, D., J. Soding, P. Desbiolles, and J. Dalibard 

(1998), Europhysics Letters 44, 25. 
Guzman, J., G. B. Jo, A. N. Wenz, K. W. Murch, C. K. 

Thomas, and D. M. Stamper-Kurn (2011), Phys. Rev. A 

84, 063625. 

Hadzibabic, Z., and J. Dalibard (2011), Rivista Del Nuovo 

Cimento 34, 389. 
Hall, D., M. Matthews, J. Ensher, C. Wieman, and E. Cornell 

(1998a), Phys. Rev. Lett. 81, 4531. 
Hall, D., M. Matthews, C. Wieman, and E. Cornell (1998b), 

Phys. Rev. Lett. 81, 1543. 
Hamley, C. D., C. S. Gerving, T. M. Hoang, E. M. Bookjans, 

and M. S. Chapman (2012), Nat Phys 8, 305. 
Hammerer, K., A. S. Sorensen, and E. S. Polzik (2010), Rev. 

Mod. Phys. 82, 1041. 
Harber, D. M., H. Lewandowski, J. McGuirk, and E. Cornell 

(2002), Phys. Rev. A 66, 053616. 
Hecht, E. (1989), Optics, 2nd ed. (Addison- Wesley, Reading). 
Hendry, P. C, N. S. Lawson, R. A. M. Lee, P. V. E. McClin- 

tock, and C. D. H. Williams (1994), Nature 368, 315. 
Hensler, S., J. Werner, A. Griesmaier, P. O. Schmidt, A. Gor- 
litz, T. Pfau, S. Giovanazzi, and K. Rzazewski (2003), 

App. Phys. B 77, 765. 
Higbie, J., L. Sadler, S. Inouye, A. P. Chikkatur, S. R. Leslie, 

K. L. Moore, V. Savalli, and D. M. Stamper-Kurn (2005), 

Phys. Rev. Lett. 95, 050401. 
Ho, T.-L. (1998), Phys. Rev. Lett. 81, 742. 
Ho, T.-L., and V. Shenoy (1996), Phys. Rev. Lett. 77, 2595. 



Ho, T.-L., and S. Yip (2000), Phys. Rev. Lett. 84, 4031. 
Ho, T. L., and Q. Zhou (2010), Nature Physics 6, 66. 
Hohenberg, P. C, and B. I. Halperin (1977), Rev. Mod. Phys. 
49, 435. 

Horikoshi, M., S. Nakajima, M. Ueda, and T. Mukaiyama 

(2010), Science 327, 442. 
Hou, J. M., and M. L. Ge (2003), Phys. Rev. A 67, 063607. 
Huang, K., and C. N. Yang (1957), Phys. Rev. 105, 767. 
Huhtamaki, J. A. M., M. Takahashi, T. P. Simula, 

T. Mizushima, and K. Machida (2010), Phys. Rev. A 81, 

063623. 

Hunger, D., S. Camerer, T. W. Hnsch, D. Knig, J. P. Kot- 
thaus, J. Reichel, and P. Treutlein (2010), Phys. Rev. Lett. 
104, 143002. 

Imambekov, A., M. Lukin, and E. Dernier (2003), Phys. Rev. 
A 68, 063602. 

Imambekov, A., M. Lukin, and E. Demler (2004), Phys. Rev. 

Lett. 93, 120405. 
Isacsson, A., and S. M. Girvin (2005), Phys. Rev. A 72, 

53604. 

Isoshima, T., and K. Machida (2002), Phys. Rev. A 66, 
023602. 

Isoshima, T., K. Machida, and T. Ohmi (1999), Phys. Rev. 
A 60, 4857. 

Isoshima, T., K. Machida, and T. Ohmi (2001), J. Phys. Soc. 
Jpn. 70, 1604. 

Isoshima, T., T. Ohmi, and K. Machida (2000), J. Phys. Soc. 
Jpn. 69, 3864. 

Jae-yoon, C, K. Woo Jin, L. Moonjoo, J. Hyunseok, 

A. Kyungwon, and S. Yong-il (2012), New Journal of 

Physics 14, 053013. 
Jaksch, D., C. Bruder, J. Cirac, C. Gardiner, and P. Zoller 

(1998), Phys. Rev. Lett. 81, 3108. 
James, A. J. A., and A. Lamacraft (2011), Phys. Rev. Lett. 

106, 140402. 

Jensen, K., V. M. Acosta, J. M. Higbie, M. P. Ledbetter, 
S. M. Rochester, and D. Budker (2009), Phys. Rev. A 79, 
023406. 

Jing, H., Y. Jiang, and P. Meystre (2009), Phys. Rev. A 80, 
063618. 

Johnson, B. R., J. S. Denker, N. Bigelow, L. P. Levy, J. H. 

Freed, and D. M. Lee (1984), Phys. Rev. Lett. 52, 1508. 
Kalman, O., T. Kiss, J. Fortagh, and P. Domokos (2011), 

Nano Letters 12, 435. 
Kasamatsu, K., M. Tsubota, and M. Ueda (2005), Interna- 
tional Journal of Modern Physics B 19, 1835. 
Kawaguchi, Y., M. Nitta, and M. Ueda (2008), Phys. Rev. 

Lett. 100, 180403. 
Kawaguchi, Y., H. Saito, K. Kudo, and M. Ueda (2010), 

Phys. Rev. A 82, 043627. 
Kawaguchi, Y., H. Saito, and M. Ueda (2006a), Phys. Rev. 

Lett. 96, 080405. 
Kawaguchi, Y., H. Saito, and M. Ueda (2006b), Phys. Rev. 

Lett. 97, 130404. 
Kawaguchi, Y., H. Saito, and M. Ueda (2007), Phys. Rev. 

Lett. 98, 110406. 
Kawaguchi, Y., and M. Ueda (2011), Phys. Rev. A 84, 

053616. 

Kawaguchi, Y., and M. Ueda (2012), Physics Reports 520, 
253. 

van Kempen, E., S. Kokkelmans, D. Heinzen, and B. Verhaar 

(2002), Phys. Rev. Lett. 88, 093201. 
Ketterle, W., D. Durfee, and D. Stamper-Kurn (1999), 

"Making, probing and understanding Bose-Einstein con- 



57 



densates," in Bose-Einstein condensation in atomic gases, 
Proceedings of the International School of Physics "Enrico 
Fermi," Course CXL, edited by M. Inguscio, S. Stringari, 
and C. Wieman (IOS Press, Amsterdam) pp. 67-176. 

Kibble, T. W. B. (1976), J. Phys. A 9, 1387. 

Kimura, T., S. Tsuchiya, and S. Kurihara (2005), Phys. Rev. 
Lett. 94, 110403. 

Kirtley, J. R., M. B. Ketchen, K. G. Stawiasz, J. Z. Sun, 
W. J. Gallagher, S. H. Blanton, and S. J. Wind (1995), 
App. Phys. Lett. 66, 1138. 

Kis-Szabo, K., P. Szepfalusy, and G. Szirmai (2005), Phys. 
Rev. A 72. 

Kis-Szabo, K., P. Szepfalusy, and G. Szirmai (2007), Phys. 

Lett. A 364, 362. 
Kita, T., T. Mizushima, and K. Machida (2002), Phys. Rev. 

A 66, 061601. 

Kitagawa, M., and M. Ueda (1993), Phys. Rev. A 47, 5138. 
Kjall, J. A., A. M. Essin, and J. E. Moore (2009), Phys. Rev. 
B 80, 224502. 

Klausen, N., J. Bohn, and C. Greene (2001), Phys. Rev. A 
64, 053602. 

Kleine Biining, G., J. Will, W. Ertmer, E. Rasel, J. Arlt, 
C. Klempt, F. Ramirez-Martinez, F. Piechon, and 
P. Rosenbusch (2011), Phys. Rev. Lett. 106, 240801. 

Klempt, C., O. Topic, G. Gebreyesus, M. Scherer, T. Hen- 
ninger, P. Hyllus, W. Ertmer, L. Santos, and J. J. Arlt 

(2009) , Phys. Rev. Lett. 103, 195302. 

Klempt, C., O. Topic, G. Gebreyesus, M. Scherer, T. Hen- 
ninger, P. Hyllus, W. Ertmer, L. Santos, and J. J. Arlt 

(2010) , Phys. Rev. Lett. 104, 195303. 

Koashi, M., and M. Ueda (2000), Phys. Rev. Lett. 84, 1066. 
Kobayashi, M., Y. Kawaguchi, M. Nitta, and M. Ueda (2009), 

Phys. Rev. Lett. 103, 115301. 
Kobayashi, S., M. Kobayashi, Y. Kawaguchi, M. Nitta, and 

M. Ueda (2012), Nuclear Physics B 856, 577. 
Kokkelmans, S. J. J. M. F., B. J. Verhaar, and K. Gibble 

(1998), Phys. Rev. Lett. 81, 951. 
Korshunov, S. (1985), Sov. Phys. JETP 62, 301. 
Koschorreck, M., M. Napolitano, B. Dubost, and M. W. 

Mitchell (2010), Phys. Rev. Lett. 104, 093602. 
Koschorreck, M., M. Napolitano, B. Dubost, and M. W. 

Mitchell (2011), App. Phys. Lett. 98. 
Kronjager, J., C. Becker, M. Brinkmann, R. Walser, P. Navez, 

K. Bongs, and K. Sengstock (2005), Phys. Rev. A 72, 

063619. 

Kronjager, J., C. Becker, P. Navez, K. Bongs, and K. Seng- 
stock (2006), Phys. Rev. Lett. 97, 110404. 

Kronjager, J., C. Becker, P. Soltan-Panahi, K. Bongs, and 
K. Sengstock (2010), Phys. Rev. Lett. 105, 090402. 

Kronjager, J., K. Sengstock, and K. Bongs (2008), New Jour- 
nal of Physics 10, 045028. 

Krutitsky, K. V., M. Timmer, and R. Graham (2005), Phys. 
Rev. A 71, 033623. 

Kudo, K., and Y. Kawaguchi (2010), Phys. Rev. A 82, 
053614. 

Kuwamoto, T., K. Araki, T. Eno, and T. Hirano (2004), 
Phys. Rev. A 69, 063604. 

Kuzmich, A., N. P. Bigelow, and L. Mandel (1998), Euro- 
physics Letters 42, 481. 

Kuzmich, A., L. Mandel, and N. Bigelow (2000), Phys. Rev. 
Lett. 85, 1594. 

Lahaye, T., T. Koch, B. Frohlich, M. Fattori, J. Metz, 
A. Griesmaier, S. Giovanazzi, and T. Pfau (2007), Nature 
448, 672. 



Lahaye, T., C. Menotti, L. Santos, M. Lewenstein, and 
T. Pfau (2009), Reports on Progress in Physics 72, 126401. 

Lahaye, T., J. Metz, B. Froehlich, T. Koch, M. Meister, 
A. Griesmaier, T. Pfau, H. Saito, Y. Kawaguchi, and 
M. Ueda (2008), Phys. Rev. Lett. 101, 080401. 

Lamacraft, A. (2007), Phys. Rev. Lett. 98, 160404. 

Lamacraft, A. (2008), Phys. Rev. A 77, 063622. 

Lamacraft, A. (2010), Phys. Rev. B 81, 184526. 

Landau, L., and E. Lifshitz (1977), Quantum Mechanics: 
Non-Relativistic Theory (Pergamon Press, New York). 

Landau, L. D., E. M. Lifshits, and L. P. Pitaevski'i (1984), 
Electrodynamics of continuous media, 2nd ed., Course of 
theoretical physics / L D Landau and E M Lifshitz (Perg- 
amon, Oxford Oxfordshire ; New York). 

Lauchli, A., G. Schmid, and S. Trebst (2006), Phys. Rev. B 
74, 144426. 

Law, C, H. Pu, and N. Bigelow (1998), Phys. Rev. Lett. 81, 
5257. 

Lacki, M., S. Paganelli, V. Ahufinger, A. Sanpera, and J. Za- 

krzewski (2011), Phys. Rev. A 83, 013605. 
Leanhardt, A. E., A. Grlitz, A. P. Chikkatur, D. Kielpinski, 

Y. Shin, D. E. Pritchard, and W. Ketterle (2002), Phys. 

Rev. Lett. 89, 190403. 
Lee, J. Y., X. W. Guan, M. T. Batchelor, and C. Lee (2009), 

Phys. Rev. A 80. 
Lee, T. S., E. Dantsker, and J. Clarke (1996), Rev. Sci. In- 

strum. 67, 4208. 
Leonhardt, U., and G. E. Volovik (2000), Jetp Lett+ 72, 46. 
Leslie, L. S., A. Hansen, K. C. Wright, B. M. Deutsch, and 

N. P. Bigelow (2009a), Phys. Rev. Lett. 103, 250401. 
Leslie, S., J. Guzman, M. Vengalattore, J. D. Sau, M. L. 

Cohen, and D. Stamper-Kurn (2009b), Phys. Rev. A 79, 

043631. 

Lewandowski, H. J., D. M. Harber, D. L. Whitaker, and E. A. 

Cornell (2002), Phys. Rev. Lett. 88, 070403. 
Lewenstein, M., A. Sanpera, V. Ahufinger, B. Damski, A. Sen, 

and U. Sen (2007), Adv Phys 56, 243. 
Lewenstein, M., and L. You (1996), Phys. Rev. Lett. 77, 

3489. 

Lhuiller, C, and F. Laloe (1982a), J. Physique 43, 197. 
Lhuiller, C, and F. Laloe (1982b), J. Physique 43, 225. 
Li, Y. Q., S. J. Gu, and Z. J. Ying (2003), Journal of Physics 

a-Mathematical and General 36, 2821. 
Li, Z., A. X. Zhang, J. Ma, and J. K. Xue (2010), Phys. Lett. 

A 374, 1401. 

Li, Z. D., P. B. He, L. Li, J. Q. Liang, and W. M. Liu (2005), 

Phys. Rev. A 71, 053611. 
Lin, Y. J., R. L. Compton, K. Jimenez-Garcia, W. D. Phillips, 

J. V. Porto, and I. B. Spielman (2011), Nat Phys 7, 531. 
Lin, Y. J., R. L. Compton, K. Jimenez-Garcia, J. V. Porto, 

and I. B. Spielman (2009), Nature 462, 628. 
Liu, W. V., and C. Wu (2006), Phys. Rev. A 74, 013607. 
Liu, Y., E. Gomez, S. E. Maxwell, L. D. Turner, E. Tiesinga, 

and P. D. Lett (2009a), Phys. Rev. Lett. 102, 225301. 
Liu, Y., S. Jung, S. E. Maxwell, L. D. Turner, E. Tiesinga, 

and P. D. Lett (2009b), Phys. Rev. Lett. 102, 125301. 
Lu, M., N. Q. Burdick, S. H. Youn, and B. L. Lev (2011), 

Phys. Rev. Lett. 107, 190401. 
Liicke, B., M. Scherer, J. Kruse, L. Pezze, F. Deuretzbacher, 

P. Hyllus, O. Topic, J. Peise, W. Ertmer, J. Arlt, L. Santos, 

A. Smerzi, and C. Klempt (2011), Science 334, 773. 
Lysebo, M., and L. Veseth (2010), Phys. Rev. A 81, 032702. 
Majorana, E. (1932), Nuovo Cimento 9, 43. 
Makela, H., and K. A. Suominen (2007), Phys. Rev. Lett. 



58 



99, 190408. 

Makela, H., Y. Zhang, and K. A. Suominen (2003), Journal 
of Physics a-Mathematical and General 36, 8555. 

Marcin, F., and S. Krzysztof (2012), Journal of Physics A: 
Mathematical and Theoretical 45, 045103. 

Matthews, M., D. Hall, D. Jin, J. Ensher, C. Wieman, E. Cor- 
nell, F. Dalfovo, C. Minniti, and S. Stringari (1998), Phys. 
Rev. Lett. 81, 243. 

Matuszewski, M. (2010), Phys. Rev. Lett. 105, 020405. 

Matuszewski, M., T. J. Alexander, and Y. S. Kivshar (2008), 
Phys. Rev. A 78, 023632. 

McClelland, J. J., and J. L. Hanssen (2006), Phys. Rev. Lett. 
96, 143005. 

McGuirk, J., H. Lewandowski, D. Harber, T. Nikuni, 
J. Williams, and E. Cornell (2002), Phys. Rev. Lett. 89, 
090402. 

McGuirk, J. M., D. M. Harber, H. J. Lewandowski, and E. A. 

Cornell (2003), Phys. Rev. Lett. 91, 150402. 
Medley, P., D. M. Weld, H. Miyake, D. E. Pritchard, and 

W. Ketterle (2011), Phys. Rev. Lett. 106, 195301. 
Mermin, N. D. (1979), Rev. Mod. Phys. 51, 591. 
Mermin, N. D., and T.-L. Ho (1976), Phys. Rev. Lett. 36, 

594. 

Metcalf, H. J., and P. v. d. Straten (1999), Laser cooling and 

trapping (Springer- Verlag, New York). 
Mias, G. I., N. R. Cooper, and S. M. Girvin (2008), Phys. 

Rev. A 77, 023616. 
Michel, L. (1980), Rev. Mod. Phys. 52, 617. 
Miesner, H.-J., D. Stamper-Kurn, J. Stenger, S. Inouye, 

A. Chikktur, and W. Ketterle (1999), Phys. Rev. Lett. 

82, 2228. 

Mizushima, T., N. Kobayashi, and K. Machida (2004), Phys. 
Rev. A 70. 

Mizushima, T., K. Machida, and T. Kita (2002a), Phys. Rev. 
A 66, 053610. 

Mizushima, T., K. Machida, and T. Kita (2002b), Phys. Rev. 

Lett. 89, 030401. 
Monaco, R., J. Mygind, and R. J. Rivers (2002), Phys. Rev. 

Lett. 89, 080603. 
Mueller, E. J. (2004), Phys. Rev. A 69, 033606. 
Mueller, E. J., T. L. Ho, M. Ueda, and G. Baym (2006), 

Phys. Rev. A 74, 033612. 
Mukerjee, S., C. Xu, and J. E. Moore (2006), Phys. Rev. 

Lett. 97, 120406. 
Mukerjee, S., C. Xu, and J. E. Moore (2007), Phys. Rev. B 

76. 

Miiller, T., S. Foiling, A. Widera, and I. Bloch (2007), Phys. 

Rev. Lett. 99, 200405. 
Mur-Petit, J., M. Guilleumas, A. Polls, A. Sanpera, 

M. Lewenstein, K. Bongs, and K. Sengstock (2006), Phys. 

Rev. A 73, 13629. 
Myatt, C, E. Burt, R. Ghrist, E. Cornell, and C. Wieman 

(1997), Phys. Rev. Lett. 78, 586. 
Nakahara, M., T. Isoshima, K. Machida, S.-i. Ogawa, and 

T. Ohmi (2000), Physica B: Condensed Matter 284, 17. 
Nakamura, K. (1994), Quantum Chaos: A New Paradigm of 

Nonlinear Dynamics, Cambridge Nonlinear Science (Cam- 
bridge University Press, Cambridge). 
Natu, S. S., and E. J. Mueller (2010), Phys. Rev. A 81, 

053617. 

Natu, S. S., and E. J. Mueller (2011), Phys. Rev. A 84, 
053625. 

Navon, N., S. Nascimbene, F. Chevy, and C. Salomon (2010), 
Science 328, 729. 



Nikuni, T., and Y. Endo (2008), Journal of Low Temperature 

Physics 152, 21. 
Nozieres, P., and D. Saint James (1982), J Phys-Paris 43, 

1133. 

Ohmi, T., and K. Machida (1998), J. Phys. Soc. Jpn. 67, 
1822. 

Oktel, M. ., and L. S. Levitov (2002), Phys. Rev. Lett. 88, 
230403. 

Partridge, G. B., J. C. Jaskula, M. Bonneau, D. Boiron, and 

C. I. Westbrook (2010), Phys. Rev. A 81, 053631. 
Pasquiou, B., G. Bismut, Q. Beauflls, A. Crubellier, 

E. Marechal, P. Pedri, L. Vernac, O. Gorceix, and 

B. Laburthe-Tolra (2010), Phys. Rev. A 81, 042716. 
Pasquiou, B., E. Marchal, L. Vernac, O. Gorceix, and 

B. Laburthe-Tolra (2012), Phys. Rev. Lett. 108, 045307. 
Pasquiou, B., E. Marechal, G. Bismut, P. Pedri, L. Vernac, 

O. Gorceix, and B. Laburthe-Tolra (2011), Phys. Rev. 

Lett. 106, 255303. 
Penrose, O., and L. Onsager (1956), Phys. Rev. 104, 576. 
Phuc, N. T., Y. Kawaguchi, and M. Ueda (2011), Phys. Rev. 

A 84, 043645. 

Pietila, V., T. P. Simula, and M. Mottonen (2010), Phys. 

Rev. A 81, 033616. 
Pitaevskii, L. (1997), Phys. Lett. A 229, 406. 
Podolsky, D., S. Chandrasekharan, and A. Vishwanath 

(2009), Phys. Rev. B 80, 214513. 
Polkovnikov, A. (2005), Phys. Rev. B 72, 161201. 
Polkovnikov, A., K. Sengupta, A. Silva, and M. Vengalattore 

(2011), Rev. Mod. Phys. 83, 863. 
Pollack, S. E., D. Dries, M. Junker, Y. P. Chen, T. A. Cor- 

covilos, and R. G. Hulet (2009), Phys. Rev. Lett. 102, 

090402. 

Pu, H., C. Law, S. Raghavan, J. Eberly, and N. Bigelow 

(1999), Phys. Rev. A 60, 1463. 
Pu, H., and P. Meystre (2000), Phys. Rev. Lett. 85, 3987. 
Pu, H., S. Raghavan, and N. Bigelow (2000), Phys. Rev. A 

61, 023602. 

Pu, H., W. P. Zhang, and P. Meystre (2001), Phys. Rev. 

Lett. 87, 140405. 
Rajaraman, R. (1987), An Introduction to Solitons and In- 

stantons in Quantum Field Theory, Solitons and Instan- 

tons, Vol. 15 (North Holland Personal Library, Elsevier 

Science, Amsterdam). 
Riedel, M. F., P. Bohi, Y. Li, T. W. Hansch, A. Sinatra, and 

P. Treutlein (2010), Nature 464, 1170. 
Rizzi, M., D. Rossini, G. De Chiara, S. Montangero, and 

R. Fazio (2005), Phys. Rev. Lett. 95, 240404. 
Robins, N. P., W. Zhang, E. A. Ostrovskaya, and Y. S. 

Kivshar (2001), Phys. Rev. A 64, 021601. 
Rodriguez, K., A. Arguelles, A. K. Kolezhuk, L. Santos, and 

T. Vekua (2011), Phys. Rev. Lett. 106, 105302. 
Romano, D. R., and E. J. V. de Passos (2004), Phys. Rev. A 

70, 043614. 

Rong, C, J. Q. Liang, and Z. Yunbo (2005), Journal of 
Physics B: Atomic, Molecular and Optical Physics 38, 
2569. 

Ruostekoski, J., and J. R. Anglin (2003), Phys. Rev. Lett. 
91, 190402. 

Ruutu, V. M. H., V. B. Eltsov, A. J. Gill, T. W. B. Kibble, 
M. Krusius, Y. G. Makhlin, B. Placais, G. E. Volovik, and 
W. Xu (1996), Nature 382, 334. 

Sackett, C, H. Stoof, and R. Hulet (1998), Phys. Rev. Lett. 
80, 2031. 

Sadler, L., J. Higbie, S. Leslie, M. Vengalattore, and 



59 



D. Stamper-Kurn (2006), Nature 443, 312. 
Saito, H., Y. Kawaguchi, and M. Ueda (2006), Phys. Rev. 

Lett. 96, 065302. 
Saito, H., Y. Kawaguchi, and M. Ueda (2007a), Phys. Rev. 

A 76, 043613. 

Saito, H., Y. Kawaguchi, and M. Ueda (2007b), Phys. Rev. 
A 75, 013621. 

Saito, H., and M. Ueda (2005), Phys. Rev. A 72, 053628. 

Salim, E., J. DeNatale, D. Farkas, K. Hudek, S. McBride, 
J. Michalchuk, R. Mihailovich, and D. Anderson (2011), 
Quantum Information Processing 10, 975. 

Santos, L., M. Fattori, J. Stuhler, and T. Pfau (2007), Phys. 
Rev. A 75, 053606. 

Santos, L., and T. Pfau (2006), Phys. Rev. Lett. 96, 190404. 

Sau, J., S. Leslie, M. Cohen, and D. Stamper-Kurn (2009), 
Phys. Rev. A 80, 023622. 

Sau, J. D., S. R. Leslie, M. D. Cohen, and D. M. Stamper- 
Kurn (2010), New Journal of Physics 12, 085011. 

Scheel, S., R. Fermani, and E. A. Hinds (2007), Phys. Rev. 
A 75. 

Scherer, M., B. Lucke, G. Gebreyesus, O. Topic, F. Deuret- 
zbacher, W. Ertmer, L. Santos, J. J. Arlt, and C. Klempt 
(2010), Phys. Rev. Lett. 105, 135302. 

Schleier-Smith, M. H., I. D. Leroux, and V. Vuletic (2010), 
Phys. Rev. Lett. 104, 073604. 

Schmaljohann, H., M. Erhard, J. Kronjager, M. Kottke, S. v. 
Staa, L. Cacciapuoti, J. J. Arlt, K. Bongs, and K. Seng- 
stock (2004), Phys. Rev. Lett. 92, 040402. 

Semenoff, G. W., and F. Zhou (2007), Phys. Rev. Lett. 98, 
100401. 

Shankar, R. (1977), J Phys-Paris 38, 1405. 

Shin, Y. I., C. H. Schunck, A. Schirotzek, and W. Ketterle 

(2008), Nature 451, 689. 
Siggia, E., and E. Ruckenstein (1980), Phys. Rev. Lett. 44, 

1423. 

Simkin, M. V., and E. G. D. Cohen (1999), Phys. Rev. A 59, 
1528. 

Smith, G. A., S. Chaudhury, A. Silberfarb, I. H. Deutsch, and 
P. S. Jessen (2004), Phys. Rev. Lett. 93, 163602. 

Snoek, M., J. L. Song, and F. Zhou (2009), Phys. Rev. A 80, 
053618. 

Soding, J., D. Guery-Odelin, P. Desbiolles, G. Ferrari, and 
J. Dalibard (1998), Phys. Rev. Lett. 80, 1869. 

Song, J. L., G. W. Semenoff, and F. Zhou (2007), Phys. Rev. 
Lett. 98, 160408. 

Sorensen, A., L.-M. Duan, J. Cirac, and P. Zoller (2001), 
Nature 409, 63. 

Stamper-Kurn, D., and W. Ketterle (2001), "Spinor conden- 
sates and light scattering from Bose-Einstein condensates," 
in Coherent Matter Waves, edited by R. Kaiser, C. West- 
brook, and F. David, Chap. 2 (Springer- Verlag, New York) 
pp. 137 - 218. 

Stamper-Kurn, D., H.-J. Miesner, A. Chikkatur, S. Inouye, 
J. Stenger, and W. Ketterle (1999), Phys. Rev. Lett. 83, 
661. 

Stamper-Kurn, D. M., M. Andrews, A. Chikkatur, S. Inouye, 

H.-J. Miesner, J. Stenger, and W. Ketterle (1998), Phys. 

Rev. Lett. 80, 2072. 
Steinke, S. K., S. Singh, M. E. Tasgin, P. Meystre, K. C. 

Schwab, and M. Vengalattore (2011), Phys. Rev. A 84, 

023841. 

Stenger, J., S. Inouye, D. Stamper-Kurn, H.-J. Miesner, 

A. Chikkatur, and W. Ketterle (1998), Nature 396, 345. 
Stoof, H. T. C, E. Vliegen, and U. Al Khawaja (2001), Phys. 



Rev. Lett. 87, 120407. 

Sugawa, S., K. Inaba, S. Taie, R. Yamazaki, M. Yamashita, 
and Y. Takahashi (2011), Nature Physics 7, 642. 

Sukachev, D., A. Sokolov, K. Chebakov, A. Akimov, 
S. Kanorsky, N. Kolachevsky, and V. Sorokin (2010a), 
Phys. Rev. A 82, 011405. 

Sukachev, D., A. Sokolov, K. Chebakov, A. Akimov, N. Ko- 
lachevsky, and V. Sorokin (2010b), Jetp Lett+ 92, 703. 

Sun, B., W. X. Zhang, S. Yi, M. S. Chapman, and L. You 

(2006) , Phys. Rev. Lett. 97, 123201. 

Suter, D. (2005), The Physics of Laser- Atom Interactions, 
Cambridge Studies in Modern Optics (Cambridge Univer- 
sity Press, Cambridge). 

Takahashi, M., S. Ghosh, T. Mizushima, and K. Machida 

(2007) , Phys. Rev. Lett. 98, 260403. 

Takano, T., M. Fuyama, R. Namiki, and Y. Takahashi (2009), 

Phys. Rev. Lett. 102, 033601. 
Tao, C. J., P. L. Wang, J. H. Qin, and Q. Gu (2008), Phys. 

Rev. B 78. 

Timmermans, E. (1998), Phys. Rev. Lett. 81, 5718. 
Treutlein, P., P. Hommelhoff, T. Steinmetz, T. W. Hansen, 

and J. Reichel (2004), Phys. Rev. Lett. 92, 203005. 
Treutlein, P., D. Hunger, S. Camerer, T. W. Hansch, and 

J. Reichel (2007), Phys. Rev. Lett. 99, 140403. 
Tsuchiya, S., S. Kurihara, and T. Kimura (2004), Phys. Rev. 

A 70, 043628. 

Turner, A. M., R. Barnett, E. Demler, and A. Vishwanath 

(2007), Phys. Rev. Lett. 98, 190404. 
Turner, A. M., and E. Demler (2009), Phys. Rev. B 79, 

214522. 

Uchino, S., M. Kobayashi, and M. Ueda (2010), Phys. Rev. 
A 81, 063632. 

Ueda, M. (2012), Annual Review of Condensed Matter 
Physics 3, 263. 

Ueda, M., and M. Koashi (2002), Phys. Rev. A 65, 063602. 
Uhlmann, M., R. Schutzhold, and U. R. Fischer (2007), Phys. 

Rev. Lett. 99, 120407. 
Vengalattore, M., J. Guzman, S. R. Leslie, F. Serwane, and 

D. M. Stamper-Kurn (2010), Phys. Rev. A 81, 053612. 
Vengalattore, M., J. M. Higbie, S. R. Leslie, J. Guzman, L. E. 

Sadler, and D. M. Stamper-Kurn (2007), Phys. Rev. Lett. 

98, 200801. 

Vengalattore, M., S. Leslie, J. Guzman, and D. Stamper- 
Kurn (2008), Phys. Rev. Lett. 100, 170403. 

Vitagliano, G., P. Hyllus, I. n. L. Egusquiza, and G. Toth 
(2011), Phys. Rev. Lett. 107, 240502. 

Vollhardt, D., and P. Wolfle (1990), The Superfluid Phases 
of Helium 3 (Taylor and Francis, New York). 

Volovik, G. E. (1987), J Phys C Solid State 20, L83. 

Volovik, G. E. (2003), The Universe in a Helium Droplet, In- 
ternational Series of Monographs on Physics (Oxford Uni- 
versity Press, Oxford). 

Volovik, G. E., and V. P. Mineyev (1977), Zhurnal Eksperi- 
mentalnoi I Teoreticheskoi Fiziki 73, 767. 

Wagner, A., C. Bruder, and E. Demler (2011), Phys. Rev. A 
84, 063636. 

Weinstein, J. D., R. deCarvalho, J. Kim, D. Patterson, 
B. Friedrich, and J. M. Doyle (1998), Phys. Rev. A 57, 
R3173. 

Werner, J., A. Griesmaier, S. Hensler, J. Stuhler, T. Pfau, 
A. Simoni, and E. Tiesinga (2005), Phys. Rev. Lett. 94, 
183201. 

Wheeler, M. H., K. M. Mertes, J. D. Erwin, and D. S. Hall 
(2004), Phys. Rev. Lett. 93, 170402. 



60 



Widera, A., F. Gerbier, S. Foiling, T. Gericke, O. Mandel, 

and I. Bloch (2005), Phys. Rev. Lett. 95, 190405. 
Widera, A., F. Gerbier, S. Foiling, T. Gericke, O. Mandel, 

and I. Bloch (2006), New Journal of Physics 8, 152. 
Wildermuth, S., S. HofTerberth, I. Lesanovsky, E. Haller, 

L. M. Andersson, S. Groth, I. Bar- Joseph, P. Kruger, and 

J. Schmiedmayer (2005), Nature 435, 440. 
Williams, J. E., T. Nikuni, and C. W. Clark (2002), Phys. 

Rev. Lett. 88, 230405. 
Wineland, D., J. Bollinger, W. Itano, F. Moore, and 

D. Heinzen (1992), Phys. Rev. A 46, R6797. 
Wirth, G., M. Olschlager, and A. Hemmerich (2011), Nat 

Phys 7, 147. 

Wright, E., D. Walls, and J. Garrison (1996), Phys. Rev. 
Lett. 77, 2158. 

Wright, K. C., L. S. Leslie, and N. P. Bigelow (2008), Phys. 
Rev. A 78, 053412. 

Wright, K. C., L. S. Leslie, A. Hansen, and N. P. Bigelow 
(2009), Phys. Rev. Lett. 102, 030405. 

Xie, Z.-W., W. Zhang, S. T. Chui, and W. M. Liu (2004), 
Phys. Rev. A 69, 053609. 

Yamada, K. (1982), Progress of Theoretical Physics 67, 443. 

Yang, K., and Y. Q. Li (2003), International Journal of Mod- 
ern Physics B 17, 1027. 

Yasunaga, M., and M. Tsubota (2008), Phys. Rev. Lett. 101, 
220401. 

Yasunaga, M., and M. Tsubota (2010a), Phys. Rev. A 81, 
023624. 

Yasunaga, M., and M. Tsubota (2010b), Journal of Low Tem- 
perature Physics 158, 51. 

Yasunaga, M., and M. Tsubota (2011), Phys. Rev. A 83, 
013618. 

Yi, S., and H. Pu (2006), Phys. Rev. Lett. 97, 020401. 
Yi, S., L. You, and H. Pu (2004), Phys. Rev. Lett. 93, 040403. 
Yip, S.-K. (1999), Phys. Rev. Lett. 83, 4677. 
Yip, S. K. (2003), Phys. Rev. Lett. 90, 250402. 



Yip, S. K. (2007), Phys. Rev. A 75, 023625. 

Youn, S. H., M. W. Lu, U. Ray, and B. L. Lev (2010), Phys. 
Rev. A 82, 043425. 

Yukawa, E., and M. Ueda (2012), "Hydrodynamic de- 
scription of spin-1 bose-einstein condensates," Preprint, 
arXiv:1209.5509. 

Zhang, J. A., and T. L. Ho (2010), Journal of Low Temper- 
ature Physics 161, 325. 

Zhang, J. M., S. Cui, H. Jing, D. L. Zhou, and W. M. Liu 

(2009) , Phys. Rev. A 80, 043623. 

Zhang, K., L. Zhou, H. Y. Ling, H. Pu, and W. Zhang (2011), 

Phys. Rev. A 83, 063624. 
Zhang, P., C. K. Chan, X. G. Li, Q. K. Xue, and X. G. Zhao 

(2002), Phys. Rev. A 66, 043606. 
Zhang, P., A. Z. Zhang, S. Q. Duan, and X. G. Zhao (2001), 

Phys. Lett. A 286, 107. 
Zhang, W., D. L. Zhou, M. S. Chang, M. S. Chapman, and 

L. You (2005), Phys. Rev. A 72, 13602. 
Zhang, W. X., B. Sun, M. S. Chapman, and L. You (2010), 

Phys. Rev. A 81, 033602. 
Zhang, W. X., S. Yi, and L. You (2004), Phys. Rev. A 70, 

043611. 

Zhou, F. (2001), Phys. Rev. Lett. 87, 080401. 
Zhou, F., and G. W. Semenoff (2006), Phys. Rev. Lett. 97, 
180411. 

Zhou, F., and M. Snoek (2003), Annals of Physics 308, 692. 
Zhou, L., H. Pu, H. Y. Ling, K. Y. Zhang, and W. P. Zhang 

(2010) , Phys. Rev. A 81, 063641. 

Zhou, L., H. Pu, H. Y. Ling, and W. P. Zhang (2009), Phys. 

Rev. Lett. 103, 160403. 
Zurek, W. H. (1985), Nature 317, 505. 

Zurek, W. H. (1996), Physics Reports- Review Section of 

Physics Letters 276, 177. 
Zurek, W. H., U. Dorner, and P. Zoller (2005), Phys. Rev. 

Lett. 95, 105701. 



